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Supersonic Flow About Slender Bodies of 
Elliptic Cross Section 


A. KAHANE* ann A. SOLARSKI* 
Republic Aviation Corporation 


SUMMARY 


A theory for slender bodies of elliptic cross section in super- 
sonic flow, based on the generalized cross-section slender body 
theory of Ward,’ is presented. The theory is applied to the zero 
angle of attack pressure distribution on elliptic cones and para- 
bolic arc bodies of elliptic cross section; results of several ex 
ample computations are presented. Comparison is made with 
the nonlinear elliptic cone calculations and experimental results 
of Ferri;®> good agreement is obtained. The wave drag of bodies 
with elliptic bases is also treated. Application of the theory to 


bodies at angle of attack is indicated. 


INTRODUCTION 


— BODY OF REVOLUTION in supersonic flow has 
been treated extensively in the literature (see Van 
Dyke! for bibliography), since the appearance of the 
initial paper on the subject by von K Arman and Moore.’ 
Bodies of noncircular cross section, on the other hand, 


have received but brief attention. 


Maslen* has presented a linearized solution method 
for cones of arbitrary cross section in which a finite 
number of oblique supersonic line sources extending 
downstream from the cone vertex are positioned so as 
to satisfy the boundary conditions for the body. A 
linearized conical flow treatment of the elliptic cone, 
which yields an infinite series solution, was given by 
Laporte and Bartels.‘ Graham? has treated with lin 
ear theory, using a continuous distribution of quadri- 
poles along the body axis, a body that undergoes a 
transition from a circular cross section to a two-lobed 
cross section. A (the so-called 
linearized method of characteristics) was originated by 

Presented at the Aerodynamics Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953 

* Principal Aerodynamicist, Development Division 

t Senior Aerodynamicist, Development Division; now with 
Cook Research Laboratories 


nonlinear method, 


Ferri’ for treatment of bodies of arbitrary shape and 
cross section. This method should yield more accurate 
results than any of the linearized methods but, on the 
other hand, requires the lengthiest calculations. The 
method entails, first, a calculation by the nonlinear 
method of characteristics of the flow field about the 
equivalent body of revolution (same area distribution 
as given body); a new solution is then obtained by per- 
turbing the boundary Finally, Ward’ 
has presented a general first-order linear solution for 
slender bodies of arbitrary cross section, the application 
of which involves the determination of the two-dimen- 
sional incompressible potential for the expanding cross- 
sectional boundary of the body. 


conditions. 


This paper treats the slender body of elliptic cross 
section, moving at supersonic velocity, according to the 
method of Ward. The theory is applicable to bodies 
with either pointed noses or infinite cylindrical fore- 
bodies and of arbitrary shape except as subject to re- 
strictions on slope and curvature as detailed in the foot- 
note on page 514. Solutions of particular cases are 
given in relatively simple analytical closed forms, which 
are more amenable to computation than those of any 
of the existing methods. The theory is applied to 
cones, and comparisons are made with the experimental 
and linearized characteristics method results of Ferri.° 
A further example calculated is that of a parabolic arc 
body, and the resulting pressure distribution is com- 
pared with that for the equivalent parabolic body of 
revolution. The wave drag of bodies having elliptic 
bases is considered, and wave drag coefficients for ellip- 
tic cones are calculated. Application of the method 
to bodies at angle of attack is also indicated. 

t Consider a plane fixed in space normal to the flight direction 
The trace of a body passing through the plane comprises a 
boundary, which, in general, expands (or contracts) and trans 


lates. 
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A comparison of the results of several theories for 
cones and parabolic bodies of revolution is given in 
Appendix (B) to indicate the order of accuracy of the 
slender theory. 


SYMBOLS 


a = ellipse semimajor axis 

b = ellipse semiminor axis 

B = Mach Number function, B? = M,? — 1 

c = one-half distance between foci of ellipses, c? = 
a? — b? 

Gn, Da, Cn 

d,, An, Bn? = series coefficients, n = 0, 1, 2,... 

Cu Do 

Cp = drag coefficient based on base area 

Cp = static pressure coefficients, Cp = (p — p:)/q 

d = effective thickness ratio (thickness ratio of 
equivalent body of revolution ) 

oe” = maximum thickness of body of unit length 

D = drag force 

F(c/r, 0) = function defined by Eq. (28) 


G(c/r, @) function defined by Eq. (A-1) 
M = Mach Number 


p = static pressure 

q = dynamic pressure 

r,0,s = polar coordinates 

R = radius of circular body boundary 
= axial coordinate 

S = cross-sectional area 

l = time 

Uz, Uy, Us = perturbation velocity components 

U = free-stream velocity 

a() = differential of ( ) along outward normal to 

boundary 

€ = ratio of ellipse minor and major axes, e = b/a 

¥ = ratio of specific heats 

m = circular cone semivertex angle 

tn = elliptic coordinates 

p = air mass density 

o = dummy variable for s 

T = tangent to boundary 

v = outward normal to boundary 

77) = potential function 

Subscripts 

c = circular cross section 

e = elliptic cross section 

0 = boundary value 

1 = free-stream value 


Primes denote differentiation with respect to s. 


THE WARD SLENDER Bopy THEORY 


Consider a body in supersonic flight at velocity U and 
fix a cartesian (x, y, s) system of coordinates, with 
origin on the body nose, so that the s-axis is in the 
stream direction and increasing toward the rear of the 
body. On assumption of vanishing fluid viscosity and 
heat conductivity, a potential g exists, and a second- 
order nonlinear hyperbolic equation for g as dependent 
variable can be deduced from the equations of motion. 
The perturbation velocity components u,, u,, u, are 
expressible in terms of the potential by 
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_O0¢ _O¢ , 
Us = [ > uw = l : mie f ~ (1) 
Ox oy Os 


The nonlinear potential equation is linearized by neg- 
lecting squares and products of the velocity perturba- 
tions, giving® 


B’¢,, — Grr ~ Ow = 0 (2) 


where B? = M,? — 1 and M, is the flight Mach Num- 
ber. Ward,’ using operational methods, obtains a 
general series solution of Eq. (2) in polar coordinates 
and, after further approximation of this solution in the 
neighborhood of the surface of the body, obtains the 
so-called slender approximation to ¢, 


S’(s) - 
g= Inr + bo + > a,r-" cos (nO + B,) (3) 
a7 1 


where the coefficients a, and bp are functions of s; by is 
given by 


l B ; 
bos) = = S"(s) a — — f In (s — o)S’’(e) do (4 
24 2 ; 


It will be noted from Eq. (3) that, in any plane s = 
constant, ¢ is a harmonic function 


i.e., a solution of 
the two-dimensional Laplace equation, 


Prr + Sy = 0 (5) 


which governs the two-dimensional flow of an inviscid 
incompressible fluid. It is apparent that Eq. (5) also 
can be obtained directly from Eq. (2) by assuming 
that the flight Mach Number is close to unity or ¢,, is 
small compared with ¢,, and ¢,,. The latter assump 
tion is equivalent to the assumption that the body is 
slender, since it has been shown’ that near the body 
¢ss is O(d*? In d*), whereas ¢,, and ¢,, are O(1). From 
Eq. (5) the solution for the flow about such a body may 
be obtained by treatment of the relevant two-dimen- 
sional incompressible flow problem. This approach 
was used in 1924 by Munk’ to obtain the lift and mom- 
ent of slender airship hulls, by Jones’? with regard to 
the lift of thin pointed low aspect-ratio wings, and 
more recently by many other investigators. The 
solution of Ward is an improvement over the method 
of Munk and Jones (whose “‘cross flow’’ approach is 
applicable only to systems at angle of attack), since it 
furnishes the relation for bo, which is required for zero 


angle of attack problems—i.e., thickness problems. 


t The following restrictions apply to the Ward? slender body 
If the body is of unit length and has a maximum thick 
ness of d*, then, since body is assumed slender, d* must be small 


theory. 
compared to unity. The angle that any tangent plane to the 
boundary makes with the free-stream direction must be small 
and 0(d*); 
0(d*) 
radius of curvature of the body in any plane normal to the free- 


the rate of change of this angle must be small and 
If d** is the maximum thickness of any section, then the 


stream direction must be 0(1/d**) for all points where the section 
is convex outward, with no restriction if the curve is convex in 
wards. With these restrictions, the error in Eq. (3) is a factor 
1 + O(d*) 
remainder of the text are from Ward.’ 


This error term and the error terms given in the 
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SLENDER BODIES OF 


PoTENTIAL FOR EXPANDING TWwo-DIMENSIONAL 
ELLIPTIC BOUNDARY 


Consider a slender body of elliptic cross section at 
zero angle of attack moving at supersonic speed. As 
the body passes through a plane fixed in space normal 
to the flight direction, the locus of the intersection of 
the body boundary and the plane comprises an ex- 
panding (or contracting) elliptic curve (see Fig. 1). 
Since the partial differential equation of the flow is the 
two-dimensional Laplace equation, Eq. (5), which 
governs inviscid, incompressible, and steady or un- 
steady flows, the problem is reduced to the determina- 
tion of the potential of an elliptic curve (which may be 
considered the boundary of a two-dimensional elastic 
body) which expands in a fixed plane. Since the fluid 
surrounding the elliptic boundary is considered in- 
compressible, lateral disturbances in the plane are prop- 
agated at infinite velocity, and the flow pattern at any 
instant is dependent, therefore, only on the boundary 
conditions (the normal velocity distribution on the 
boundary) at that instant. This type of 
boundary value problem is known in the mathematical 
literature as the Neumann problem or the second bound- 
The problem 


classical 


ary value problem of potential theory. 
may be stated as follows: Let the closed elliptic 
boundary be denoted by C and the region of fluid ex- 
terior to the boundary by R. The partial differential 


equation and boundary conditions are 


Grr + $y = 0,in R 


¢, = ¢,(x, y),on C 


where ¢, (x, y) is the prescribed normal velocity on 
points on the boundary. [The particular functional 
form of ¢, (x, y) is dependent on the manner of expan- 
sion of the boundary.] Since elliptic bodies are treated 
herein, a restriction on the manner of expansion of the 
boundary is that the new boundary be an ellipse. The 
center of the ellipses are assumed to lie at x = y = 0; 
translation of the boundary, therefore, need not be 
considered. 

Inasmuch as the boundary conditions are to be satis- 
fied on an ellipse, the appropriate coordinate system 
for this problem is elliptic coordinates (see, for example, 
Milne-Thompson''). If 


x = ccosh écos 7] 
(6) 
y = csinh é sin »\ 


it can be shown that the curves — = constant are confo- 
cal ellipses and that the curves 7 = constant are confo- 
cal hyperbolae orthogonal to the ellipses, the foci of 
the ellipses and hyperbolae being at x = + c, y = 0. 
(see Fig. 2). The semimajor axis, a, and semiminor 
axis, b, of each ellipse are related by 


9 ~ 


c? = a’? — 5? (7) 
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The Laplace equation, Eq. (5), in elliptic coordinates 
transforms to (see, for example, Bateman!*) 


Pee + nn = 0 (8) 
and the boundary conditions become 
Og Og 
= | (n) on ~£&= & (9) 
O& o& 0 


Eq. (8) is solved by 
13) 


where & is the boundary ellipse. 
separation of variables (see, for example, Pipes 


yielding 
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+ = nt né - Ox pe e 
g=at+t+ort > (A,e™ + B,e~™) X és = = 6¢ b cos n bé 
n l O§& 0 
(C, cos nn + D, sin nn) (10) 4 i a (17 
6y = 6& = a sin n 6& 
O0E/ 0 


Since the velocity is zero at = o, the A,’s are zero 


and Eq. (10) becomes 


e¢=a+tat dX e (a, cos nn + b, sin mn) (11) 
n 1 


On the elliptic boundary, & = £, 


Og ng : 
( = A — =. ne ""(a, cosnn + b, sin nn) (12) 
OE 0 n= 1 


and the coefficients a, and 0}, in the Fourier series are 
found in the usual manner (by multiplying either by 
sin #7 or cos ny and integrating between the limits 0 and 


a = Lam 
27 Jo | 


27) as 


fie 
a, = ~ } cos nn dn (13) 
nT o£ 
Og 
b, = ry ( °) sin nn dn | 
nw OE / ) 


The coefficient co is evaluated in Appendix (A). 

Consider now the boundary conditions of the prob- 
lem. Since the boundary, which is at station s and 
intersects the fixed plane at time /, is denoted by C, let 
the boundary at s + ds, which intersects the plane at 
time ¢ + dt, be denoted by C; (see Fig. 1b). The 
equation of C is 


(x?/a?) + (y?/b2) = 1 


If vy is a normal to C directed outward and if dvt is the 
distance along this normal between C and Cj, then 


(14) 


(O¢g/Ov)o = (1/U)(Ov/dt) (15) 
Thus, 
(O0¢/0E)o) = (O¢G/Ov)o(dv/0E)o 
1 6p 
7 a (°°), (16) 


7 1H) (“) 
yi U 6& 0 dt 0 


where the subscript 0 indicates the value of the deriva- 


tive on the boundary. From Fig. 1b it is evident that 


(dv)? = (6x)? + (dy)? 
‘ eat 
and since} 


+t The symbol 6 preceding a quantity denotes the differential of 
that quantity along the normal between C; and C—i.e., the value 
at C; minus the value at C 


t From Eq. (6) 


(Ox/Of)o = c sinh {cos n = bcos 7 


(Oy/0F)o = c cosh fp sin » = a sin n 


(6n being zero) 
then, 
(6v/5£)o = a? sin? n + b? cos? n (18 


Differentiation of Eq. (14), considering a and 6 vari 


able, yields 


bx dx + a*y dy = a(b? — y?) da + bd(a’ x*)db 


This equation relates dx and dy from C to C, given 
(x, vy) on Cand da and db. Along the normal »y [using 


Eq. (17)], 


dx = 6x = bcos n 6 
dy = 6y = a sin n 6é 
or 
Pi ty 4 db. . 
- »— cos? n + a— sin? 7 
6£& dt dt , 
= (19 
dt a* sin? n + 6? cos? 7 


Substitution of Eqs. (18) and (19) into Eq. (16) yields 


@). 
of 0 


between 


l (0 ik oe ) . 
) ‘OS* — ea 2 
TV oy CO * a sin? ( 


The relation derivatives with respect to 


time ¢ and axial coordinate determined by 


body. 


s may be 


considering the boundary conditions of the 


Since the body surface is a stream surface, 
ov Oy /Ov (1/U)(6v/dt) 
+ (O¢/¢s) 


ds 1 + (d0¢/Os) 


aa + O(d*? In d*) (21 


20) can be written 


(52) da to 4 db... = 
= 0 cos* a sin- (22 
OE /o ds ds ’ 


Thus Eq. 


Eq. (22) may now be used in the evaluation of the coef- 
ficients of Eq. (13), which results in 


1 422 Sea i”) = S’(s) 
ied Ae d Qn 


ee 9: 
aa = ro ; (22a) 
4 ds\a 


a, = 0, ns TS 
Db, = @, w= 1.2.3.4 


ne 
wt 


Substitution of the 5) in 


Eq. (11) gives, for the potential, 


above relations and Eq. (A 
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given 
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(19 
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SLENDER BODIES OF EL 
ts c 
a (s+ Ve wt 
ei Qn 2 
aad (‘) ete a . . 
er" "CO (23) 
4 ds\a ’ 
PRESSURE DISTRIBUTION ON ELLIPTIC CONE 
The linearized form of the pressure equation is’ 
Ap 7 b— Pr os veil O¢ 7 (o°) 7 (Se) + 
g oP Os Ox oy 
O(d*4 In d*) (24) 


The potential for an elliptic cone is, from Eq. (23), 
since for a cone and for any other body with }/a con- 


stant [d(b a) ds = OJ, 


¢ = [S’(s)/2mr][E + In (c/2)] + do (25) 
The above equation indicates that, for elliptic bodies 
expanding so that b/a is constant, the instantaneous 
flow pattern is such that the ellipses confocal to the 
boundary (£ = constant lines) are equipotential lines, 
with the set of confocal hyperbolae (n = constant lines) 


being the stream lines. 


Of 
os cds 


= (1/c)(dc/ds)F(c/r, 0) 
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LIPTIC CROSS SECTION 

Clearly, the s = constant plane stream line and 
potential field patterns alter with time (or axial dis- 
tance s), so that at each new instant of time the bound- 
ary ellipse £ is an equipotential line of a new potential 
field. Since the value of £ on the boundary is con- 
stant [f depends, as shown by Eq. (31), only on the 
ratio b/a], the flow pattern in each plane is identical 
except for scale. This conclusion is not surprising, 
since the flow is conical—.e., the perturbation veloci- 
ties and pressures are constant along rays extending 
downstream from the cone vertex. 

Since the boundary is an equipotential line, the tan- 


gential velocity on the body in any plane s = constant 


is zero, and Eq. (24) becomes 
Ap C2) 
Ov 


oe 
C, = « — 9 — (26 
71 Os 
From Eq. (25), 
Oe 3 (s) Cc 
Os ~ Qn (« tine ‘) t 


S’(s) (= 4 ] “) 4b! = 
' Wis) (24) 
Os cds 


2a 


From Eq. (A-1) 


where F(c 7, 6) is defined by the above equations, and c/r is found on the boundary using the equation of the 


ellipse. 


From Eq. (4), or Ward,’ 


ho'(s) = 


Boundary values of the quantity F(c’r, @) are tabulated for several ellipse eccentricities in Table | 


I ves B ve 
S’’(s) In — S”(0) Ins — / 
2a 2 Jo 


(29 


In (s — oa) S’’"(e) ie 


TABLE | 


Boundary Values of the Quantity F(« 


r, 0) for Several Ellipse Eccentricities 


F(c/r, 0) 

Angle 6 e= 1 e= 7/; e = !/. e= !/; e='/, e=! 
0 1.0000 1.5000 2 0000 3.0000 4.0000 5.0000 
2 1.0000 1.4946 1.9690 2.7578 3.1060 2 9234 
} 1.0000 1.4800 1 8916 2 2436 1.9163 1.3835 
6 1.0000 1.4559 1.7747 1.7874 1.2325 0.80719 
& 1.0000 1.4244 1. 6401 1.3626 0.86925 0.55974 
10 1.0000 1. 3867 1.5005 1.0915 0.66852 0.43491 
12 1.0000 1.3449 1.3745 0.90581 0.54841 0.36417 
14 1.0000 1.3005 1] .2524 0.77529 0.47170 0.32045 
16 1.0000 1.2551 1.1478 0.68154 0. 42008 0.29162 
18 1.0000 1.2097 1.0906 0.61259 0.38380 0.27165 

20) 1.0000 1.1655 0. 98082 0. 56070 0.36616 0 2AT2Z8 
25 1. 0000 1.0633 0.81649 0.47662 0.31616 0.23506 
30 1.0000 0.97672 0.73685 0.42858 0.29348 0.22293 
4() 1.0000 0.84922 0.62229 0.38061 0.27069 0.21133 
50 1.0000 0.76841 0. 563238 0.35632 0.26029 0.20569 
60 1.0000 0.71814 0.53061 0.34741 0). 25487 0.20275 
70) 1.0000 0. 687938 0.51232 0.33768 0.25195 0.20117 
80 1.0000 0.67175 0. 50292 0.334385 0.25046 0.20042 
0 1 0000 0.66667 0. 50000 0.33333 0.25000 0. 20000 
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For a cone, S’’’(s) = 0 and S’’(s) = S’’(0), so that 
bo’(s) = [S’’(s)/2x] In [B/2s] 
If the ratio of minor to major axes of the cone is called 
e = b/a 


then, on the boundary, 


dc Cc da 
=“={Vi-@ 
ds s ds 


[from Eq. (7)] 
(1/c)(dc/ds) = (1/a)(da/ds) 


Now, for an elliptic cone, 


S’(s)/2ma = ea(da/ds) 30) 
vw? , e (oU ) 
S’'(s)/29 = &(da/ds)? 
From Eqs. (6), 
a l 
£ = cosh-!- = cosh! — 
c Vi1-é (31) 


=Inv( + «)/(1 — «) 
The above relations when substituted into Eq. (27) 
vield 
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Cc B Le 
- F(0) +i “wee 
r 4 ds { 


On the boundary, from Eqs. (14) and (15), 


Og (““) "| 
= € 41 
Os ds} \ 


Og da .je? + tan? 6 
=<“¥ : 
Ov ds Ye' + tan? @ 


The pressure coeflicient, Eq. (26), may now be ex. 
pressed, using Eqs. (82) and (33), as 


da\*{ AG B da 
(7) 41 eo) int 
as r é as 
oe tan @\| 
‘ 
e? + tan? 6/| 


where (cr, 0) is given by Eq. (2S) or Table 1. 


C, = 


Pressure distributions on two elliptic cones at zero 
angle of attack have been calculated for J/, 1.81; 
The cones are of ratio 
1.SS and 1/3. 
responding equivalent area circular cones have semi- 
and 11.0°, respectively. The 


the results are plotted in Fig. 3. 
of minor to major axes, ¢, of | The cor- 
vertex angles of 10.5 
geometry of the cones is given in Fig. 3. In that 
figure is also indicated the pressure distributions cal- 
culated by Ferri® with the linearized method of char- 
acteristics, as well as his experimentally obtained 
values. Agreement between the two theoretical meth- 
ods and experimental results is good, even in the case of 
the « '’; cone, whose major axis semivertex angle 
(18.4°) would seem to be greater than that permitted 
under the assumpt:ons of the slender theory. (See 


footnote on page 514.) 


PRESSURE DISTRIBUTION ON A PARABOLIC BODY OF 


ELLIPTIC CROSS SECTION 


Consider a body of unit length of elliptic cross section 
whose semimajor and semiminor axes are characterized 
by the following parabolic are relations: 


a(s) = 4Q@os(1 — s) 
b(s) = 4@os(1 — s) 


where 2@) and 2@,) are the maximum values of the 
major and minor axes, respectively. The ratio, ¢, of 
the elliptic sections is then constant, since 
b(s) Bo 
e= = = const. 
a(s) (to 


so that the potential is given by Eq. (25). If an effec- 
tive thickness ratio for the body is defined as 


d = 2V Q@ 


(i.e., the thickness ratio is the same as that for the body 


of revolution with equivalent area distribution); then, 


Y (28), and (29), 


using Eqs. (27), 
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Ov _ 4d2((1 — 6s + 6s?) In Bd(1 — s) + 3s(8s — 2)] + 
") 


, ' 
129 (1 — 6s + 6s*) In ne, = oe © ies 2s) 1 — r*, a) | - 
} 2Ve r 


O¢y 1 ; 
( a on 4a249 (1 — 6s + 6s?) In = : + (1 — 2s)7} 1 — Fé, 0) (34) 
Os { 2Ve r f 


( “ 


¢ 


where (Ov/Os), is the value of O0g/0s for the equivalent parabolic body of revolution [see Eq. (B-5)], and F(c/r, 6) 


is given by Eq. (28). Using Eq. (33), the normal velocity is found to be 


Oy/Ov = 2d(1 — 2s) V e(e® + tan? @)/(e4 + tan? 8) (35) 


The pressure distribution may now be computed using Eqs. (26), (34), and (35). 

The zero angle of attack pressure distribution on a parabolic arc body is shown in Fig. 4 for a Mach Number of 
2. The body has an effective thickness ratio, d, of 0.10; the ratio, ¢, of minor to major axes of the elliptic cross 
section is '/;. The dashed lines in the figure represent pressure coefficient values (these are independent of 6) of 
the equivalent parabolic arc body of revolution, which are obtained from Eqs. (B-S8) and (B-9) of Appendix (B). 
The variation of pressure coefficient with angle 6 at stations close to the nose and tail of the body is similar to that 
for the elliptic cone—.e., the pressure peaks at 6 = O° and is fairly constant for values of 6 between 30° and 90°. 
At the midpoint of the body (s = 0.5), however, C, is independent of angle 6 and varies only slightly with @ in the 
region close by on either side. The slight drop of C, near 6 = 0° at s = 0.4 and 0.6 is attributed to the predomi- 


nance of the O0g/Ov term in this region. 


WAVE DRAG OF BopIEs WITH BASES OF ELLIPTIC CROSS SECTION 


Ward’ has shown that the wave drag of a body of unit length is expressed by 


D l oe ate l = = S’(t) *! l sa 
— = | In S’’'(s)S''(a) da ds — In S""(a) do — 
(2/5)p, 2x Jo Jo -=— ~€ 2a 0 l—g@ 


~ me 
( / ¢g 7 ir) + 0 (d*§ In? d*) (36) 
JC Ov s 1 


where the contour C, for the third integral is the boundary at the base of the body. Since the first two integrals 
are independent of the cross-sectional shape of the body, it is apparent that the wave drag of arbitrary bodies dif- 
fers from the wave drag of the equivalent body of revolution by the difference of the value for both bodies of the 
third integral.t The third integral vanishes if the rear of the body is either pointed or cylindrical; thus, for 
these cases, the drag of a body of arbitrary cross section is equal to the drag of the equivalent body of revolution. 





It where R(1) and (1) are the base radius and elliptic 
[= Fa g(O¢ Ov) dr] , coordinate for the circular and elliptic body, respec- 
, . tively. Since the bodies are assumed equivalent [i.e., 
then the drag difference between bodies of equivalent Sos ; 
erage gia S.(s) = .38,(s)1, then 
area distribution having elliptic and circular bases is 
D, — D S’(1)}? 2R(1) 
(D. — D.)/q = Ile — Le ‘= ! In -— & (1) | + 
q 24 c(1) 4 
where the subscripts ¢ and e denote circular and elliptic le\? 
: T 4 ade 
bases, respectively. g 4 (| ) 
" ° i ‘ . . GS/s = 
Expressions for J, and /, are obtained from Eqs. ‘ae Eas. (30 1 (31 
& or, using Eqs. (30) and (31), 
B-1), (B-3), (18), (20), and (25) as , ‘ 
’ ; [S’(1)}? 2V €(1) 
S’(1) Coe — Ode = “orgy ™ 
I, = S'(1) In R(1) + b(1) 2mS(1) 1 + (1) 
or 9 9 
= fa(1)]*? /de\? : : te 
a : = ACp, + ACp (37) 
a“ fO'(1) c(1) | S8e(1) \ds/, 1 
Ie = S%(1)4 - &)(1) + In rs + bo(1) | — 
on as The difference in drag between equivalent elliptic and 


w a(“*) circular base bodies may be considered, therefore, as 

the sum of two quantities, ACp, and ACp,, as given by 
the first and second terms, respectively, of Eqs. (37). 
The quantity ACp, arises from the rate of change of 
area at the base, while ACp, results from the rate of 


ds) 


+ It will be noted that —(1/2)p S.¢ O0¢/Ov)dr] 1 is the lin 
earized value of the planar component of the kinetic energy at 


the station s = 1 
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change of eccentricity of the base elliptic boundary. 
Since 


9 J 
——- £aeees 1 
l+e 


Then, it is also evident that 


lt 


then ACp, < 0. 
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AC», > 0 


Therefore, the wave drag of an elliptic-base body with 
zero rate of change of eccentricity at the base |ie, 
(de/ds), ~; = 0] is lower than the wave drag of the 
equivalent body of revolution; a nonzero rate of change 
of eccentricity at the body base increases the drag of 
the body. 

The drag difference between equivalent elliptic and 
circular cones is plotted in Fig. 5 for a range of elliptic 
cone eccentricities.t The limiting case of an ellipse, 
the flat plate (« = 0), has no meaning as far as Fig. 5 js 
concerned, since the curves give the drag difference be- 
tween equivalent (in cross-sectional area) cones. The 
wave drag for « = 0 is, of course, zero. 

The curves of Fig. 5 may be used to determine the 
difference in wave drag coefficient of two arbitrary 
equivalent area bodies, one having a circular base and 
the other having an elliptical base, with zero rate of 
change of eccentricity at the base. For this applica- 
tion, the abscissa of the figure corresponds to the angle 
between the meridian plane body curve at the base of 
the circular body and the body axis. 

The wave drags of the two elliptic cones of Figs. 3a 
and 3b have been calculated by adding to the drags of 
the equivalent circular cones, Eq. (B-4), the differences 
given by Eq. (36). Drag coefficient values of 0.097 
and 0.098, respectively, are obtained as compared to 
those calculated by Ferri® of 0.103 and 0.099 for the 
same cases. This good agreement is not surprising in 
view of the good agreement obtained between the two 
theories for the pressure distribution. If the percent- 
age drag difference between equivalent elliptic and cir- 
cular cones is defined by the following equation 


, . Cp, _ Cb 
“> drag difference = ( C 100 
D 


then for the elliptic cone of Fig. 3b, for example, the 
drag reduction given by the results of the slender 
theory is about 10 per cent, while that according to the 
results of Ferri’s nonlinear theory is about 17 per cent. 
Since the elliptic cone drag given by both theories is 
nearly the same, the apparent wide discrepancy in per- 
centage drag reduction is due to the difference of values 
of drag of the equivalent circular cone given by the 
nonlinear and slender theories, the nonlinear drag 
being about 10 per cent higher than the slender drag. 
Hence, the drag differences (the numerator of the above 
equation) are widely different, being small differences 
between large numbers, and account, therefore, for the 
lack of close agreement between the two theories. It 1s 
suggested, then, that the drag of an elliptic cross-sec- 
tion body be calculated directly from the slender 
theory [Eq. (36)], rather than by adding the drag dii- 
ference between equivalent elliptic and circular cross- 

+ Eccentricity is defined in this paper as 1/e; the formal defini 


tion in analytic geometry is YY 1 — é 
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section bodies, Eq. (37), to the drag of the equivalent 
circular body as calculated by more accurate methods. 
REMARKS ON THE CALCULATION OF THE PRESSURE 
DISTRIBUTION ON SLENDER BODIES OF ELLIpTic CROSS 


SECTION AT ANGLE OF ATTACK 


In order to calculate the pressure coefficient on slender 
elliptic bodies at angle of attack, a, the additional pres- 
sure coefficient due to angle of attack may be super- 
posed upon the zero angle of attack pressure coeffi- 


cient as calculated from the present theory—.e., 


Ap/q = (Ap/q)a=0 + (AP/Qa 


The additional potential, g,, due to angle of attack 
(to order of the slender theory) is the two-dimensional 
incompressible potential of an elliptic boundary (which 
maintains constant area and eccentricity with time) 
translating in the reference plane with velocity La. 
This potential is a classical result and is given, for ex- 


ample, by Ward.’ 


CONCLUDING REMARKS 


A theory, based on the first-order solution of Ward’ 
for bodies of generalized cross section, has been pre- 
sented for the supersonic flow about bodies of elliptic 
cross section. The theory has been applied to the cal- 
culation of zero angle of attack cones and parabolic 
arc bodies of elliptic cross section; several example 
pressure distributions for these bodies are presented. 
The cone pressure distributions are compared to the 
nonlinear theoretical calculations and experimental 
results given by Ferri.6 Agreement is found to be good 
even for a cone that has a larger slope than that for 
which a slender theory appears applicable. 

The wave drag of bodies having elliptic bases has 
been derived. It is shown that the wave drag of an 
elliptic base body with zero rate of change of eccentric- 
ity at the base is lower than the wave drag of the equiva- 
lent area body of revolution. A nonzero rate of change of 
eccentricity at the base of the body increases the wave 
drag of the body. The wave drag of the two elliptic 
cones calculated also agrees very well with the nonlinear 
results of Ferri.* 

The Ward’ slender theory is also applied to cones and 
parabolic arc bodies of revolution. Comparisons are 
made with the results of other linear theories and (for 
the cone) the nonlinear theory. 

The application of the present elliptic cross-section 
body theory to bodies at angle of attack is also briefly 


indicated. 


APPENDIX (A) 


The coefficient co, which occurs in Eq. (11), will be evaluated herein. 


fon [OT Od 
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EVALUATION OF COEFFICIENT (Co 


Using Eq. (6), it may be determined that 
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~ 


~ 


= In (r/c) + ('/2) In G(c/r, 0) (A-J 
where G(c/r, 0) is defined by Eqs. (A-1). Eq. (11) may now be written [using Eqs. (22a) | as 
atde » : 
be ~ é€” Cos 27 
S"(s) * (S) | se t ds 
go = In r + Co + In . G : 6 = P 
2r 4a c r ra” s£ 
(") o(S0) 
C/ r 
S’(s) ‘ 
7 Inr + co + filr, 0) + fo(r, @) (A-9 
aT 
From the general solution of Ward's Eq. (3), the Thus, on the surface of the cone, 
above equation must have the following form: 
O¢ = If o\ Jp f , 
S's = [S’’(s)/2r] In r + dy 
n (S) — - Os 
= Inr +oo + do + > ar cos (n@ + B,) : dR\? Br 
27 n=1 = ( ‘) In . [using Eq. (29) | B-? 
(A-3) ds oo 
O0¢/Ov = dR/ds (B-3 


where the last two terms of Eq. (A-3) are equal to the 


last two terms of Eq. (A-2). Thus, since 
0 < |cos 27} < 1 


and since e~*’ is finite for finite &, then 


: : ; Sts) 2 
do = lim (fi + fe) = In (A-4) 
r= i 2a ( 
Comparing Eqs. (3) and (A-3) and using Eq. (A-4), 
there results 
oo ts) c : 
Co = —dy + bo = In. + dy (A-5) 
27 2 


with dp) given by Eq. (4). 


APPENDIX (B)-—— PRESSURE COEFFICIENT ON CIRCULAR 
CONES AND PARABOLIC BODIES OF REVOLUTION 


The pressure coefficient for circular cones and para- 
bolic bodies of revolution at zero angle of attack will be 
derived here from the Ward-type slender theory, and 
comparisons will be made with more exact theories. 

Since the potential ¢ for an axisymmetrical body is 
independent of the polar coordinate @ when the body is 
at zero angle of attack, the coefficients a, of Eq. (3) 
are zero and 


g = [S’'(s)/2r] Inr + 3b (B-1) 
For a cone of boundary radius r = R(s 
S(s) = wR? 
os) = 2rR(dR/ds) 
S’'(s) = 2x(dR/ds)? = const. 
S’’(s) = 0 


Since dR/ds = tan yp, where yu is the cone semivertex 
angle, then from Eqs. (B-2), (B-3), and (24) 


C, = —tan? p[2 In (B tan w/2) + 1] 
” Cp (B-4 


It will be noted that the above equation is equivalent 
to the first two terms of Broderick’s"” series solution 
(Broderick expands the potential, which satisfies the 
nonlinear potential equation, in a power series in the 
thickness and logarithms of the thickness) for the cone 
The Ward-type slender body solution, therefore, gives 
essentially the first-order solution to the problem. It 
is of interest to compare the circular cone pressure coef 
ficient given by this slender theory with the results of 
higher order linear theories and the nonlinear theory. 
This comparison is made in Fig. 6 for cone semivertex 
angles of 5, 10, 15, and 20 deg. The other theories 
used are described below and are listed in order of esti 
mated accuracy. 

(a) Nonlinear Taylor-Maccoll Theory.'*—The non 
linear theory is, of course, the most accurate of the 
theories listed. 
obtained from the M.I.T. Tables.'’ 

(b) von Karman-Moore Linearized Theory. 


The values plotted in the figure were 


von Kar 
man and Moore? demonstrate that a solution to the 
linearized potential equation, Eq. (2), which satisfies 
the boundary conditions for a body of revolution at 


zero angle of attack, is 
‘ - fia) da 
¢g(s,7) = = (B-5 
. V (S — oa)? — Br? 


The potential, thus, is that due to a supersonic source 
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distribution, of local strength f(s), extending down- 
stream from the origin along the body axis. The upper 
limit of integration indicates that only that portion of 
the line source in the fore Mach cone of the point (s, 7) 
affects the flow at the point. In the von Karman- 
Moore theory, the boundary conditions are satisfied 


exactly at the body surface. This fact requires that a 
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stepwise calculation be carried out for an arbitrary 
body of revolution; the line source is brokcn up into a 
number of discrete sources of constant value of f’(s), the 
integral [Eq. (B-5)] being replaced by a corresponding 
summation formula. The source strength for each 
step can be calculated once the strengths of the sources 
upstream of the given source are known. The solution 
for a cone can be obtained in closed analytical form; the 


perturbation velocities are” 


¢, = —K, cosh™! (cot w/B) 


¢g, = K,BV (cot u/B)? — 1 
where 


; tan yu 
kK, = 


V cot? « — B* + tan uw cosh™! (cot u/B) 


Some controversy exists at the present time as to 
the pressure coefficient equation to be used with a linear 
theory once the perturbation velocities are obtained 
Van Dyke! advo- 


cates the use of the exact pressure equation, 


from the equation for the potential. 


Ap 2 
= , x 
q 7M," 
j = 1 ~ e 
) l + > MO aed Yr" - © | = 7 (B-6) 


where y is the ratio of the specific heats. He recog- 
nizes, however, the inconsistency of its use, since the 
order of approximation of the linearized pressure 
equation [Eq. (24)] is not lower than that of the lin- 
earized potential equation. Results obtained using 
both pressure equations are plotted in Fig. 6. 

(c) von Karman Asymptotic Approximation to the 
Solution of the Linearized Potential Equation.°\—This 
solution is obtained from the von Karman-Moore? in- 
tegral formula for the potential, Eq. (B-1), by satisfying 
the linearized body conditions, Eq. (22), rather than 
the exact boundary conditions, on the body axis rather 
than on the body surface. With this assumption, the 


source strength can be shown to be 
f(s) = (1/2r)S’(s) 
This relation substituted into Eq. (B-5) gives, for the 


potential, 


] Ws Br S'(a) de 
Q\S,7) =; : a 
2 , V (s — ao)? — B?r? 


Thus, if the body area is given in analytical form, the 
above equation may yield the potential in closed 
It has been shown by Lighthill'® that 


(B-7) 


analytical form. 
Eq. (B-7) gives ¢ with an error of O(d). 
Eq. (3) has an error term of O0(d*), it is seen that the 


Thus, since 


von Karman asymptotic approximation leads to errors 

of the same order as the Ward-type slender theory. 
The pressure coefficient for a cone given by the von 

Karman theory, using the linearized pressure equation 
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{Eq. (24) ], can be shown to be 


4 ! 
P = tan? yp [2 cosh ( = ) = | 
B tan p 


In Fig. 6 it is indicated that the von Karman-Moore 
solution with the exact pressure equation agrees most 
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closely with the nonlinear results. As is expected the 
agreement of the linear theories is best for smaller cone 
angles; the agreement is also best at the lower Mach 
Numbers. The agreement between the von Karman 
asymptotic approximation and the Ward-type or 
Broderick first-order solution is close, as expected ex- 
cept at the higher Mach Number range and for the 
Here, the von Karman solution 
appears to be inferior. 

The slender theory pressure coefficient on a parabolic 


larger cone angles. 


body of revolution of unit length, whose surface is 
characterized by the equation 


R = 2d(1 — s)s 


(d is the fineness ratio of the body), may be derived 
using Eqs. (B-2) and (29), with the result that 


O¢/Os = 4d*[(1 — 6s + 6s?) In Bd(1 — s) + 


35(3s — 2)] (B-8) 


O¢/Ov = 2d(1 — 2s) (B-9) 


The pressure coefficient calculated from the Eqs 
(B-8), (B-9), and (26) for a body of fineness ratio of 
0.10 is presented in Fig. 7 for flight Mach Numbers 
of 1/2; 2, and 3. In that figure is also indicated results 
of calculations using the von Karman" asymptotic 
solution to the linear theory integral equation for ¢ 
[Eq. (B-7)], which was applied to the parabolic body of 
revolution by Jones and Margolis." The agreement 
between the two methods is, as is expected, very good. 
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Three-Dimensional Laminar 
Boundary-Laver Flow 


FRANKLIN K. MOORE* 
Lewts Flight Propulsion Laboratory, NACA 


SUMMARY 


The possibilities for predicting steady three-dimensional 
boundary-layer behavior are inferred from detailed consideration 
of the case of a cone at angle of attack to a supersonic stream. 

Qualitative features of the cone boundary layer are discussed, 
and solutions are presented in which (a) angle of attack is as 
sumed small (a method of perturbations), thickening of the 
boundary layer on the lee side of the cone is indicated but no 
separation of the cross flow occurs; and (b) consideration is re 
stricted to the plane of symmetry of the flow, with no restriction 
on angle of attack. On the lee side of the cone, unique solutions 
are limited to small angle of attack, and beyond a certain angle of 
attack, the boundary layer cannot remain thin 

The concepts of displacement thickness and flow separation are 
generalized to three dimensions, and failure of the thin boundary 
layer assumptions is related to the formation of strong vortices in 
a separated region, for laminar or, more tentatively, turbulent 


flow 


INTRODUCTION 


— TECHNIQUES FOR PREDICTING the steady plane 
laminar flow of a viscous fluid are rather well de 
veloped. Interest in aircraft applications has centered 
attention on the boundary layer, which is a thin layer of 
fluid enveloping a body, to which the frictional and dissi- 
pative effects of viscosity are confined and beyond which 
the fluid behaves in a nonviscous manner, substantially 
unaffected by the presence of the boundary layer. 

It is the purpose of the present paper to discuss the 
steady three-dimensional boundary layer. It is possible 
to attempt a general treatment, introducing particular 
problems for illustrative purposes only. This approach 
is largely frustrated by the great variety of aerodynamic 
effects which may arise in three-dimensional flows. The 
reverse approach is therefore adopted. Detailed con 
sideration is given to a specific problem, which is in- 
trinsically of interest, and from the results are drawn 
whatever general conclusions seem warranted. The par- 
ticular problem treated concerns the laminar boundary 
layer of a semi-infinite circular cone at angle of attack 
to a supersonic stream (see Fig. 1). This choice is made 
for the following reasons: 

(a) The semi-infinite cone is a useful idealization of 
the nose portion of a spike inlet or of the nose portion of 
the fuselage of a supersonic aircraft or missile. The 
aerodynamic performance of these configurations at 
angle of attack is of current interest. 

Presented at the Aerodynamics Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953 

* Aeronautical Research Scientist. 


(b) The limiting cases of small and large angle of 
attack illustrate available three-dimensional boundary- 
layer theories to which plane flow considerations are at 
least partially applicable. 

(c) Cases intermediate between these two limits are 
in no way characteristic of plane flow and yet are 


amenable to analysis. 


NOTATION 


a es functions of 6 and ,, see Eqs. (1) 

c = constant, see Eq. (3) 

c = specific heat at constant pressure 

f, g = flow functions, see Eqs. (5) 

h = overall local boundary-layer thickness 

= gas constant 

K = function related to 6*, see Eq. (12 

k = parameter, see Eq. (16) 

MU, M, = Mach Number and free-stream Mach Number, 
respectively 

p = static pressure 

gq = velocity vector in the boundary layer 

Rk,, R, = Reynolds Number (= p,;,x/p,, prttr/mr) 

i = absolute static temperature 

u,v, WwW = velocity components in direction of x, y, and ¢, 
respectively 

ch = coordinates along a generator of the cone, nor- 
mal to the surface, and circumferentially 
around cone, respectively, see Fig. 1 

a = angle of attack, see Fig. 1 

1 = ratio of specific heats 

oa = lengths based on mass-flow defects, see Eqs. 
(13) 

6* = displacement thickness of boundary layer 

6 = semivertex angle of cone 

r = similarity coordinate, see Eq. (6 

m = coefficient of viscosity 

p — density 

?, Vv, y = flow functions 


Subscripts 
1 = conditions in flow at outer edge of boundary layer 


r = reference condition 


Subscript notation for partial differentiation is used where 


convenient, 


Superscripts 
(n) denotes order of approximation in perturbation method 
Bar over quantity denotes evaluation of nonviscous flow at cone 
surface when a = 0. 
Primes denote ordinary differentiation. 
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Fic. 1. Circular cone at angle of attack-—notation. 


GENERAL FEATURES OF THE PROBLEM 


(1) Nonviscous Flow Outside the Boundary Layer 


A circular cone at angle of attack to a steady super- 
sonic stream is shown in Fig. 1. The nonviscous flow 
outside the boundary layer is assumed to be unaffected 
by the presence of the thin boundary layer and is 
conical in the Busemann sense (the shock is attached) ; 
all physical quantities are constant along rays from the 
apex and, in particular, along generators of the cone. 
Therefore, physical quantities at the outer edge of the 
boundary layer depend on ¢g but not on x. These circum- 
ferential variations govern the behavior of the boundary 
layer at its outer edge and therefore require specifica- 
tion. The most convenient representation of the non- 
viscous flow is provided by the M.I.T. tables.'~* Young 
and Siska‘ have elucidated the application of these 
tables at the cone surface (alternatively, at the outer 
edge of a thin boundary layer) and have made com- 
parisons with experiment. At the cone surface, 


u,/u4 = 1 — aA, cos ¢ + a*(B, + 

D, cos 2¢) + 
w/a = aA, sing + aD. sin2¢g+... » €D) 
p/p = 1 + aA; cos ¢ + a?(B; + 


D; cos 2¢) + -9 ete. | 


The A's, B's, and D's are functions only of 6 and .\/, and 
are obtainable from references 1, 2, and 3 via reference 4. 
Being the result of expansion in powers of a, Eqs. (1) 
may properly be used only for small or moderate a. A 
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Fic. 2. Typical circumferential static pressure distributions on 
cone at angle of attack. 


further limitation has been pointed out by Ferri,’ who 
shows that Eqs. (1) are correct only if the outer flow 
may be considered isentropic (neither 6 nor J/, too 
large). 

At zero angle of attack, the pressure is constant 
around the cone. For small positive angle of attack, the 
pressure is slightly higher at the bottom of the cone and 
slightly lower at the top, as shown by Eqs. (1) and Fig. 2. 
As the angle of attack is further increased, there is a 
region of rising pressure (‘‘adverse pressure gradient” 
near the top of the cone. Finally, in the limit [beyond 
the scope of Eqs. (1)] of a/@— ~, the pressure distribu- 
tion approaches the form appropriate to the plane cross 
flow around a circular cylinder. 


(2) Secondary Flow in the Boundary Layer 


A stream line at the outer edge of the boundary layer 
would follow a course around the cone as shown in Fig. 
3. Static pressure is constant across a thin boundary 
layer. Therefore, particles following a stream line 
within the boundary layer are subject to the same cir- 
cumferential pressure gradient as are those following the 
outer stream line but, because of their lower inertia, 
tend to take a course conforming more closely to the 
direction of the circumferential pressure gradient, as 
shown by the dotted stream line in Fig. 3. This is the 
well-known phenomenon of ‘secondary flow’’ in the 
boundary-layer transverse to the outer flow and is 
characteristic of the three-dimensional boundary layer. 


(3) Boundary-Layer Similarity 


It may be shown®* that the equations of motion, to- 
gether with the fact that outer flow quantities are inde- 
pendent of x, permit the supposition of parabolic simi- 
larity in the laminar boundary layer—that is, physical 
quantities in the boundary layer depend only on the two 
variables y/+/x and ¢. Accordingly, if the boundary 
layer is viewed in any meridional plane, flow quantities 
are constant along the family of parabolas originating 
at the apex of the cone and the boundary layer may be 
said to grow parabolically along generators of the cone, 


as indicated in Fig. 1. 


BOUNDARY-LAYER EQUATIONS FOR CONE 


In plane or axisymmetric boundary-layer problems, 
it is customary to introduce a stream function in order 
to satisfy the continuity equation identically. In 
three-dimensional problems, a similar device may be 
employed:’ In the coordinate system of the present 
problem, the steady continuity equation, 


+ This is a special case of the independent conclusions of 
Hayes® and the author’ that parabolic boundary-layer similarity 
may be assumed along generators of any ‘‘conical’’ body (e.g., a 
cone of elliptic cross section). For this class of flows, therefore, an 
analysis (either experimental or theoretical) at any surface 
spherical about the apex may be used to infer the behavior of the 
complete laminar boundary layer. 
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1 Oo ( 4 Opv ‘ ] Opw 0 
ux) = 
x OX P oy x sin@ Og 


js satisfied by the definition of new functions WV and 9, 


pux sin 0 = OW/dy 

; ov l vl 
pwr sin 6 = — = : (2 
or yx sin 6 " 


pwx sin 0 = OP/dy 


Definitions (2) are substituted into the remaining 
equations (momentum, energy, and state), under the 
assumption of a thin laminar boundary layer and with 
the use of certain transformations’— viz. : 

(a) A transformation due to Howarth confers the ad- 
vantage that the transformed equations closely re- 
semble the corresponding equations of incompressible 
flow. This transformation requires the assumption of 
the linear temperature-viscosity relation of Chapman 
and Rubesin :° 


nla = O/T, (3) 


where subscript y represents a reference condition and C 
is a constant chosen to match Eq. (3) to the more pre- 
cise Sutherland formula at the body surface. 

(b) A transformation based on that of Mangler® 
vields equations in nearly Cartesian form, 

(c) A transformation expressing parabolic similarity 
along cone generators. 

Further, it is assumed that the Prandtl Number is I, 
that the ratio of specific heats is constant, and that there 
is no heat transfer through the cone surface. These 
assumptions permit the equation of constant enthalpy to 
be used in the boundary layer in place of the more com- 
plicated differential energy equation. 

The following boundary-layer equations result:’ 


10 


momentum, 


p'(¢) 2 
[ 5 ms s+ Par IE a 
3p sin @ 3 sin @ 


5) 5) 
3 sin 6 Er hyo + = (g.)> + 2fra = 0 (a) 
» Si ) 


— p'(g) 2 2 
S24 =~ 51s =O ~ 
3p sin @ 3 sin 6 3 sin 6 


a ; z p’(¢) 
~fhih T tua * =>; (4b) 
3 3u*,sin@ p 


enthalpy, 


u,* 


= [(fy)? + (g)?] = C71 + 


State, 


p(= pr) = pJT (4d) 


Where f(\, ¢) and g(A, ¢) are related to the functions 
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a 
Qos 


—— > TYPICAL STREAMLINE AT OUTER 
EDGE OF BOUNDARY LAYER 

——- TYPICAL STREAMLINE INSIDE 
BOUNDARY 


Fic. 3. Secondary flow on cone at angle of attack 


WV and ® of Eqs. (2), 
u/u,=f,; w/uUr =2y (5) 


The similarity coordinate \ is a function of y/ Vx, 


ee : *Ry/¢ 
- (37 ) P d(R,/C) (6) 
p R, 0 Pr 


The boundary conditions are: At the outer edge of the 
boundary layer (A ~ ©), the u and w velocity compon- 
ents should take on their corresponding nonviscous 


r 


values, 
(¢) 


fy ( o,g) = : { > 
u, u, 


At the cone surface (A = 0), u« and w vanish, 
f,(0, ¢) = g,(0, ¢) = 0 (7b) 
and the normal component v vanishes, t 
f(0, ¢) = g(0, ¢) = O (7c) 


Eqs. (4) may be reduced to a pair of simultaneous 
equations with two independent and two dependent 


t From Eqs. (2), the simplest way to make v = 0 at the cone 
surface is to require V and ®to vanish separately. This point is 


discussed more fully in reference 7. 





KARMAN VORTEX 
STREET - 





BOUNDARY LAYER - K / 


SIDE VIEW END VIEW 


Fic. 4. Slender cone at large angle of attack. 
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variables and, despite the simplifications already made, 
present a rather formidable problem in general. There- 
fore, attention will be given to the special simplifications 
available in limiting cases. 


LIMITING CASES 


1) Large Angle of Attack, Small Cone Angle 


In the special limiting case illustrated in Fig. 4 for 


which a/6— ~ (in effect, @— 0 in order that the shock 
may be considered attached), if the transformation 


A /V3sin6d; F f/V3sin6; G=e/V3siné 


is substituted into Eqs. (4) and if terms of order 0/a are 
neglected, Eqs. (4) reduce" precisely to the equations 
for the laminar boundary layer on an infinite cylinder 
at angle of attack. Similarity exists in terms of the 
square root of the local radius of the cone (x sin @), 
Thus, if the cone is sufficiently slender and the angle of 
attack is sufficiently large, it is correct to suppose that 
the boundary layer at any point is identical with the 
boundary layer that would exist on an infinite cylinder 
of constant radius equal to the local radius of the cone 
at that point. For the case of the infinite cylinder at 
angle of attack, Prandtl,'’ Sears,'? and Jones'* have 
shown that the profiles of velocity components lying in a 
cross-sectional plane may be obtained by plane-flow 
methods directly, as though the motion were that about 
a cylinder in a plane flow with a stream velocity equal to 
the cross-flow component of the actual stream (see Fig. 
4). Although this result is strictly correct only for in 
compressible flow, the approach has been successfully 
exploited by Allen and Perkins!‘ for the estimation of 
viscous cross-forces on slender bodies of revolution at 
large angle of attack in supersonic flow. 

Thus, when a 6— ©, the circumferential cross flow 
on a cone at angle of attack may be expected to exhibit 


oye ”" + 2fe srr 
zero order: 
[fo"(@) = 1; 


Pe ere ea 
o 
g(D'(o) = 1; g)'(0) = g (0) 
first order: and 
2A» 


f'(0) = f(0) = O04 
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plane-flow characteristics in particular, separation in 


the form of a von Karman vortex street f (Fig. 4). 


(2) Small Angle of Attack, Large Cone Angle 


Mangler® has shown that there is a transformation of 
variables which relates the laminar boundary layer on a 
body of revolution at angle of attack to an equivalent 
plane flow about (say) an airfoil section. A special case 
relates the laminar boundary layer on a circular cone at 
zero angle of attack in supersonic flow with the plane 
In effect, 
Mangler’s transformation has been included in Eqs. (4), 


boundary layer on a semi-infinite flat plate. 


which therefore reduce to the Blasius flat-plate equation 
in the limit a/@—> 0. 


A MeEtTHop OF PERTURBATIONS 


1) Differential Equations 


In certain problems of the three-dimensional boundary 
layer, solutions may be obtained by supposing that the 
In the 
present instance, a solution is sought for small angle of 
attack (a/@ < 
symmetric flow at zero angle of attack. 


flow is nearly plane or nearly axisymmetric. 


1) as a perturbation of the known axi- 
The functions 
fand g.are written in a form chosen to be consistent (to 
first order in a) with the outer flow as represented in 
Eqs. (1), 


f f() — aA; cos gf P(A) +...) 

g = aA, s1n gg™(A) +... \ 
Incorporation of Eqs. (8), (4c), and (4d) into Eqs. (4a 
and (4b) and into boundary conditions (7) and collect 


ing terms of zero and first order in a /@ separately yields 


7 An experimental study of shed vortices on the lee side of 


slender bodies of revolution is provided by Allen and Perkins 
0 } 
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(The reference condition has here been taken to 
correspond to conditions in the nonviscous flow at the 


cone surface when a = 0.) 


2) Results 


Successive solution? of the ordinary boundary-value 
problems (9), (10), and (11) provide a complete de 
scription of the laminar boundary layer on a cone at 
small angle of attack. When a = 0, the boundary-layer 
thickness is constant in any cross section. When there 
isa small angle of attack, the boundary layer 1s slightly 
thicker on the top and thinner on the bottom, because of 
the circumferential component of flow illustrated in Fig. 
3. The formula for displacement thickness is 


(a/A)K(0, MT.) cose +...) (12) 


The analysis yields A, which is shown in Fig. 5, to 
gether with a cross-sectional sketch of the circumferen 
tial variation of displacement thickness due to angle of 
attack. Roughly, the change in 6* due to angle of 
attack is in the same proportion to (6*),, 9 as a@ is to @. 
Further examination of the results of this perturba 
tion procedure yields no indication whatever of cross 
flow separation; this result might be anticipated in view 
of the favorable nature of the circumferential pressure 


gradient indicated for small angle of attack in Fig. 2. 


DISPLACEMENT THICKNESS 


The definition of displacement thickness [Eq. (12) | of 
a three-dimensional boundary layer requires clarifica 
tion. 

The boundary layer is normally considered as a region 
wherein viscosity effects predominate and outside of 
which nonviscous laws apply. When the Reynolds 
Number is large, the boundary layer is assumed to be 
so thin that the nonviscous portion of the flow occurs as 
if there were no boundary layer. This assumption is 
strictly correct for infinite Reynolds Number. 
but finite Reynolds Number, however, the growth of the 


For large 


boundary layer causes the nonviscous flow to be de 
flected (“‘displaced"’) away from the body surface. 

This displacement effect of a plane boundary layer is 
analyzed by equating the actual mass flow in the bound- 


arv laver 


| pu dy (h outer edge of boundary layer 
J/0 
to the mass flow associated with a (fictitious) com- 
pletely nonviscous flow terminated at 6* rather than at 
the body surface, 
pilty(h — 6*) 

vielding the usual definition of displacement thickness 
in terms of the ‘“‘mass-flow defect’’ of the boundary 
laver, 

+ Details are given in reference 10. The solution of Eqs. (9 
is the well-known Blasius function. 


BOUNDARY-LAYER FLOW 529 


*h 

. pu 

6* = | (1 - ) dy 
0 pil 


In three-dimensional flow, mass-flow defects occur 
with reference to both velocity components tangential 
to the surface; in the case of the cone at angle of attack, 


these are 


. sa pw 
6, / (1 ) dy (13 
J0 piv 


and it is not clear that either describes the displacemeyit 
effects of the boundary layer 

A unique definition can be given" for the displace 
ment thickness of a three-dimensional boundary layer 
At the far outer edge (4) of any boundary layer, there is 


(9 


a certain distribution of the velocity component (¢ 
normal to the body surface, determined by the charac 
teristics of the boundary layer. The displacement thick 
ness 6* may then be defined as the local height of a 
fictitious impermeable surface that would deflect a 
completely nonviscous flow in just such a way as to 
produce at y = /a distribution of normal velocity equal 
to that actually produced by the boundary layer. The 
continuity equation and the thin boundary-layer as- 


sumptions provide the equation for 6*, 


div | pigi6* — / (pig: — pq dy 0 (14 


(The divergence operator and the vector g have two 

components, both tangential to the surface.) Eq. (14 

applies to turbulent, as well as laminar, boundary layers 
In the case of a cone at small angle of attack, Eq. (14 


simplifies to 
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2A» 


6* = 6,—a : (6, — 6,) cos¢ + 
3 si 6 
which is equivalent to Eq. (12). 
Thus, 6* depends on the two mass-flow defects defined 
in Eqs. (13) but is not equal to either. 


BOUNDARY LAYER IN PLANE OF SYMMETRY 


(1) Differential Equations 


With reference to the important question of cross-flow 
separation, the analysis for small angle of attack pro- 
vides only the information that the phenomenon must 
be sought at higher angles of attack than are contem- 
plated in the perturbation analysis. A complete analysis 
would be difficult to carry out for large angle of attack. 
However, in common with many three-dimensional 
flows, the boundary layer on a cone at angle of attack 
has a plane of symmetry; the vertical meridional plane 
¢ = 0, m. (See Fig. 1.) Eqs. (4) and (7) simplify" to 
ordinary differential equations if consideration is re- 
stricted to this plane, where p’(¢) = f/,, = g, = 0, 


(f + ky)f" + 2f'" =0 
er + 


2 — 1 5) 
(« + ait “ Met — "if — 


(f + kv" + 20" — RI)? — 5 


» 


Nl o)= y’( o) = |] 
f'(0) = y’(0) = f(0) = WO) = O 


whére y is a new dependent variable and k is a parame- 
ter, 

[wi (O, w)/us(O, wr) |W(A) = g (A; 0, 2) 
k = [2/3 sin 6)(w,,(0, )/u,(0, 7) } (16) 


(The reference condition is taken at the outer edge of 
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the boundary layer at g = 0 or 7, depending on whether 
bottom or top of the cone is under consideration.) The 
parameter k depends essentially on a; the greater the 
angle of attack the greater the magnitude of w,,, and, 
hence, k. At the bottom of the cone k is positive; at 
the top, negative. 


(2) Results at Bottom of Cone 

Solution of Eqs. (15) yields, in the plane of sym- 
metry, the profile of the velocity component along gen- 
erators 


u/u, = f'(A) 
and the circumferential gradient of the circumferential 
velocity in the boundary layer, as a fraction of the 
corresponding gradient just outside the boundary layer, 
W/W = W'(X) 


At the bottom of the cone, velocity profiles have been 
obtained'*® for various angles of attack, 1/7, = 3.1, and 
C= 15". 


10. The results (Fig. 6) show the expected thinning of 


The outer flow was represented as in reference 


the boundary layer as angle of attack is increased. The 
fact that the 
boundary layer exceeds its stream value is presumably a 


circumferential velocity within the 


) 


consequence of the secondary flow illustrated in Fig. 3. 


(3) Limitations of Method When Applied at Top of Cone 


(a) Uniqueness.— At the top of the cone, solutions 
can be obtained only for rather small values of a. Be- 


yond the angle of attack for which 


Vin ; 
k=—- ( ,) (a = 0.5° when M, = 
a y mar 


solutions exist but are not unique.'® The following is 
offered in explanation of this fact: 

At the bottom of the cone, fluid is entrained by the 
thin boundary layer and is conveyed around the cone 
Here 
the condition of the boundary layer is uniquely deter- 


away from the plane of symmetry (see Fig. 7a). 


mined by the characteristics of the outer flow at the 
plane of symmetry, and a thin boundary layer solution 
is unique, for all angles of attack. 

At the top of the cone, for large a, fluid arrives in the 
vicinity of the plane of symmetry after having traveled 
around the body. Therefore, the boundary layer at the 
top depends not only on the outer flow in the plane of 
symmetry but on the influences to which the boundary 
layer has been subject during its passage around the 
body. A plane-of-symmetry analysis fails to yield a 
unique solution under these circumstances. 

However, the true lateral region of influence of the 
plane of symmetry grows parabolically (Fig. 7b) along 
the top generator according to the law of simple molecu- 
lar diffusion applicable in the present case. Therefore, 
at the top of the cone, while the flow is toward the plane 





of sv 
tram: 
angl 
regi¢ 
goill 
ax 
plan 
of tl 

h 
atta 


bk 


no s 
ary 

Hay 
bou 
Whi 
con 
sim 
[Eq 
tion 
cal : 
the: 
tha 
ince 
crit 
tiot 


bo 
Vo 
th 
di 


in 


ether 

The 

‘'r the 
, and, 


syin- 
gen- 


ntial 
~ the 
iver, 


been 
and 
nce 
g of 
The 
the 











LAMINAR 


of symmetry for all angles of attack, in the sense of being 
transverse to cone generators, for a sufficiently small 
angle of attack the flow is outward from the parabolic 
region described, and uniqueness is expected. The fore- 
going argument yields Eq. (17), which gives the 
maximum angle of attack for which conditions in the 
plane of symmetry uniquely describe the flow at the top 
of the cone. 

(b) /:xistence. At the top of the cone, for angles of 
attack greater than that for which 


k = lia = 6.2° when M, = 3.1,0 = 7.5") (18) 


no solution of Eqs. (15) exists which satisfies the bound- 
ary conditions. This fact was first poimted out by 
Hayes® in a discussion of the character of the laminar 
boundary layer on conical bodies in supersonic flow. 
Whereas the question of uniqueness discussed above is 
concerned only with the usefulness of a mathematical 
simplification, the existence of a critical angle of attack 
[Eq. (18)] beyond which a thin boundary-layer solu- 
tion fails to exist would seem to be of more direct physi- 
cal significance. If a thin boundary layer does not exist, 
then it follows that the actual boundary layer is so thick 
that the assumption of a thin laminar boundary layer is 
incorrect. Following Hayes,® one suspects that the 
critical condition (18) is connected with laminar separa- 


tion of the circumferential flow. 
SEPARATION 


1) General Considerations 


To form a proper explanation of the critical condition 
(1S), one must have a general view of the sort of in- 
formation concerning three-dimensional separation 
which thin boundary-layer theory can provide. 

In plane flow, the ‘‘separation point” is taken as the 
point at which a reverse-flow velocity profile appears, as 
This criterion of separation has proved to be 
However, 


in Fig. Sa. 
difficult to generalize to three dimensions. 
without specifying the velocity field inside the separated 
region, it may be observed that at the separation point 
the wall stream surface bifurcates, and at a point of re- 
attachment (if such occurs) the two stream surfaces 
join again at the wall. Thus, as in Fig. Sa, there is em- 
bedded at the base of the boundary layer a distinct 
bubble that does not exchange fluid with the rest of the 
flow. Generalization to three dimensions may therefore 
be accomplished by defining a separated region as a 
bubble of fluid embedded in the boundary layer be- 
tween the solid body and a stream surface meeting the 
body in a closed curve and containing a sheet pattern of 
vorticity the precise nature of which would depend on 
the particular problem. This sort of separation is pre 
dictable by thin boundary-layer analysis. 

Of course, the sort of separation of greatest engineer 
ing importance occurs when the embedded vortex sheet 
coalesces to form strong concentrated vortices, with the 
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(b) MOTION OF VORTICES IN THICK BOUNDARY LAYER. 


Fic. 8. Plane separation 


consequence that the outer flow is greatly disturbed and 
large pressure effects occur. This process is not ame- 
nable to an analysis that presumes a thin boundary layer. 
However, thin boundary-layer theory may predict the 
circumstances under which coalescence occurs; if such 
an analysis predicts a thin boundary layer throughout 
the separated region, the associated vortex sheet would 
presumably tend to adhere to the body. If, however, 
the theory indicated infinite boundary-layer thickness 
anywhere in the separated region, there the coalescence 
process would most especially begin, finally engulfing 
most of the separated region, 

The process of coalescence may perhaps occur by the 
following mechanism: Let the vortex sheet of Fig. Sa be 
represented by a series of individual vortices as in Fig. 
Sb and represent the effect of the wall by the image 


vortices below the wall. Then, if the fluid above the wall 
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imposes no additional constraint (i.e., a boundary layer 
of locally infinite thickness), each vortex would move to 
the left toward the separation point under the influence 
of the induced field of its image.+ Vortices initially near 
the separation point tend to remain fixed, however, be- 
cause there the boundary layer is supposed to remain 
thin. Therefore, each vortex moving upstream tends to 
overtake the vortex ahead of it, and coalescence into a 
single strong vortex ensues. 

A general criterion for the existence of a thin boundary 
layer would express the requirement that the boundary 
layer entrain fluid at its outer edge. If h represents the 
local overall thickness of the boundary layer and if it is 
recalled that the flow in the vicinity of 4 may be repre- 
sented as the nonviscous flow over a fictitious displace- 
ment surface 6*, the criterion may be written to specify 
that the fluid flux into the region between the surfaces 
h and 6* in the fictitious nonviscous flow be positive. 


div[(h — 6*)piq] > 0 (19) 


2) Examples 


(a) Plane Cylinder.—Taking x along the body, in- 
equality (19) may be written 


(d/dx)(h — 6*) — Opi, /Ox 


(20) 
h — 6* pill 


As the rear stagnation point is approached, u, —~ 0 and 
—du,/dx is finite and positive. Since h — 6* is by 
definition greater than zero, inequality (20) can be satis- 


fied only if 

d(h — 6*)/dx ~ 
in violation of the thin boundary-layer assumptions. 
This conclusion is reached with no detailed knowledge of 
the boundary layer and is, of course, obvious on the 


grounds that the boundary-layer fluid conveyed around 


+ A familiar example of this induced motion is the spreading of 
a smoke ring that impinges on the wall of a room. 
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Fic. 9. Separation on cone at angle of attack. 
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the body must finally leave the body in the vicinity oj 
the rear stagnation point and proceed downstream. 
vortex sheet would be predicted to cover the after-por 
tion of the body and, because of the infinitely thick 
boundary layer at the rear stagnation point, would 
coalesce to form a vigorous vortex system. 

In three-dimensional flow, consideration of the outer 
flow is, in general, not sufficient to predict the certain 
“blow-up” of a thin boundary layer; information con 
In the 
present instance, knowledge of parabolic similarity 


cerning the boundary layer itself is required. 


along generators permits inequality (19) to be written, 
near the top of the cone, 


(3/9 


—(dw,/dg) — (3/2), sin 0 


(0/O¢)(h — 6*) ” 
h — 6* Wi 


indicating infinite boundary-layer thickness only ii 
w, = 0 and, even then, only if dw,/d¢ is less than zero 
(true at the top of the cone) and is larger in magnitude 
than 3/2, sin @ or if k< —1 [ef. Eq. (18) ]. 
may be explained as follows: For large angles of attack 


This result 


at the top of the cone, the boundary-layer fluid must 
leave the vicinity of the body and proceed downstream 
(Fig. 9a), just as in the case of a plane evlinder. How 
ever, the boundary layer is presumed to grow parabolic- 
ally in the streamwise direction along generators of the 
cone. Thus, below some critical angle of attack (k = 
—1, as it turns out), the fluid conveyed around the 
cone toward the top may simply augment the growing 
boundary layer (Fig. 9c). 

As the angle of attack is increased from zero, at first 
there is no vortex sheet and the boundary layer is every 
where thin. It may be shown that an adverse pressur 
gradient first appears at the top of the cone when k 
—2/3 (a = 4.3° when M, = 3.1, 6 = 7.5 


then, an imbedded vortex sheet appears in the range 


Possibly 


—1 < k < —2/3 for which the boundary layer may be 
regarded as thin. This possible (not proved) state of 
affairs is illustrated in Fig. 9d. 

When the angle of attack is greater than that for 
which k = —1, the prediction of infinite boundary-layer 
thickness implies the freedom sufficient for coalescence 
of the vortex sheet. In the flow about a cone at moder- 
ate angle of attack, this coalescence gives rise to a 
steady symmetric vortex pair appearing as lobes at the 
top of the cone (Fig. 9b). Allen and Perkins'' provide 
experimental demonstration of this type of separation, 
which is an essentially three-dimensional phenomenon 
to be distinguished from the unsteady vortex street 
(Fig. 4) produced in plane flow or on the cone when 
ad> 1]. 

Fig. 10 shows how the critical angle of attack (k = 
—1) varies with M/, and 9. The critical angle of attack 
is seen to be of the same order as the semivertex angle of 
the cone. The figure indicates the possibility of rather 
profound qualitative differences in the flows at high 
Mach Number about cones of different vertex angles. 
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Of course, the argument given herein for the signifi 


cance of the critical point Rk = —1 is idealized. Physical 
significance should be attached not so much to the 
mathematical singularity at k = —1 as to the tendency 


to infinite thickness and the consequent tendency of a 
vortex sheet to coalesce when k approaches —1. Fur- 
thermore, possible effects connected with laminar in 
stability are not considered. 

(b) Conical Bodies in General. Inequality (19) and 
the condition of parabolic boundary-layer similarity may 
be used to provide a criterion [equivalent to Eq. (18) | 
for the maximum angle of attack consistent with a thin 
boundary layer on any smooth conically symmetric 
body (e.g., a cone of elliptic cross section) in supersonic 
flow. 

(c) Turbulent Boundary Layer on a Cone. 
larity also holds for the turbulent boundary layer on a 
cone and if the similarity law is nearly linear because of 


If simi 


the strong mechanical diffusion, vigorous separation 
would be indicated to occur at an angle of attack for 
which k —4/3 (a higher angle than in the laminat 
case, for which the critical value is k = —1). 


CONCLUDING REMARKS 


1; Summary Description of Cone Boundary Layer 


At zero angle of attack, the laminar boundary-layer 
growth along generators may be studied by plane-flow 
methods. As the angle of attack is increased to a small 
positive value, the boundary layer tends to thicken at 
the top of the cone and to thin at the bottom in conse- 
quence of the circumferential flow due to angle of attack. 
Also, a circumferential component of skin friction pro- 
vides a small viscous lift. There is no separation of this 
circumferential flow. 

When the angle of attack is increased to the value for 
which k = —2/3 (about 4.3° when @ = 7.5°, M, = 3.1), 
an adverse circumferential pressure gradient first ap- 
pears at the top of the cone. Presumably, at somewhat 
higher angle of attack, there appears a separated bubble 
embedded at the base of the boundary layer, growing 
in extent as the angle of attack is further increased. 

The boundary layer may no longer be regarded as thin 
when the angle of attack is reached for which k = —1 
(about 6.2° when @ = 7.5°, W/, = 3.1). At and above 
this angle of attack, whatever flat vortex bubble exists 
at the top of the cone must be in the process of coales- 
cence into a symmetric pair of strong steady vor- 
tices. 

If the cone is sufficiently slender, at an extremely high 
angle of attack (a/@ > 1), the circumferential flow is 
qualitatively the same as the plane flow about a cylinder, 
and a von Karman vortex street is shed at the top of the 
cone. 

Similar considerations would apply for the general 
class of bodies with conical symmetry and for the fully 


turbulent boundary layer. 
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2) The Use of Boundary-Layer Theory 


In the treatment of three-dimensional problems, r¢ 
duction in the number of variables is of first importance. 
The number of dependent variables may be reduced by 
making various physical approximations and by defining 
a pair of functions such that continuity is satisfied 
identically. The number of independent variables (co- 
ordinates of the solution) are reduced by 

(a) Invoking similarity when possible (e.g., parabolic 
similarity in the cone case). 

(b) Perturbing a known flow (e.g., a 6< 1, in the 
cone case). This results in linear equations to which 
the method of separation of variables may often be 
applied. 

(c) Restricting consideration to a plane of symmetry 
of the flow. As in the cone case, this method succeeds 
only when the boundary layer “‘starts”’ at the plane of 
symmetry. 

It is felt that laminar boundary-layer analysis may 
prove more useful and realistic for many three-dimen 
sional flows than for plane flow. For example, in cases of 
plane flow about cylinders, the question usually ad- 
dressed to theory is: “At what point on the body does 
separation always occur?” The actual separation often 
so distorts the outer flow that boundary-layer calcula- 
tions based on a theoretical pressure distribution do not 
On the other hand, in the case 
“At what angle of 


give a reliable answer. 
of the cone, the question may be: 
attack does separation first appear anywhere on the 
cone?’’ The answer to this question may reliably be 
given by theory because, at its first appearance, separa 


tion is apparently not violent. 
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Critical Bending Stress for Flat Rectangular 
Plates Supported Along All Edges and 
Flastically Restrained Against Rotation 
Along the Unloaded Compression Edge 


JAMES H. JOHNSON, Jr.,* anb ROBERT G. NOELT 


North American Aviation, Ine. 


ABSTRACT 


A theoretical analysis of the critical bending stress for flat rec- 
tangular all edges and elastically re- 
strained along the unloaded compressive edge is presented. Criti 
cal buckling constants for four positions of the axis of zero stress 


plates supported along 


are computed for all values of edge restraint and wave length 
These four positions are locations of the axis of zero stress, vary- 
ing from mid-depth of the plate to the unloaded tension edge 
The latter location of the zero stress axis is that of triangular com 


pressive loading 
A discussion of the effects of elastic restraint against rotation 


along the unloaded tension edge is given 


INTRODUCTION 


I THE DESIGN OF aircraft structural components, 
consideration must sometimes be given to the local 
instability of plate structures that are subjected to un 
symmetrical bending loads. For example, unsym 
metrical channel, I- and Z-section beams, fuselage rings, 
and multiweb box beams frequently buckle locally be 
cause of unsymmetrical distribution of bending loads. 
The local instability stress of a structure composed 
of two or more plates joined at the unloaded edges can 
* In 
order to use this method, critical stress coefficients k 
must be known as a function of the edge-restraint co 
efficients « and the half-length of buckle waves X for all 
These data are already 


be found by the moment distribution method.* 


loaded plates in the systein. 
given for plates in compression® ~? and for plates sub 
jected to symmetrical bending in the plane of the plate.” 
Also, the moment distribution method has been used 
to determine the local instability stresses for columns 
of various cross-sectional shapes* * and for multiweb box 
beams subjected to symmetrical bending loads.” 

In the present paper, critical stress coefficients are 
given for all values of \ and e for flat plates subjected 
to the general bending load of Fig. 1. The values pre 
sented herein can be used with the results of references 
2,5, 6, and 7 and the moment distribution method* 
to determine the critical stresses of flat-plate structures 
subjected to loads resulting in unsymmetrical bending 
stress distributions. 

Received July 15, 1952 

* Research Engineer, Aerophysics Laboratory 

T Stress Analyst, Aerophysics Laboratory 
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SYMBOLS 


flexural rigidity of plates, in.Ibs. |#t°/12(1 — yu?) 

Young’s modulus of elasticity, lbs. per sq.in. 

intensity of compressive force in x direction at edge 
of plate x = Oand x = X, lbs. per in 

intensity of distributed forces in x direction at edges 
of plate x = Oand x = X, lbs. per in 

stiffness per unit length of elastic restraining medium 
or amplitude of sinusoidal applied moment divided 
by amplitude of resulting sinusoidal rotation of 
elastic medium in quarter radians, Ibs.-in 

external work of applied forces during buckling 

internal energy of deformation during buckling 

length of plate in x direction, in 

Fourier coefficients of sin (mmy/b 

width of plate in y direction, in 

integers 

number of half waves in x direction 

nondimensional critical-stress coeflicient that de- 
pends upon conditions of edge restraint and shape 
of plate 

thickness of plate element, in 

plate-buckling deformation, normal to plane of the 
plate 

plate coordinate in direction of stress 

plate coordinate, perpendicular to direction of stress 

numerical factor that depends upon the location of a 
plate axis of zero stress 

ratio of half-wave length to plate width (A/b 

Lagrangian multiplier 


Poisson's ratio 








FIGURE |. 


EDGE 


INFINITELY LONG FLAT PLATE UNDER BENDING FORCES-TENSION 


HINGED, COMPRESSION EDGE ELASTICALLY RESTRAINED 
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d = half-wave length of buckles in longitudinal direction, 
in. 
€ = nondimensional restraint coefficient given by (459) + 
D) 
a = longitudinal compressive stress, lbs. per sq.in. 
6 = amplitude of sinusoidal edge rotation 
Subscripts 
p = plate 


y 


restraint 


ll 


cr = critical 


THEORETICAL ANALYSIS 


The energy method was used for the calculation of the 
critical stress of an infinitely long flat plate under bend- 
ing forces in the plane of the plate. The plate was con- 
sidered to be simply supported along the tension edge 
and elastically restrained by a restraining medium along 
the compression edge. It was assumed that the re- 
straining medium was such that a sinusoidal restraint 
produced in-phase sinusoidal rotation of the restrained 
edge. This condition is satisfied by a rigid joint be- 
tween two or more flat plates that buckle into the same 
wave length. The critical stress of the plate was deter- 
mined from the condition of neutral stability in which 
the work of applied forces must equal the sum of the 
strain energy of the buckled plate and the energy of the 
restraining medium, 


T=V,+V, (1) 


If w is the deflection normal to the plate at any point 
(x, y) in the plane of the plate shown in Fig. 1, 7’ and 
lV’, are given by the following equations [reference 1, 
Eqs. (201) and (199) }: 


t s* “| (~)’ ow? 
i = N,- N ( ) + 
si / Ox il Oy 
_ Ow Ow 
2Ney dx dy (2) 
Ox OV 
—— D f [ (oe ‘ omy 
= 2 0 JO l ox Oy" 


O’w O’w o'w \* 1) 
a1 = p@) e = oe dx dy (3) 
Ox? Oy" Ox OY f : 


The stiffness S) of a plate has been defined‘ as the ratio 
of moment per unit length at any point along the re- 


a 
ad l 7 n=1 


n 


(n + 4 = odd). 
the constraining relationship of Eq. (10). 


grangian function 
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strained edge to the rotation in quarter radians at that 
point when the moment is distributed sinusoidally, 
Since it was assumed that the restraining moment was 
in-phase with the rotation of the edge, plate stiffness § 
With this definition of S,, 
I’, is given by the following equation [reference 5, Eq 
(B-4) ]: 


—— I I) Jet 
= £0 ax 
0 ov \ ai 
[ | (2) | 
2 So dx } 
J oy \ hb 


The deflected surface may be represented by the 


is constant long the x axis. 


Fourier series 


TX nTy 


> a, sin (5 
b 


n=1 


w= sin 

The intensity of distributed forces acting in the middle 
plane of the plate is given by the equation 

N, = Na, [1 — Oy b) | (6 


Substitution of wand .V, from Eqs. (5) and (6) in Eqs. 
(2), (3), and (4) yields 


Nor’ 
a, T (1 = *) > Ay? 7 


a 2 
Noab > ,Qynt -) 
\ ig w=+<4) 
(n + 1 = odd), 
‘ iii 
V, = S0*A 9) 


where 6 is the amplitude of the sinusoidal edge rotation 
and is given by 


4= : B é. (N17 h) (10) 
n=1 


Substituting these expressions in Eq. (1) and making 
the substitutions 


No = k (x°D bd") 


Sy = De’'4b, B = X/b 


gives the following expression in k: 


Bb? 6? 
4 


= > a,2(1 + B'n?)? + 


sa 17)? 1 T 


n 


The coefficients a, must be so adjusted as to make the expression above for & a minimum with 
To minimize k with its constraining relationship, construct the La- 





£(d,,9) 


nx1 
follow! 


n= 1 
The 
= (), 
of the 
equati 


[o 


Ne 


In. 
14) t 
minal 
value 
desire 
stress 
calcul 
value 
and |} 


Cn 
lengt 
wher 
chose 
given 
simp! 





le 





























~RITICAL BENDING STRESS FOR FLAT PLATES 537 


28D?" SRB? ee: 1,Q,nt 
L(a,0) = Do an2(1 + B2n?)? + cee ss _ > s — kp? (1 — *) > a,? - 


n=] a4 T n=1 i (w* = -¢*)° -/ n=] 
nT 
Y (« — > a ) = 0 
aa b 
OL 16RB? a ant a nT 
| = 2a,(1 + Bn”)? — — > ———. — bf? (1 — (2a,) + 7 = (0) 
OU, 1" ; (nm? — 27)? 2 b (12) 
las 00 = (484b?6/r*) — y = O 
n+i= odd). The equations of expression (12) with the constraining relation (10) give, after simplification, the 


following system of linear homogeneous equations: 


Qa Ska 8° i 
\«, t + Bn?)? — ea (1 _ ) _ = > h- ~8 site nw (*) — 
2 " ; (nm? — 17)? 2 b 


ia ¥ 
nan — = 0) 
| x “ teB4 (;) | 
n= 1 = odd). 


The results presented in this paper were determined by neglecting all even a, greater than d¢ or aj,3 = 0, where 
j= 0, 2,4, 6,..., ©. The work of solving Eqs. (13) may be greatly reduced by evaluating the odd a’s in terms 
of the even a’s. Performing these operations and substituting in Eqs. (13) leads to the following set of stability 





equations: 
24,2\2 , a 64a°Btk?n : 272 
»| (1 + B2n?)? — RB? (1 ~~ *) a ; > - aie — - a — 
= Tv i 1,3,5 (m- — 1°-)° A |(7? _ })2 
64a2B*h2n > 17? 64a2Btk2n S ait 
4 aud o PIE ‘. eye i - — = — ios il 
ail 13,5 (m* — 2°)®[A] (2? — 16)* n* i-1.3,5 (nm? — 17)? [A] (2? — 36) 


nt SaB 7k 1° (*) 
l = 0 
» | > x? ay ; (n? — 72)? [A | b (14)* 


<3! Sap k i” 16B4 SaB?k 1? | 
1 + : . ly + [ 4+ 1 
T° wT” ;21.3.5 (17 — 4)? [A | a wy “ 5 (i? — 16)? [A] “ 


2484 Sap7k ' :? | | l 264 3 ¢? | ( 1 ) 
| so ieee + v = 0 
| . i 2 5 (17 — 36)? [A] . € ew ; a ; [A] b 


n = even), where 
[A] = {(1 + 6%?)? — kp? [1 — (a 2)]f 


* Note that Eqs. (14) reduce to Eqs. (Al1) of reference 2 by setting a = 2. 


In order for the system of linear homogeneous Eqs. unloaded compressive edge. The curves given in Figs. 
14) to have other than the trivial solution, the deter- 6 and 7 for intermediate values of a were obtained by 
minant of the coefficients must vanish. The minimum  cross-plotting the results calculated in the present paper 
value of k for which this determinant vanishes is the and those of references | and 2. 
desired critical stress coefficient. Values of critical 
stress coefficients k satisfying this condition have been EFFECTS OF RESTRAINT ALONG TENSION EDGE 
calculated and are given in Figs. 2, 3, 4, and 5 for all 


9 


4 : vig a The results determined for this paper are for a long, 
values of « and 8 and for values of @ of 7/4, 3/2, 5/4, pal 


flat rectangular plate under loading and restraint con 
ditions of Fig. 1. This case was considered because it 
was assumed that restraint along the tension edge 
would have little effect on plate stiffness. To check 

Critical stress coefficients for flat plates of’ finite the validity of this assumption, the critical stress co 
length may be obtained by substituting a/mb for \,b, — efficient for a plate simply supported along the com 
where the number of half-wave lengths m must be _ pression edge and built in along the tension edge was 
chosen so that k is a minimum. Stress coefficients are determined for various a/b ratios. These results were 
given in Figs. 6 and 7 for finite plates having all edges compared with critical stress coefficients for a plate 
simply supported and for finite plates having a fixed simply supported on all edges. This comparison 1s 


and |. 
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FIGURE 2. CRITICAL STRESS COEFFICIENTS FOR A FLAT PLATE IN 
BENDING IN THE PLANE OF THE PLATE WITH VARIOUS RESTRAINT 
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FIGURE 3. CRITICAL STRESS COEFFICIENTS FOR A FLAT PLATE IN 
BENDING IN THE PLANE OF THE PLATE WITH VARIOUS RESTRAINT 
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IGURE 4. CRITICAL STRESS COEFFICIENTS FOR A FLAT PLATE IN 


BENDING IN THE PLANE OF THE PLATE WITH VARIOUS RESTRAINT 
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FIGURE 5. CRITICAL STRESS COEFFICIENTS FOR A FLAT PLATE IN 
BENDING IN THE PLANE OF THE PLATE WITH VARIOUS RESTRAINT 
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FIGURE 6. CRITICAL STRESS COEFFICIENTS FOR A FLAT PLATE IN 
BENDING IN THE PLANE OF THE PLATE-ALL EOGES SIMPLY SUPPORTED 


shown in Fig. 8. For all a, it can be seen that restraint 
along the tension edge increases the plate stiffness and, 
The amount 
of this increase varies with the a/b ratios and a. An 
assumption of simple support along the tension edge 
gives conservative values for the critical stress coefficient 


thereby, the critical stress coefficient. 


for values of a> 1. The degree of conservatism is not 
so great as indicated in Fig. 8, because the condition 
of built-in edge cannot be attained in any practical 
structural design. Therefore, an assumption of simple 
support along the unloaded tension edge of a plate in 
bending would lead to results slightly conservative for 


practical designs. 


CONCLUSIONS 
1) The curves of Figs. 2, 3, 
cient data with the results of references 2, 6, and 7 for 
use in computing the critical bending stresses of beams 
by use of the moment distribution method of reference 
4, 
(2) The stress coefficients given in Figs. 6 and 7 should 
be used to compute the critical bending stresses of flat 
rectangular plates of finite length having the particular 


1, and 5 provide suffi- 


edge conditions shown. 

(3) In the practical case, the assumption of simple 
support along the unloaded tension edge of a plate in 
bending gives slightly conservative critical stresses. 
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FIGURE 7. CRITICAL STRESS COEFFICIENTS FOR A PLATE IN BENDING 
IN THE PLANE OF THE PLATE-UNLOADED EDGE ON COMPRESSION SIDE 
FIXED, ALL OTHER EDGES SIMPLY SUPPORTED 
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FIGURE 8. COMPARISON OF CRITICAL STRESS COEFFICIENTS 
FOR HINGED AND FIXED TENSION EOGE 
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The Aerodynamic Characteristics of Low 
Aspect Ratio Wing-Body Combinations 
in Steady Subsonic Flow 


H. R. LAWRENCE? 


Cornell Aeronautical Laboratory, Inc. 


SUMMARY 


The linearized flow field about a thin wing of arbitrary twist 
mounted on an infinite inclined circular cylinder in a uniform in 
compressible free stream is formulated as the solution of a pair of 
simultaneous integral equations in two variables 

By making use of an extended form of the slender-body ap 
proximation, the simultaneous integral equations are shown to be 
equivalent to a single integral equation valid on the wing surface 

The method employed by the writer in references 9 and 15 is 
used to approximate this two-variable integral equation by a 
single-variable integral equation for the chordwise lift distribu 
tion. The chordwise lift distribution on the wing-body configura- 
tion is shown to be identical with that of a wing with the same ex 
posed plan form but altered twist distribution. A procedure for 
computing the distribution of lifting pressures is presented 

Methods for extending this result to bodies of finite length and 
to account for the effects of compressibility are discussed. 


List OF SYMBOLS 


body radius 


{A = local 


semispan of wing-body configuration 
measured from body centerline 

= maximum semispan of wing-body configura 
ti0n 


= semichord of wing at wing-body intersection 


F(x, ¢ = @ + rv |see Eq 54 

N(x = [see Eq. (28 

g(x = {see Eq. (36 

H(¢ = |see Eq. (56 

Ux, 4 = lift per unit area on wing 

L(x = lift per unit chord on body 

Ly (x = lift per unit chord on wing 

L(x = lift per unit chord on wing-body configuration 
P(n) = {see Eq. (33) and Appendix | 

t = free-stream static pressure 

p = static pressure 

Ox, 9 = (see Eq. (31 

Ry = projection of wing on xy plane 

Rx = projection of body on xy plane 

r, 0 = polar coordinates in ys plane 

l = free-stream velocity 

u(x, y) = I(x, y)/2pU? 

Wal xX, ¥ aerodynamic angle of attack, see Eq. (49 
W(X, ¥ = geometric angle of attack, see Eq. (48) 
wi(X, 9 = induced angle of attack, see Eq. (50 
x,y, = cartesian coordinates 


Cartesian coordinates corresponding to x, y, 


c. = 


respective ly 
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angle of attack distribution on wing 


a(x, J 

a = angle of attack of body 

B(x) = h(x) a, local exposed semispan of wing 

By = b a, Maximum exposed semispan of wing 
Vv? = Laplacian operator 

oo (‘@ = velocity potentials 

m = see Eq. (53 

c =y+is 


INTRODUCTION 


he problem of calculating the aerodynamic charac 

teristics of a wing-body configuration was first con- 
sidered by Lennertz.' The case of a high aspect ratio 
wing mounted on an infinite circular cylinder was in- 
vestigated by making use of the Prandtl? lifting-line 
theory in conjunction with the concept of image vortices 
in the Trefftz plane. A formula was obtained for the lift 
on a wing with uniform loading, and the twist distribu- 
Mul 


thopp* investigated the case of a high aspect ratio wing 


tion for minimum induced drag was computed. 


mounted on an infinite cylinder of arbitrary cross sec- 
Following Lennertz, he assumed that the wing 
Multhopp then 


tion. 
can be considered as a lifting line. 
proved that the wing-body problem can be reduced to 
the solution of the Prandtl lifting-line integral equation 
by conformal transformation of the wing-body con- 
figuration in the Trefftz plane. Although the wing- 
body problem for cases where the Prandtl wing theory 
is valid is solved in principle by this method, analytic 
difficulties and laborious numerical calculations are en- 
countered in practical application of the Multhopp pro- 
cedure. Flax and Treanor‘ have developed a variational 
procedure that appears to avoid these computational 
problems. 

The aerodynamic theory of high aspect ratio wings 
mounted on infinite cylindrical bodies appears to be 
reasonably complete for cases where the Prandtl lifting- 
Development of a method for 
high 


line theory is valid. 
computing the aerodynamic characteristics of 
aspect ratio swept-wing-body configurations remains 
as a subject for further research. 

A method for computing the aerodynamic character- 
istics of low aspect ratio wing-body configurations was 
first presented by Spreiter.° The method is restricted 
to cases in which accelerations in the stream direction 
can be neglected in the Laplace equation. ‘This type of 


solution will be called the ‘‘very low aspect ratio’’ case, 
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since the same assumption, when applied to a wing 
alone, leads to the Jones® theory, valid for aspect ratios 
below one-half. Spreiter solves the wing-body problem 
by conformal transformation, thus reducing the problem 
to the Jones case. Although the transformation used is 
identical with that employed by Multhopp, the analysis 
proceeds with little difficulty because of the simple 
character of the Jones solution. Relatively simple 
analytic solutions are obtained for many cases of in- 
terest. Ward’ has derived the Spreiter results for the 
supersonic case by a more rigorous procedure based on 
the Laplace transform and has cleared up certain dis- 
crepancies in the Spreiter theory. Adams and Sears® 
have clarified the discrepancies of the Spreiter theory in 
the subsonic case by an adaptation of the Ward tech- 
nique in which the Fourier transform replaces the La- 
place transform. The objective of the present paper is 
to develop a method for computing the aerodynamic 
characteristics of wing-body configurations at small 
angles of attack in a subsonic flow field. The results will 
be restricted to low aspect ratio wings (say 0 to 4) with 
straight trailing edges and arbitrary twist mounted on 
finite bodies of circular cross section at arbitrary angles 
of incidence. The body radius will be assumed constant 
along the wing-body intersection. 

The problem will be attacked by considering the case 
of a wing mounted on an infinite circular cylinder in an 
incompressible flow. Slender-body theory may then be 
used to generalize the result to bodies of finite length, 
and the Prandtl-Glauert transformation may be applied 
to remove the restriction to incompressible flow. 


FORMULATION AS A BOUNDARY-VALUE PROBLEM 


Consider a thin wing with straight trailing edge 
mounted on an infinite circular cylinder embedded in a 
uniform incompressible free stream of velocity U/ in the 
positive x direction. The wing and body axes lie 
approximately in the plane z = O with angles of attack 
a(x, vy) and ap», respectively, relative to the free stream. 
The wing-body configuration is assumed to possess a 
plane of symmetry y = 0. 

The flow field about the infinite cylinder of radius a, 
when the wing is not present, may be described by the 
velocity potential ¢ defined below. 


d(x, r, 0) = Ux + Uae cos 6(a?/r) (1) 


where (7, @), the polar coordinates in the yz planes, are 


defined by the formulas 


s = recos8, y = rsin 0 (Z) 


“ 


The velocity potential ¢", which describes the flow 
field about the wing-body configuration, is considered 
to be the sum of the potential ¢ due to the body and a 
perturbation potential [7 due to the presence of the 
wing. The angles a(x, y) and ay are assumed to be small 
so that the approximations of linearized lifting-surface 


theory may be employed. 
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The conditions that (1) the potential function $") pe 
harmonic outside rigid bodies and their wakes, (2) the 
flow through rigid surfaces vanish, and (3) the effects oj 
the wings vanish far upstream may be stated in the 
following manner: 


A*7o@™ = 0 (in flow field) (3 
¢:"! = —Ua(x, vy) (on wing surface) (4) 
¢‘? = —Uaycosé@ (on body surface) (5 
@ = ¢ (atx = —o) (6 
Upon making use of the definition ¢” = ¢° + U¢ 


and Eq. (1) for ¢®, the above conditions may be re 


written in terms of ¢ as follows: 


A*g = 0 (in flow field) (7) 
a- , 

o. = —a(x, y) — a (on wing surface) (8) 
y 

¢ =0 (on body surface) (9) 

p = (0 (atx = —o) (10) 


Since only the diiferences in pressure between the 
upper and lower surfaces of the wing or body will be 
computed, the nonlinear terms in the Bernoulli equa- 
tion may be neglected. 

The linearized Bernoulli equation reduces to the 
following relationship between the static pressure p and 
the perturbation velocity potential l’¢: 


(in flow field) (11) 


bo — p = pl*gz 


where /» is defined as the static pressure in the free 
stream. 

The boundary problem is uniquely specified by the 
conditions (7), (8), (9), and (10) when supplemented by 
the Kutta condition. The forces acting on rigid sur 


faces may be computed from Eq. (11). 


FORMULATION AS AN INTEGRAL EQUATION 


Consider the function $(x, vy, z) defined in terms of the 
undetermined function u(x, y) by the formula 


| ~ £ dé dn sulé, n) 
o(x, ¥, 3) = = = Aas 
: 2a SeNF = +s" 


x—& 
t= (12 
Yo ~- +e - oF + a , 


where the integration is taken over the entire plane 
z= 0. Itis evident that ¢(x, y, 2) satisfies Eqs. (7) and 
(10). Upon differentiating Eq. (12) with respect to x 
and allowing z to approach zero from above and below, 
there results the formulas 


o:(x, y,0 +) = u(x, y) (13) 


(x, v, O—) = —u(x, y) (14) 


By making use of Eq. (11), the lift per unit area 
I(x, y) acting on a thin wing located in the plane z = 0 


may be expressed in the form 


I(x, vy) = pU*[ox(x, v, O +) — o2(x, y, O—)] (15) 
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LOW ASPECT RATIO WING-BODY COMBINATIONS 543 
per unit area acting on a thin wing in the plane z = 0. 
Let Ry and R, be defined as the projections of the wing 
Since the 


or, with the aid of Eqs. (15) and (14), 


I(x, y) = 2pU7u(x, y) (16) and body, respectively, on the plane z = 0. 
lift and, hence, u(x, y) must vanish at all points off 
Rw + Rz in the plane z = 0, the region of integration in 


Eq. (16) shows that the undetermined function 
Eq. (12) can be at most Rw + Re. 


(x, y) may be interpreted as proportional to the lift 


[he boundary condition (8) may be reformulated by differentiating Eq. (12) with respect to z and, after some 


manipulation, allowing z to approach zero. The final result of this process is the integral equation 


a? 1 oO - i dé dn u(é, n) V (2 — &)° + (vy — @)° 
a(x, y) + a = = * 1+ = (17) 
. 27 OY . Ry + R, I-99 x—é€ 


y 


valid on the wing surface Ry. 
12) may be rewritten in polar coordinates defined by Eqs. (2) as shown, 


reos? f ff dé dn ulé, n) x—eé 
d(x, 7,0) = ~~ / / - - -11+ : | (18) 
2T JSR, +R, r* — 2ra sin @ + n° aS (pm t)? + r? — 2rn sin 6 + 7? 


Upon differentiating with respect to 7, simplifying the resultant expression by partial integration, and substitut- 


Eq. 


ing r = a, there results the formula 


cos 8 0 a a dé dn y(&, n)(a* — 7?) x—é 
(x, a, 0) = 5 If * = cs " : 1 + s = 
24 OY, Ry +R, <0" — “Yn + a’) Vv (x g*) + 9° — 2yq + * 


cos 6 O -_ @ dé dyn u(é, n) x—¢€ 
> . (19) 
2r Oy Ry + Rp =n V(x — &)? + 7° 2yn + a? 


where y is defined as a sin @. An integral equation valid on the region R, is obtained from the boundary condition 
9) by setting ¢,(x, a, @) equal to zero in Eq. (19). The integral equations (17) and (19) represent a complete 
formulation of the wing-body problem when combined with the Kutta condition and Eq. (11) for the pressure. 


integral equation formulation of the wing-body problem, 


In order to simplify the analysis, an approximate 
which will be employed in subsequent developments. 


integral equation valid on Ry will be derived to replace 
Eq. (19). It will be assumed that the radical occurring 
in the integrands of Eq. (19) may be replaced by 
x—&. This is equivalent to the assumption that the 
wing-body configuration is slender insofar as satisfying 


THE CHORDWISE LIFT DISTRIBUTION 


The lift per unit area on the wing /(x, y) is expressed 
acd : am in terms of the function u(x, y) by Eq. (16). The lift 
the boundary condition on the body is concerned. The : ; j ; 
; ; Me per unit chord Lyw(x) acting on the wing is therefore 
final result, however, is not restricted to slender con 
At large values of » where the inequality 
no longer holds, the integrand ap 


: : given by the formula 
figurations. : 

5 "dix 
z=~—~ & ae n = 
i ' ; : Lw(x) = 4pU” dy u (x, ¥) (21) 
proaches zero as 1/n, and all wing-body configurations ‘ 








appear to be slender with respect to the boundary con- 
This conclusion is restricted, of 
The approxima- 


dition on the body. 
course, to the case of slender bodies. 
tion employed will be further justified by comparison 
with known results valid for high aspect ratio wings in 
combination with cylindrical bodies. 

Upon defining $(x, vy) as the integral from — ~ to x of 
u(g, y), the approximate integral equation valid on the 


region R,z may be written in the form 


re) *b(x) dn o(x, n)(a*? — 7?) 
0 = ne Bast DA - (20) 
Oy biz) 20(9? — 2yn + @*) 
where }(x) is the wing semispan when 6(x) > a and 


b(x) = a otherwise. Eqs. (17) and (20) represent the 


where the integration region has been halved by con- 
sideration of symmetry about the plane y = 0. By 
using the Eq. (11), the lift per unit chord acting on the 
body Lz(x) may be written in the form 


* 


Lax) = pl 3 J d6 a cos 6 (x, a, 8) (22) 
0 


The function ¢,(x, a, 0) may be evaluated in terms of 
u(x, y) by making use of Eq. (18). As in the case of the 
boundary condition on the body surface, the radical in 
Eq. (18) will be approximated by |x — &. Then, upon 
differentiating Eq. (18) with respect to x, substituting 
a for r, and integrating over the variable £, there results 


the formula 
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acos@ {x : 
z(X, a, 0) = s - (23) 
T _b(x) @* — 2an sin 6 + 7? 


Eq. (25) for $,(x, a, 6) is now substituted in Eq. (22), 


and the order of integration 
result 


pU-a? rou 
L(x) = i dyn u( 
T b(x) 


AUGUST, 1953 
Fourier series. The following formula is obtained: 


2s do’ sin? 6’ T 
> 7, = Te oe 
0 a — 2an cos 0’ + 7° a 


T 
=— (i9| > a) 
” 


dn u (x, 7) 


is interchanged to give the 


x, 1) X and substituting in Eq. (24) the lift per unit chord act- 


ing on the body may be written in the form 


2n d6 cos? 6 ™ 
s. © sin 0 2 (24) 
o a 2an sin 6 + 7 


The integration over 8 can be performed by substitut- 
ing 0’ + (m 2) for 6 and expanding the integrand in a 


a 
La(x) = 2pU? Lf dn u(x, n) + 
J0 
*hix) a’ ’ 
dn — u(x, n) (25 
a 


DEVELOPMENT OF A SINGLE-VARIABLE INTEGRAL EQUATION 


The first step in developing a single-variable integral equation will be to reduce the simultaneous integral 


equations (17) and (20) to a single integral equation. It will be shown that by suitably relating the values of 


(x, y) on Rg and Ry. Eq. (20) is satisfied identically. The method employed is closely related to the use of 


Eq. (20) is rewritten in the form 


ra) | “a dn o(x, n)(a? — v) my a . *b dn o(x, n)(a? — 4 a 
oy _q 2n(n? — 2yn + a?) b a 2n(n? — 2yn + a?) 


The variable 7’ = a*/n is introduced in the second term to give the result 
dn (x, a*,/n) 


oO “a dn o(x, n)(a® — n°) *+a?/b re 2 
| 7 $( n v) 4 . 2 
Oy _g 2n(n? — 2yn + a?) we a2/b 29(n”? — 2yn + a?) 


It is evident that the above equation will be satisfied identically if the following relations hold 


images. 


d(x, vy) = O (\y| < a*/b) ( 
d(x, a?/y) = (x, y) (a2/b < jy < b)§ 


(26 


Eqs. (26) may be used to reduce the region of integration of Eq. (17) from Rw + Rz to Ry alone. There results 


the integral equation shown below in which the region of integration and the region of validity are identical. 


a” 1 oO ~ £ dé dn u(é, n) V (wx — £)? + (vy — n)? 
+ a(x, ¥) = en [ + i + 
¥ . 2m Oy Ru ~—y x—& 


1o0 {f F£ dédn us, n) a y + V (x — &)? + [y — (a? n)]*) = 
2a ot | J. y — (a*/n) 2° | x—é , = 


In order to obtain a single-variable integral equation, Eq. (27) will be multiplied by a weighting factor and in- 
tegrated over the wing span.’ Following reference 10, a solution for the case of very low aspect ratio, 
V [b + (a*/b))? — Ly + (a*/y)]? 


will be chosen as the appropriate weighting factor for the case on hand. Upon multiplying (27) by this weighting 


factor and making use of the definitions 


Pe ”b aya” a*\ a7\* 
f(x) = / +f dy a Ca I | y(0 + ) (s + ) (2S) 
=" a ‘ e b \ y 


*¢ P5 *b dn u(&, n) V (x £)? + (y — n)* 
Q(x, ¥) = / a f + | yo [ + ; | + 
b a lial = 2 
s *—a “b> dn ult, n) a? 
dé + oe 21 
Eee b Fe 3 Clee | 


there results the relationship 
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| *—¢ b a*\? a*\* Oo ; 
f(x) = oe : + | dy y(2 — 4 ae (» + a iy Q(x, ¥) (30) 


where c is the wing semichord at the intersection of the wing and body. 

In evaluating the integral expression (29) for Q(x, y), the radicals will be approximated by the expression 

(1 2)[|Iv — & + Vix — &)? + (6 — a)?] 

It is evident that when |x — £ is large compared with the span, this approximation is satisfactory. The limiting 
case in which the aspect ratio approaches infinity will be discussed later. It will be shown that, for a rectangular 
wing, the integral equation to be derived by this approximation approaches a proper limiting form. 

When the approximate value for the radicals in Eq. (29) is introduced and the spanwise symmetry is taken into 
account, the formula for Q(x, vy) may be written in the form 


"4 d x —& + Vix — &)° + (6 — a)? y a* y 
O(x, vy) = 2 dé dn u (&, 7) [ te = aa 2 ‘ 
. ™_ se 2(x — §) fag fF C/G 


The integral equation (30) will now be transformed by integrating the right-hand side by parts with respect to 


y. Since Q(x, +a) vanishes, as seen from Eq. (31), and the radical vanishes at y = +4, there results the form 
. 1 f*®dvQ(x, y) (y* — a‘) 
f(x) = - (32) 
wT Sa ¥ V [b + (a?/b)?] — [vy + (a?/y)?] 
Upon substituting Eq. (31) for Q in Eq. (32), interchanging the order of integration, and making use of the 


definition 


“b dy yi - ai 4 a’ y 
P(n) = 7 ee ae os oe (33) 
a YY V [(0+ (a?/d)? — [y+ (a? yy)? ] LX -— er WY — (aD) 


Eq. (32) may be written as follows 


l * x- + V(x — &)°+ (b -— a) 7 
f(x) = / dé |2 + : | / dn u(é, n)P(m) (34) 
T se — - ; 


The integral expression (33) for P() is evaluated in the Appendix. There results the formula 


Str 


P(n) = (w/2)[1 + (a?/n?)] (35) 


It is convenient to define the function g(x) as follows: 


ms *b a” 
g(x) = 2 / dé / dn u(&, n) (1 + :) (36) 
‘ ( ee n~ 


The integral equation (34) may be written in terms of g(x) in the form 


| l “ p V (x — &)? + (6 — a)? 
f(c) = — p(x) + dé g’(—) | 1 + : 
2 Pie, = -§ 


The function g’(x) may be interpreted physically in terms of the lift per unit chord of the wing-body configuration. 


37) 


The conditions of Eq. (26) are employed to transform the expression (25) for the body lift per unit chord to the 


“bis a? 
L(x) = ou? f dn — u(x, 7) 
a n- 


Upon adding the above formula for Ly to Eq. (21) for the wing lift per unit chord, the total lift per unit chord L(x) 


°(b)s a 
L(x) = 4pU* / dyn u(x, n) (1 + ) (38) 
Ja i 


form 


is given by the expression 


and making use of Eq. (56) 


L(x) = 2pt/*e"(x) (39) 
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The Kutta condition is applied as a spanwise average by the condition that g’(x) vanish at the trailing edge of the 
wing. The integral equation (37) together with the Kutta condition g’(c) = 0 and Eq. (28) for the function f(x) 
are sufficient to determine the function g(x) when the plan form and angle of attack of the wing-body configuration 
are specified. The lift per unit chord is then determined from Eq. (39). 


COMPARISON WITH OTHER THEORIES FOR LIMITING CASES 


When the aspect ratio of the exposed wing is small, the integral equation (37) is seen to reduce to the form 


*b(x) oft? 2\2 2\ 2 
g(x) = f(x) = 2 | dy ac y) + = - | y(° + ) — (y + _ ) (40) 
A y’ b y 


This result is identical with the solution obtained by Flax! for vanishing aspect ratio. The special case when 
ay = Ois given by Puckett and Morikawa." If the angle of attack of the wing a(x, vy) and body ap are equal, the 
well-known result of Spreiter® and Ward’ is obtained. 

The limiting form of the integral equation (37) as the span becomes large will be evaluated. In order to simplify 
the analysis, a rectangular untwisted wing will be specified. Furthermore, the wing and body will be assumed to 
be at the same angle of attack. The integral defining the function f(«) may be evaluated explicitly in this case 


> 


so that the integral equation (37) may be written as follows: 


z\2 |] 1 V (x — £)? + de? [1 — (a de) |? 
= beta (: ) = 9 g(x) + | / dé g’(€) )! + ; : | 41) 


bo? x—é 


where is the semispan of the configuration. When the span becomes large, the above integral equation reduces 


l © dé g’(&) a \* 
— g : = 2 ao( bo — a) (: > ) (42) 
TJo.%—€ bo 


Eq. (42) is the well-known integral equation of two-dimensional airfoil theory. The solution for g(¢), when the 


to the limiting form 


Kutta condition g’(c) = 0 is specified, is given in reference 12 as 
g(c) = 2nr[1 + (a/be) |?ao(bo — a)e 43) 
and the lift coefficient based on exposed wing area is given by the formula 
Cy = 2rao{1 + (a/bo)]° (44) 


This result agrees with the results of reference 4 to within + per cent for a, bp ratios up to 0-4. 


SOLUTION OF THE SINGLE-VARIABLE INTEGRAL EQUATION 


The solution of the integral equation (37) will be shown to depend in a simple manner upon the solution of the 
problem for a twisted wing without a body. A method for solving this wing-alone case has been presented in 
reference 9, and hence, most of the numerical work required for the solution of the wing-body problem is available. 

The exposed semispan of the wing will be denoted by 6(v) = b(x) — a. The integral equation (37) may be 


| | “6 d V (x — &)? + B(x) 
f(x) = 5 &(X) + , / dé g’(—) |} 1 + 5 (45) 
2 , x gE 


It is evident that Eq. (45) for g(x) is identical with the integral equation (43) of reference 9 when the known 
function k(x) defined therein is replaced by f(x), defined by Eq. (28). Furthermore, f(x) reduces to k(x) as the 
body radius approaches zero. Hence, the integral equation (45) may be said to contain the integral equation (43) 


rewritten in the form 


of reference 9 as a special case. 
For the case of a flat wing at angle of attack ae equal to the body angle of attack, the function f(x) may be writ 


ten as follows: 
f(x) = (#/2)aob?[1 — (a?/b*)]? (46) 


where 6 is, in general, a function of x. In the special case of a rectangular wing, f(x) is a constant, and the solution 


for the wing-alone and the wing-body problems diifers only by a constant factor. 
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CALCULATION OF THE PRESSURE DISTRIBUTION 


The distribution of lifting pressures may be computed by applying the procedure developed in reference 15. The 


integral equation (27) is rewritten in the form 


20 “b y a’ y = 
w(x, ¥) = het 4+ (40) 
w OY Ja y? — » 7” y? — (a‘t/n*) 


where the following definitions are used: 
W(x, ¥) = a(x, vy) + (a°ao/¥") (48) 


Wal(X, ¥) = w(x, V) — wil(X, V) (49) 


1 oO . ae V(x — £)?+ (vy — nn)? — [x —€& l a*/n* 
w(x, ¥) = - dé dn u(é, n) * ; ° + = (50) 
2a Ov R x— yn y — (a?/n) 


Ww 


twtr 


If the aerodynamic angle of attack w,(x, y) were known, Eq. (47) would be the governing integral equation for 
a very low aspect ratio wing. The solution of this integral equation is well known. (See, for example, reference 
11.) 

The aerodynamic angle of attack we(x, y) may be approximated in the following manner: It is evident from Eq. 
(50) that the induced angle of attack w,(x, y) vanishes as the aspect ratio approaches zero and represents a small 
correction to the geometric angle w,(x, y) for wings of low aspect ratio. Following reference 15, the induced angle 
of attack will be assumed to be constant over the span or w;(x, y) = w;(x) only. It will be noted that this assump 
tion is analogous to assuming that the downwash is constant over the wing chord in the Prandtl* lifting-line theory. 


Eq. (47) may be rewritten in the form 


20 *b y a’ y : 
w(x, Vv) — w(x) = dn o(x, n) | ~~ te, ier gar 4 (51) 
wT OY Ja = <5 9° y° — (c*/n’) 


In order to compute the induced angle w;(x), Eq. (51) is multiplied by 
V(b + (a2/b) 2 — [y + (a2/y)?? 


and integrated over the wing semispan. Upon performing certain integrations and making use of Eqs. (28) and 


(36), the following formula is obtained: 
w(x) = [f(x) — g(x)]/(m/2)b°u(a/b) (52) 


where uia ) is given by the expression 


| ” a*\* a*\? 
nb? J ’ V(° . i) 7 (» ) 
14 as 2a (: ) | (1 4 ~) — } 2( a/b) | eal 
= = = 7 _ Ss (3d 
b4 rb b? T bh? 1 + (a?/b?) 


The aerodynamic angle of attack may thus be expressed by Eqs. (48), (49), (52), and (53) in terms of the geometry 
of the configuration and the solution g(x) to the integral equation (37). 
The pressure distribution may be computed from Eq. (11) when the potential function ¢ is known. @ may be 


written as the real part of the complex potential /, or 


F(x, 6) =o@+ 1, c=yt+2z 4) 
The complex potential function F may be calculated from the formula (see, for example, reference 11) 
ee. 1 P ;, H(g) + iia) 7” 
"x, ¢) = aN WalX, n) log 2 (>) 
wT Ja IT(¢) - TT(n) 


where // is defined as follows: 


l(c) = V [b + (a*/b))? — [& + (a?/5)}? (56) 





548 JOURNAL OF THE ABRONAUTICAL SCIENCES 


RESULTS AND EXTENSIONS 


The primary result of this paper is a method for cal- 
culating the symmetric lift distribution on a configura- 
tion consisting of a wing with straight trailing edge and 
arbitrary twist mounted on an infinite circular cylinder 
at arbitrary angle of attack in steady uniform linearized 
incompressible flow. The following engineering prob- 
lems can be solved by suitable modifications of this 
result: 

(1) The three-dimensional Prandtl-Glauert trans- 
formation may be used to compute the lift distribution 
at high Mach Numbers. In the case of low aspect ratio 
thin-wing slender-body configurations, the results ap- 
pear to be applicable almost up to Mach Number unity. 
The calculation is so nearly identical to the well-known 
wing-alone case that it will not be presented. 

(2) The restriction to infinitely long bodies of con- 
stant radius may be partially removed by a slight ex- 
tension of the theory. Consider the case where the 
body radius is constant along the wing-body intersec- 
tion. If the body forward of the wing-body intersection 
is slender, the lift distribution of this forebody may be 
computed by slender-body theory and added to the re- 
sult for an infinite body-wing configuration without 
loss in accuracy. For the case of a slender body of 
variable radius aft of the wing trailing edge, the situa- 
tion is somewhat more complex because of the inter 
action of the wing wake and body. This case has been 
analyzed by Adams and Sears.* The resulting lift dis 
tribution may be added to the lift distribution for the 
forebody and the infinite body-wing configuration. 

(3) The effects of flaps and camber on the chordwise 
lift distribution may be computed in exactly the same 
way as for a wing alone. The wing-alone problem has 
been discussed in detail in reference 9. 

(4) The theory developed in this paper is strictly 
applicable to low aspect ratio wings. However, as 
noted previously, the predictions of this theory are 
satisfactory for all aspect ratios in the case of a rec 
tangular wing-body configuration. The assumptions 
made in the present analysis are essentially the same as 
those employed in reference 9 for the wing-alone case. 
Since the predictions of reference 9 for wing alone are 
found to agree with test data for aspect ratios from 0 
to 4 at all taper ratios, there is good reason to believe 
that a similar situation will apply to the wing-body 
theory presented herein. 


Appendix 


EVALUATION OF P(n) 
The integral expression for P(n) given in Eq. (33) 
may be evaluated by substituting a new variable of 
integration /, where 


(y=t+ VPP —-a?’) 
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and making use of the definitions 
2s = + (a?/n) [2r = B+ (a?/B)] 
to reduce the integral to the form 


25 "7 dit t- — a- 
P(n) = / a. % i : 
n a jek 


The substitution of a new variable zs = Vf? — a? + 
Vr? — f reduces Pin) to the form 


Pp 25 ] l l 
(n) = c - 
” n | , ] he s- ¢ — =) 


s*>— a? 


The second term of the integral vanishes when inter 
preted as a Cauchy principal value, whereas the first 
term is easily integrated to give the result 


P(n) w 2)[1 + (a?/n?)] 
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A Constant-Area, Constant-Temperature 
Combustor’ 


PHILIP ROSEN? ano R. W. HART? 
The John Hopkins Unwersity 


ABSTRACT 


The theoretical performance of a piloted constant-area, con 
stant-temperature combustor with mass addition is studied 
Perpendicular injection of mass at constant injection velocity 
and temperature is assumed. Mixing between the injected mass 
and the main stream is regarded as instantaneous, and the flow 
The effect of chemical kinetics 
The exit ve 


is regarded as one-dimensional 
is considered by the use of a simple rate equation. 
locity corresponding to maximum throughput is found to be 
V RT. The combustor length required to accelerate the gases to 
this velocity is found to depend on the reaction rate constants, 
flame temperature, input flow variables, and the temperature of 
combustion. It is determined that there exists a lower limit on 
the temperature below which the combustor is unstable. 


I) INTRODUCTION 


Sami RAM-JET COMBUSTION CHAMBERS consisted 
basically of a fuel injector and a “‘flameholder”’ 
at the upstream end of a cylindrical pipe. A fuel mix 
ture was injected into the air stream and ignited. The 
flameholder——an obstacle to airflow such as a strut(s 
or dise(s)—-provided eddy regions where motion of the 
gases was slow enough to permit rather complete burn 
ing and consequent pilot action. It was found that 
these combustion chambers are stable over a range of 
fuel and flow variables but that the stability range was 
exceedingly difficult to predict and depended on fuel 
injector geometry and flameholder geometry. 

Most research in ram-jet combustion problems has 
been concerned with a search for empirical relations 
between flow, fuel, and geometry variables, because of 
the difficulty in characterizing theoretically compressible 
turbulent flow; the effect of flameholder and injector 
geometry on the flow; the effect of mixing and heat 
transfer; and the che*nical kinetics of the reactions. 
Several investigators made simplified theoretical 
analyses of cylindrical pipe operation.'~* One-dimen 
sional flow and addition of heat were assumed. By 
using the equations of conservation of mass, energy, 
and momentum with the equation of state of an ideal 
gas, useful information, particularly limiting velocities 


and mass flows, was obtained. The effects of fuel prop 
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erties on operations of flameholder-type ram-jets 
which enter the problem through kinetics of reaction 
or mixing and geometry effects have apparently not been 
In the idealized cylindrical pipe 
theory is 


treated theoretically. 
combustors, where the above-mentioned 
valid, heat transfer and mixing are perfect and instan 
taneous throughout any cross-sectional element, but no 
transfer takes place along any direction parallel to the 
evlinder axis (ef. reference 5 

In this report we will consider the cylindrical com 
bustor from a similar point of view but will introduce 
the effect of fuel type by considering the rate of heat 
production to be described by a simple reaction rate 
equation of arbitrary order and by assuming that the 
temperature in the combustor is a fixed constant. A 
similar approach to a constant-pressure constant 
temperature can-type combustor has been used by 
Avery and Hart.’ An additional simplification is mad¢ 
by assuming that the combustor receives a known mass 
flow of stoichiometric composition from a pilot and 
through the cylindrical walls. With these assumptions 
the rate of heat production given by a rate equation 1s 
then required to equal the heat rate required (a) to warm 
the injected mixture up to the combustor operating 
temperature and (b) to increase the velocity of the in 
jected mixture, as well as the velocity of the original 
gases. Application of the conservation equations then 
vields relationships that must exist between rate of 
addition of unburned fuel to the burner and rate of 
combustion to determine the burner length and range of 


stability. 


(II) THEORETICAL ANALYSIS 
(1) The Conservation Equations..-The differential 
equations expressing conservation of energy, momen 
tum, and mass are :* 
Neglecting heat losses, diffusion, and burner drag, 
the energy conservation is 
dQ Vv? — VI, = r dm mdV? 
dx 2 ~ ax 2 dx 
rhe conservation of momentum equation is 
—A(dp/dx) = (d/dx) (pA V*) 


The equation of conservation of mass is 


(d/dx) (pAV) = dm/dx 
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where 


dO/dx = 


rate of production of heat (per unit length) 
caused by chemical reaction 


m = 

A = cross-section area of the cylindrical com- 
bustor 

P = static pressure of gases in the combustor 

ce = average heat capacity of fuel-air mixture 
and combusted products 

!% = absolute temperature = constant 

ia = velocity of injected mixture, assumed per- 
pendicular to the cylinder walls and con- 
stant 

3 = temperature of injected mixture 


With the constant values of I’,, 7°,, and A, all of the 
conservation equations can be integrated, 


. E =z Ve" 
m 
2 


where C is the constant of integration, 


(la) 


+ C,(T - r.) | +C=@0 


P+ pV? = Po + poVo? = constant = F_ (lb) 

(The subscript zero indicates evaluation at the up 
stream end of the combustor.) 

m= pAV (Ic) 


These three equations and the equation of state P = 
pRT constitute four equations relating five unknowns 
(O, P, p, V, and m). The temperature, 7, and the 
area, A, are specified. The rate equation will relate 
heat production, Q, to the other four variables and 
therefore permit us to determine the behavior of P, 7, 
p, V,m, and Q. 

(2) Limiting Exit Flow Velocity. 
of a cylindrical combustor with heat addition and with- 


Just as in the case 


out mass addition, an important limiting flow condition 
can be determined without considering details of the 
chemical reaction. It will now be shown that the con- 
stant-temperature combustor has maximum ‘“‘through- 
put’’+ when the exit Mach Number is (y)~ 

Differentiation of Eq. (le) with respect to x, and 
elimination of the density by means of Eq. (1b) and the 
equation of state yield a relation between the rate of 
mass addition and the rate of change of velocity, 

dm  AF(RT — V*)dV_— AFQ — yA) dV (20) 


dx (RT + V2)2 dx RT (1 + yM?)2 dx 


where y¥ is the specific heat ratio and .J/ is Mach Number 
of the flow. 
A second differentiation yields 


dm 


_ AF (zy 
dx? |yavar 2(RT)*\dx / |\y-vnz 


t Total mass flow through the combustor exit 


(2b) 
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rate of mass flow through the cylinder at 
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The constant F is greater than zero, so that the mag 
RT or M? = | y3 
But it is easy 
to show that the requirement that heat be generated 
(i.e., that dQ dx be positive) requires that dV’ /dy by 


flow rate has a maximum at J” = 
(unless dJ” dx is zero with d?V//dx? < 0). 


positive for I” less than R7. Using Eq. (la) with 
AFI (RT + IJ”) following from Eqs. (1b) and 


t 


nm = 
(lc), we find 


dQ dm sail ae re ; 
O< = é(f — 7.) +,(V? — V,*) is 


dx dx 2 
dV AFYV? 
dx RT+ 12 
and using dm, dx trom Eq. (2a), 
dQ diy AFyA AFy 1 — 7M 


O< xX 


dx dx la + yl?) y—10+ 7M)? 


Now, evaluating at 7? = RT, 
dQ AFdV 
O< = 
dX y= R7 2 dx | y:=R1 
The condition that the rate of heat generation dQ dv 
be greater than zero requires that d1’/dx be positive 


according to Eq. (2d). 


(The region |” greater than 
RT is of no significance in the investigation of this 
combustor, since it has now been shown that di/di 
becomes negative in this region, corresponding to mass 
subtraction.) It is easily shown, further, that dV’ di 
is positive for all velocities less than (R7)'” when 
heat is generated in the combustor. From Eq. (2d), 
this condition obviously can immediately determine 
that, for all velocities less than this critical velocity, 
the bracketed term is positive. The operating tem 
perature of the combustor (7°) is necessarily greater 
than the temperature of the injected fuel-air mixture 
(7), and the Mach Number with which the fuel-air 
mixture enters the combustion chamber (.J/,,) is, in 


any case, less than sonic, so that whenever the term 


(1 — y.Z*) is positive the bracketed term is also posi 
tive. In the limiting case when there is no unburned 


fuel beyond a certain point in the combustor, dQ ‘dx 
will be zero and so dJ"/dx will be zero. 

It is easily shown that the ratio of mass throughput 
to mass flow input from the pilot increases with de- 


creasing input Mach Number. From Eqs. (1b) and 


(Ie), 
ni M E +1)? My 1+ yM? 
mo MoLRT + V? MoL1 + yM 
t A similar phenomenon occurs at the same Mach Number in 
the analysis of a constant-area combustor given by Chambre and 


“= 


Linn (reference 1, p. 537 
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Since this ratio increases with increasing |’, we obtain 
maximum throughput by inserting the velocity, V = 


V RT, and 
m mo = (1 + Mo?)/2V y MI 


exit 


(2f) 


3) Combustion Efficiency..-We will now investigate 
the combustion efficiency 7, which the conservation 
equations require to exist in this constant-temperature 
combustor without regard to the type of heat-pro 
ducing reaction. 

Let 7 equal mass flow rate of completely combusted 
iuel-air mixture divided by total mass flow rate. Then 
the heat flow rate Q is given in terms of the heat of 
combustion by 

QO = » Ii 


where /7 is the energy released per unit mass by com 


plete combustion of the injected fuel-air mixture. 


Nitrogen in the fuel air mixture is to be included in the 
combustion products. ) 

Eq. (la) specifies the rate at which heat must be pro 
duced in order to maintain the combustor temperature 


constant. Rewriting this equation, we find 
C CAT «7, mT 
+ + (3a) 
mil Hf] 2H 


Itis evident that, downstream from the pilot where the 
mass flow rate is much greater than the input mass 
flow from the pilot, the term C,m tends toward zero, 
and the combustion efficiency tends to become inde 
pendent of the integration constant C. The combus- 
tion efficiency approaches as a limit the value 


C((T-T.) , V2 — V,' 


y 


| 
H 2H 


(= — =") - ] 
= (3b) 
r 3 ve ts 2 


The downstream limit combustion efficiency, 7*, cor 
responding to maximum throughput is obtained by 


substituting |” RT into this equation, and letting 


I] Clr ti) 
pa hate, Tomety(y— 1) 9 
T, ge T - 24 

which is relatively insensitive to changes in the tem 
perature (7°,) and Mach Number (.1/,) of the injected 
mixture, provided 7°, is much less than 7 and that the 
injection Mach Number is small. | Note that the flame 
temperature, 7’,, is here defined by C,(7 — 7.) = H = 
heat of combustion.] Since part of the heat results in 
accelerating the gas, the burner temperature can never 
reach T,. The integration constant C is determined 
by the composition, temperature, and velocity of the 
pilot output. It is easily shown that the amount of 
fuel that can be burned in a combustor operating at 
maximum throughput corresponding to the pilot tem- 
perature 7, m = FA/2~\/RT, is independent of the 
choice of C, although the combustor length does de- 


We calculate (7..i.) (nexit) — Mono = Mass 
Substituting the 
expression for efficiency derived above, and V? = R7, 
we find that the mass burned per second depends on 


pend on C. 
burned per second in the combustor. 


the pilot temperature but does not depend on the choice 

of C. The required combustor length to attain V? = 

RT does, however, depend on the choice of C. The 

burner analysis is simplified by selecting a pilot for 

which C = 0. This can be effected, for example, by 

regarding the pilot as an extension of the combustor. 
With C = 0, 


Cy. Vv? — V,? 
=: fo T.54 (3d) 
‘* 2 TH 
If» = 1, an upper limit )’* is set on the velocity, 


V* = V2H + V, 


2 — 2C,(T — Ty) 


V 2C,(T; — T) + V.? (Ge) 
The highest temperature (7*) possible with maximum 


exit velocity 1’ = V RT is found from Eq. (8e) setting 


V* = V RT,* 
T* = [2y/(3y — 1)]7, (3f) 


when J/, < 1. 


(4) Kinetics of Reaction... We now express the rate 
of heat production per unit length, dQ dx, in terms of 
the state and flow variables; then, in the following sec- 
tion, the conservation equations are solved for these 
variables, with the aid of an assumed .Vth-order reac 
tion. 

Assume a simple Nth-order reaction 
rate of change of concentration of fuel (C,) due to chem- 


i.e., the time 


ical reaction 1s 


OC, Ot Oy 4 heel * 


The concentration of fuel is readily expressed in terms 
of the air-to-fuel ratio (a/f) and the concentration of 


unreacted fuel-air mixture by 


p(1 — »)/[1 + (a/f) 


The air-to-fuel ratio will be assumed stoichiometric so 
that (a/f) is constant throughout the combustor. 

The amount of heat produced per second per unit 
length is proportional to the mass of fuel-air mixture 


reacted per unit length—4.e., 


a OC; a 
All [ 5 | 
of f 


so that under steady-state conditions 


1 — (1 — ») |}. ’ 
: (nlIm) = AK’e P af 4 1+ (a/f)|H 
dx 1 + (a/f) 


“= ( nm) = 
dx dx — 


(4a) 


This differential equation is rewritten as a differential 
equation in velocity by substituting 7 from Eq. (8a) and 


AFV/(RT + V?*). 


m from Section (2), m = 


5 5 l 





552 JOURNAL OF 


‘= fj 
4 (F + 


THE 


he 
21 RT + V2) J 


AERONAUTICAL SCIENCES 


AUGUST, 19:53 


K'’e- = “(1 2 | F yf ¢. - r ig eee Vy? \ 
— ( ian y= 
F f RT + 1? iT 21 (4b 


Eq. (+b) can be rewritten in simpler form by using Eq. (8e). 


[V4 + (8RT — 2H + V**) V? + RT(V*? — 2H)] (V? + RT) 


(v2 — ye] 
where A = K’/[1 + (a/f)]*~"'. 


A solution of this equation gives the velocity as a function of position in the combustor for an .Vth-order reac. 


We obtain 


2 ( — | \ 
iV< — KF*-1¢ 
(2H)*-! 


tion. Other flow and state variables are expressed in terms of |’, 


m = AFV/(RT + V’°), 


P = FRI) (RT + V?), 


(35) Solution of the Velocity Differential Equation for a Second-Order Reaction (N = 2). 


readily integrated. We find 


n= [C,(T — 


ep = F/(RT + V’?) 
T,)/H) + ((V? — V,,?)/2H] (4d 


bee og 4 fe + eee VM cs. 
\ 0 y*2 = V2 V2 oe Vo? ( 


y*2 





2V*44 (2RT — H)V* — = fe 4 yr ~ Ve i 
- mmf a m= &(V, Ve V..T.T..T, = 
| 2)* "Whee — V/A + 4 sath 


where x — Xo = length of combustor. 


Eq. (5a) indicates the behavior of velocity as a function of position. 
ever, to visualize readily, although it is clear that the limiting velocity |’* (corresponding to 7 = 1) is reached only 
It has already been shown that (1*)? > RT, dv/dx > 0. 


at infinite distances. 


1S? [Po + po Vo" Je mie 4 
(X — Xo) (Oa 
2H [1 + (a/f) 


The relationship is too complicated, how 


For C # 0, we obtain the following expression, which is also integrable, although somewhat more complicated: 


ny (U8 BRT — 2H + V8) + REV ADV AV ee 


Ji. (45, - vie + 


AF 


Ye 


2RCT \;* = F (x — Xp 
2TT 


AF 


In this study we concern ourselves only with the simpler case C = 0 corresponding to 


no = [(Vo? — V7) 277] + [((T — 7,)/(T; -— T,)] 


Fig. 1 shows a plot of log, (® e”®") against 7. 
Curves A and B are for E = 33,000 cal. per mole and 
Mo = O and 1/10, respectively, where Myo = Vo + 
(yRT)'*. Curves C and D are for E = 25,000 cal. per 


mole and for J) = 0 and 1,10, respectively. These 


curves show what one would expect —namely, the 
length required is larger for higher activation energies 
and slightly larger for smaller incident Mach Numbers 
The effect of Jo is not great if 1) < 
1. Consequently, for large throughput, m 


for constant 7. 
exit./Mo, the 
length is relatively insensitive to increases of Mteyi¢_/M0. 

In the above and following calculations, we have used 
values for the constants as indicated below: 


R = 2/30 = gas constant = 1/15 cal. per Gm. 
°K. for air > 2 cal. per mole °K. 


y = 1. 
T, = 2,400°K. 


C, = 0.288889 cal. per Gm. °K. 
Py) = 1/5atmosphere 

T, = 480°K. 

a/f = 16 


K’ = (a/f)B’ 
B’ = 5X 10" ce. per Gm.sec. 
V. = Vo 


(6) \Winimum Operating Temperature and Stability. 
Before we study further calculations, it is illuminating 
to consider the temperature range within which the 
We have al- 
ready seen that the maximum attainable temperature 


combustor is capable of stable operation. 


is given (for small /)) by 7* = [2y (8y — 1)]T;, and 


With V = 2, Eq. (4c) is 
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that this temperature corresponds to an exit velocity 
V* which is only attained at infinite length where the 
combustion becomes complete [Eq. (5a)]. A minimum 
stable operating temperature must exist, because the 
rate of chemical reaction decreases rapidly (essentially 
exponentially) with decreasing temperature, whereas 
the heat required to warm and accelerate the incoming 
mixture decreases relatively slowly (nearly linearly for 
the range of Mach Numbers considered). This means 
that below a critical temperature a disturbance that 
causes a drop in temperature may cause the heat gen- 
erated to fall below the heat required to maintain energy 
balance and the temperature continues to fall. 

It is easily shown that this combustor is unstable to 
changes in pilot temperature for temperatures less than 
that corresponding to the minimum of the curve of 
Fig. |. This is so because, to the left of these mini- 
mums, a drop in the combustor temperature increases 
the length required to maintain that temperature and, 
since the length is fixed, the temperature must fall fur- 
ther. To the right of the minimums, a slight drop in 
temperature decreases the length required, and this will 
tend to bring the temperature back to what it was be- 
fore. Similarly, a slight rise in temperature will in- 
crease the length required, and this will also tend to 
bring the temperature back. Thus, at first thought, we 
could consider our lower limit of temperature as corre- 
sponding to the value just beyond the minimum to the 
right. There is, however, a more severe stability cri- 
terion that yields somewhat higher minimum temper- 
ature. 

Let us consider a perturbation in which an extra 
amount of mass is added at some point in the com- 


bustor. This will lower the local temperature sud- 
denly. If this additional mass does not burn, then 
d(dnm/dx) = 0. If it does burn, then 6(dym/dx) > 0. 


Thus, from Eq. (4a), for a second-order reaction 
5[p2(1 — n)2?e 87] > 0 
If we consider that pressure and velocity remain 
(roughly) constant, then, 
bp = (p/7T)67, bn = 677/(7, — T, 


and, thus, 


or, since 67° < O, 


2 2 E 


ee = — ———- >9Q (6a) 
7 (1 — n) (T, — T,) 


If we choose 7 = no then we find that 


ae i 


| + (2RT;,/E) 
If 7, = 2,400 and FE = 33,000 cal. per mole, then 
T 2 1,860°K. 2 0.775T;,. 
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T° KELVIN x 107° 
Fic. 1. The log of the length factor de”/*" is plotted against 
temperature; Curves A and B, E = 33,000 cal. per mole, 1/) = 
0 and 0.1, respectively; Curves C and D, E = 25,000 cal. per 
mole, J, = 0 and 0.1, respectively 
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Fic. 2. The Mach Number and velocity are plotted against 
distance in the combustor; Curves A and B are Mach Number at 
T = 1,860° K. and 2,152°K., respectively; Curves C and D are 
the corresponding velocities 


If n is closer to unity, then the minimum temperature 
is below this value, as can easily be seen from Eq. (6a). 
This is as one would expect, since, if the burning gas 
were unstable at the pilot, it would be unstable further 
downstream. This lower limit was obtained graphi- 
cally by Avery and Hart for an activation energy of 
FE = 33,000 cal. per mole. 

(7) Results.—Let us consider E = 33,000 cal. per deg. 
mole and My = 1/10 and calculate the length of the 
burner at the lower limit of the stability region 7 
1,860°K. At this lower limit, we find L = 26.5 cm. 
At the upper limit 7 = 2,152°K., the length is infinite. 

In Fig. 2 Curves A and B are plots of Mach Number 
in the burner vs. x, when Py = 1/5 atmosphere, for 
T = 1,869°K. and 2,152°K., respectively. We see that 
for 7 = 1,860, the Mach Number increases monotoni- 
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Fic. 3. The ratios of mass flow rate to initial flow rate and of 
pressure to initial pressure are plotted against distance in the 


combustor; Curves A and B are m/m for T = 1,860°K. and 
2,152°K., respectively; Curves C and D are P/P) for T = 
1,860°K. and 2,152°K., respectively 
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Fic. 4. The combustion efficiency, total mass flow rate per 
unit area, and the mass flow rate of burned gases per unit area 
are plotted against distance; Curves A and B are efficiency for 
T = 1,860°K. and 2,152°K., respectively; Curves C and D are 
corresponding curves for #/A and Curves E and F are nm/A for 
the above temperatures, respectively. 


cally until it reaches M = (1/y)”° atx = 26.5 em., 
while for 7 = 2,152°K. the Mach Number approaches 
(1 y): , 
Fig. 2 are the actual velocities in feet per second. 
Curve C is for 7 = 1,860°K. and Curve D is for 7 = 
2,152°K. At T = 1,860°K. the velocity starts at 270 
ft. per sec. and rises to 2,400 ft. per sec., while for 7° = 


asymptotically at infinity. Also plotted in 


2,152°K. the velocity starts at 295 ft. per sec. and ap- 


proaches 2,575 ft. per sec. at infinity. The velocity 


difference between these two cases at the exit is rela 
tively small considering the difference in length. 

In Fig. 3 Curves A and B are plots of the ratio 
m/mo vs. x at Po = 1/5 atmosphere for T = 1,S60°K. 
and 2,152°K., respectively. On Curve A, 7/71 reaches 
a maximuin of 4.44; at Curve B mm approaches the 
Curves C and D are 
plots of the ratio of pressure at any station to entrance 


same value asymptotically at ©. 


pressure vs. x at 7) = 1,860°K. and 2,152°K., respec- 
tively. On Curve C, P/Po reaches a value of 0.505, 
and on Curve D P/ Po» approaches the same value at 
infinity. 

In Fig. 4+, Curves A and B are plots of 7 vs. x for P. 
1/5 atmosphere for 7) = 1,860°K. and 2,152°K., re 
spectively. On Curve A, 7 reaches the value of 0.83 from 
a value of 0.72 at the pilot. On Curve B, 7 approaches 
the value unity at infinity from a value of 0.S7 at the 


pilot. 

Curves C and D are plots of mass per unit area per 
unit time in pounds per foot second, for P» 1/5 
atmosphere for 7 = 1,860° and 2,152°K., respec- 


tively. On Curve C, /A increases from 0.66 Ibs. per 
sq.ft. sec. to a definite maximum of 2.9 Ibs. per sq.ft. 
sec. On Curve D, /A approaches a value of 2.7 
Ibs. per sq.ft sec. from an initial value of 0.61 Ibs. per 
sq.ft. sec. Here again we would gain nothing with a 
higher temperature. 

Curves E and F are plots of » mA for T = 1,860°K. 
and 2,152°K., respectively, at Po = 1/5 atmosphere. 
On Curve E, 7 m/A starts from a value of 0.47 Ibs. per 
sq.ft. sec. and increases to a value of 2.4 Ibs. per sq.ft. 
sec. On Curve F, nm/A approaches 2.7 Ibs. per sq.ft. 
sec. at infinity from a value of 0.58 Ibs. per sq.ft. sec. 
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Studies of Some Effects of Airplane 
Configuration on the Response 
to Longitudinal Control in 
Landing Approaches 


RALPH W. STONE, Jr.,* ann WILLIAM BIHRLE, Jr.* 
Langley Aeronautical Laboratory, NACA 


SUMMARY 


The investigation reported in this paper was undertaken to see 
whether there are inherent differences in the response of the flight 
path angle to control movements between swept-wing airplanes 
having no horizontal tail and using trailing-edge flaps for longi- 
tudinal control and conventional airplanes that are known to 
have satisfactory response characteristics. Particular emphasis 
is placed on the final few seconds of the landing approach 

This study shows that the only difference in the short-time re 
sponse of the flight path between the two types of airplanes is in 
1 time lag in the flight-path angle response of the swept-wing air 
planes. This time lag was found to result from two factors: 
first, a relatively small value of the parameter, 


(VCt,/2uc) ( V2Cm, 2uk,?) Ad, 


which is a measure of the effectiveness of the control in changing 
the flight-path angle and, second, a relatively large value of the 


parameter 
( VOL, 2uc )A6, 


which is a measure of the amount of undesirable change in lift 
accompanying the control deflection required tochange the angle of 
attack. It had been suspected that the relatively large change in 
drag with angle of attack, of low aspect ratio swept-wing air 
planes having no horizontal tails, would be an important factor 
in the flight-path angle response, but this was found not to be the 
case for the short time intervals under consideration 

The importance of the differences found in the response char 
acteristics between the two types of airplanes studied can only 
be evaluated by flight experience. Other factors such as range 
of vision, control feel, the pilot’s reaction to the relatively large 
nose-up attitudes of the low aspect ratio swept-wing airplanes, 
and psychological influences associated with new type airplanes 
may be of equal or greater importance. 


INTRODUCTION 


WwW" THE DEMAND FOR HIGH-SPEED FLIGHT, air- 
plane configurations have changed appreciably, 
and necessary compromises have led to some flight diffi- 
culties particularly at low speeds. One such difficulty 
was reported in reference 1 for free-flying airplane 
models having low aspect ratio wings and no horizontal 
tails. The models responded uncertainly to elevator 
control movements, particularly when flying beyond the 


Presented at the Aerodynamics Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953 
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angle of attack of maximum L/D ratio. It was felt 
that troubles may be encountered on similar airplanes 
at low speeds, particularly in landing approaches where 
small height corrections may be needed. A theoretical 
investigation of the responses to longitudinal control 
of airplanes having low aspect ratio wings and no hori- 
zontal tails, therefore, was undertaken at the NACA’s 
Langley Aeronautical Laboratory to evaluate the ef- 
fects of aerodynamic differences between this type of 
airplane and other more conventional airplanes that 
are known to have satisfactory response characteristics 
in landing approaches. 

Other factors such as the range of vision, control feel, 
pilots’ reaction to the relatively large nose-up attitude 
of low aspect swept-wing airplanes and psychological 
influences associated with new type airplanes also may 
have considerable influence on pilots’ opinions regard- 
ing the landing approach characteristics of such air- 
planes. The theoretical investigation, however, was 
made without regard to any such possible psycho- 
logical influences. 

Motions in response to longitudinal control move- 
ments were computed, on an analog computer, for 
two swept-wing airplanes having no horizontal tails 
and for a conventional airplane that was known to have 
good height control in landing approaches. 


SYMBOLS AND COEFFICIENTS 


The longitudinal motions presented herein were calculated 
g I 


about the stability axes. 


S = wing area, sq.ft 

( = mean aerodynamic chord, ft. 

Ww = weight of airplane 

m = mass of airplane, slugs (lV’/g 

ky = radius of gyration about }-body axis, ft 
p = air density, slugs per cu.ft 

mm = airplane relative-density coeflicient, m /pS< 
H = velocity, ft. per sec 

g = acceleration due to gravity, 32.2 ft./sec.? 
L = lift, lbs 

D = drag, lbs 

VW = pitching moment, ft.lbs 

CL = lift coefficient, L/(1/2)pV2S 

Cp = drag coefficient, D/(1/2)pV2S 

Cn = pitching-moment coefficient, 17/(1/2)p V2.8; 





JOURNAL OF 


wr 
Q 
oo 


THE AERONAUTICAL SCIENCES— 


AUGUST, 1953 


TABLE 1 


Mass, Dimensional, and Aerodynamic Parameters (Aerodynamic Parameters for Airplanes B and C Generally Given Relative 
to Those for Airplane A) 


W ky Cing Ciné, 
Airplane S é (Cm), (Cmse), 
A 19.1 0.99 
B 12.7 0.57 0.12 0.29 
Cc 26.1 0.46 0.04 0.18 
Cy. = lift coefficient at a = 0° with an elevator deflection 
that would be required to trim at landing approach 
speed 
t, = pitching-moment coefficient at a = 0° with an ele- 


vator deflection that would be required to trim at 


landing approach speed 


T = thrust, lbs. 
Zz = height, ft. ( / V sin y dt 
70 
a = angle of attack, deg. (a = 0 — y) 
Y = flight-path angle, deg. 
7] = angle of pitch, deg. 
6 = elevator deflection, deg. 
AZ = increment of height from trimmed condition 
Aa = increment of angle of attack from trimmed condition 


= increment of flight-path angle from trimmed condi- 


tion 
A6 = increment of angle of pitch from trimmed condition 
AV = increment of velocity from trimmed condition 
Aé, = increment of elevator deflection from trimmed 
condition 
éorg = pitching angular velocity, rad. per sec 
V = rate of change of velocity V with time 
CL, = OC,/06,, per deg 
Co, = OCp/06,, per deg 
Cn, = OCm/06,, per deg. 
Ci, = OCL/da, per deg. 
Cp «= coefficient of drag as a nonlinear function of a 
"$e = = OC,,/Oa, per deg 
a = CC,,/%qc/2V 


AIRPLANES INVESTIGATED 


The airplanes for which motion calculations were 
made are referred to herein as airplanes A, B, and C. 
Airplane A is a conventional straight-wing airplane 
that has good response characteristics in the landing 
approach. Airplane B is a swept-wing airplane having 
no horizontal tail, and airplane C, which also has no 
horizontal tail, is generally similar to airplane B but 
has a lower aspect ratio and different mass character- 
istics. Pertinent aerodynamic, mass, and dimensional 
information on airplanes A, B, and C is given in 
Table 1. 


PROCEDURE 


For this investigation, the three longitudinal equa 
tions of motion were used in the analog computer: 
: v2 
V = —- —(C 


2 ue 


. Co; 6, ) = 


a 


gsin y + cos a_ (ia) 
m 


Crs, Cr6, Ca Cla Cirin 
(Crg,), (Cos.), (Cma), (Cra), hie. 
0.638 
Pee 1.61 0.65 0.62 0.79 
1.50 1.70 0.39 0.57 0.65 
Pe ee , 
_ = NOp. -F Lge + Chg. 83): = 
2ué 
g cos y + sina (lb 
m 
ia = ; 
= (Cm a@ + Cmo + Cm 6 : 
») k 2 a 
2uR, 
l" 
Cote l¢ 
tuk, : 


The equations were not linearized in that the thrust 
and gravitational forces were introduced in their non- 
linear form. The lift, drag, and pitching-moment co 
efficients were introduced as functions of angle of attack 
and of elevator deflection. All of the airplanes were 
trimmed in level flight at 110 knots (185.8 ft. per sec.) 
at the beginning of the motions, and the slopes of the 
lift and pitching-moment coefficients with angle of at 
tack were taken in the vicinity of the angle of attack 
for trimmed flight at 110 knots. The drag coefficient 
[Cpy:a)] was introduced as a nonlinear function of angle 
of attack. The lift, drag, and pitching-moment co 
efficients, as functions of elevator deflection, were intro 
duced as linear functions of 6,; the slopes were taken as 
mean values between the elevator deflection for trim 
and the maximum up-elevator deflection. The ele 
vator deflections and therefore the values of Cz; 4,, 
Cp,,6,, and Cms 5. were introduced as step functions 
The thrust and the airplane damping derivative, C,,,, 
were held constant for any given calculation. 

In order to evaluate the response to longitudinal con 
trol for the three airplanes, the following procedure 
was used: As previously mentioned, the three airplanes 
were initially trimmed in steady level flight at a landing 
approach speed of 110 knots. The initial trim values 
are given in Table 2. A disturbance from this cond! 
tion was initiated by deflecting the elevator down and 
holding it down for a short time (2 sec.), after which an 
attempt to stop the ensuing descent was made by de 
flecting the elevator full-up. The amount of down 
elevator deflection used initially for airplane A resulted 
in a loss of altitude which might be desired for a final 
height correction during an approach, approximately 
10 ft. 
the initial pushover flight paths the same for airplanes 
B and C as for airplane A, so that when the final up 
elevator deflection was used in an attempt to stop the 
descent all three airplanes would be in approximately 
This was attempted by de 


For comparison purposes, it was desired to make 


the same flight condition. 
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RESPONSE TO LONGITUDINAL 


flecting the elevators down on airplanes B and C by an 
amount that would result in approximately the same 
initial rate of change of normal acceleration for all three 
airplanes. The amount that the elevator was deflected 
down from the initial trim deflection for all three air- 
planes is given in Table 3. As previously mentioned, 
the elevator on each airplane was deflected full-up 
after the initial pushover motion. The maximum 
full-up elevator was used so that it would give the 
maximum acceleration in pitch possible on each air- 
plane, and any other correction by use of up-elevator 
would be proportional to that given by full-up elevator. 
It is realized that the use of full-up elevator and of 
step input functions does not simulate the actual condi- 
tions used by a pilot, but it was felt that the method 
does reveal any difference in response due to inhcrent 
differences in stability and control characteristics 
which exist among airplanes A, B, and C. 
that the airplanes would not trim at an angle of attack 
beyond maximum lift, the elevator was returned after 
a given time from full-up to a deflection that would 


In order 


trim the airplanes at maximum lift. 


SIMPLIFIED ANALYTICAL CONSIDERATIONS 


The calculations performed in this paper, as previ- 
ously noted, are based on the three longitudinal equa- 
tions of motion and are accurate to the extent of the 
completeness and accuracy of the aerodynamic data 
used. Short-period oscillation in longitudinal motions 
may be studied, within required accuracy, by consider- 
ation of only two degrees of freedom, assuming the ve- 
locity to be constant. It is possible, therefore, that 
factors pertinent to short-time responses to longitudinal 
be obtained con- 


from 


CONTROL IN 
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TABLE 2 
Initial Trim Values for Steady Level Flight at 185.8 Ft. per Sec 
(110 Knots) 


Air- a, 6, 
plane Deg y, Deg Deg. T, lbs. be, Deg 
A $40 0 140 2.642 5 i 5° up_ 
‘ stabilizer 
B 21.85 0 21.85 1,770 —)() 0 
ts 15.97 0 15.97 2.781 ait = up : 
trimmer 


the gravitational and thrust forces in Eq. (1b) are of 
only secondary importance, and, if neglected, the 
equations when simplified to changes from the initial 
trimmed conditions may be solved simultaneously for 
the flight-path angle y, a pertinent parameter to changes 
in height, when the change in height may be expressed 
as 


AZ = / I sin y dt (2 


The solution for y gives the following result: 


VC 5, a ho2e™! — ye 
pod ) as, + 
2ue \b- 6 Va? — 4b 


( VCi :) ( aa) 
2ué 2uk,’ b? Va? — 40 
t a Vc, ICs, 
mi Aé, — ( = X 
b b 2ue 2uk,* 


» Ail » Ast 
Ao“e = A, "e 


control movement also may 

- ; ‘ wie ei where 

sideration of only two degrees of freedom. These de- 

grees of freedom are expressed in the equations of the VC VC.é 
. ° ° ° > . - G 
pitching and normal accelerations given previously as a * (= = q ) 
Eqs. (1b) and (Ic). It can be shown that changes in 2ue tuk, 

TABLE 3 


Elevator Deflections Used and Comparison of Resulting Parameters Affecting Longitudinal Motion (Parameters for Airplanes B and 
C Given Relative to Those for Airplane A) 


Pull-Up : 
Push-Down Aé,, Deg (a V2Cms. Ade VCLe (b UCrs. As 
Aé,, Deg (From : 2Quk,? Qué Que 
(From Trim Push-Down Aé. a VC rc 
Elevator Elevator (A6-), ( V2Cmbe Ade ¥ */) ( C Lie A8 ) 
(for Pull-Up) Quk,? Que ‘ Jul 4 ‘ 


Airplane Deflection ) Deflection ) 


Original Elevator Deflection 


0 


0 S84 


5S 0.25 
» 


25 1.14 
0.64 2.37 


Increased Up-Elevator Deflection 


1 86 —24 86 

B (c 1.33 — 14.33 
(c 2 45 —)(). 95 

B (c $35 (c) — 57.81 
& (c 2.45 (c) —32.71 


» 


99 

ono 
> 

° 


1.00 
1.00 


Parameter indicates total available elevator effectiveness in causing a rate of change of flight-path angle (neglecting change in 


lift due to elevator deflection 
>) Parameter indicates change in lift due to elevator deflection 


These Ade values resulted in the initial increment of y being approximately t 


due to elevator deflection 


he same as for airplane A (neglecting change in lift 
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If the velocity were constant as assumed in this sim- 
plified solution, y as expressed in Eq. (3) is of prime im- 
portance to an understanding of factors pertinent to 
the short-time longitudinal responses of airplanes. [f 
the lift due to elevator deflection, Ci s,. Were zero, then 
y, all its derivatives, and any height changes would be 
proportional to the factor 


( VC. 2uc) ( 1? Cm, 2uk,*) Ad, 


The down-elevator deflections required for airplanes B 
and C to obtain the same rate of change of normal ac- 
celeration in the initial pushover, previously discussed, 
were based on making this parameter the same as for 
airplane A. 

Another factor that may have important influence 
on the motion is the change in lift due to elevator de- 
flection (Cz;,), where an up-elevator deflection causes 
an imuiediate loss in lift proportional to the deflection. 
It should be noted that for a short coupled airplane such 
as airplanes B and C, which have no horizontal tails, a 
much larger change in lift is required in order to get a 
given pitching moment than is required on a conven- 
tional airplane with a horizontal tail. 


RESULTS AND DISCUSSION 


Comparison of Airplanes A, B, and C 


A fundamental comparison of the results of all three 
airplanes is given in Fig. 1. Here the maximum avail- 
able full-up elevator was used on each airplane to stop 
the descent following the initial pushover. In Fig. la 
is shown the change in height obtained and the elevator 
contro] movements used. The initial flight paths in 
the pushover were about the same as was desired, al- 
though airplanes B and C required about two and one- 
third and one and one-third times as much increment 
in down-elevator, respectively, as did airplane A to ac- 
complish this (see Table 3). Corrective action by the 
use of full-available up-elevator was attempted after 
the elevator was held down for 2 sec. The amount of 
the increment of up-elevator available on airplane B 
was a little more than one-half of that of airplane A, 
whereas the amount of the increment of up-elevator 
available on airplane C was about three-fourths of that 
of airplane A. After a short period of time, the ele- 
vators of airplanes A and C were moved down so as not 
to trim above maximum lift, whereas full-up elevator 
on airplane B trimmed the airplane at about maximum 
lift. 








All 
flecti 
heigh 
mucl 
or to 
as it 
abou 
furth 
plane 
ol st 
short 
were 
incre 
sent 
such 
resp! 
quic 

Fi 
ang] 
mar 
hit ; 
tack 
as d 
load 
atio 
it is 
pilo 
the 

pite 
as t 
the 
cha 
for 
pla: 

ol 

trir 

ang 
pla 
abc 

a | 

ter 


Eff 
] 


in 
air 
ot - 
pa 


pr 
tol 
eff 
ch 
ch 


pa 


fac 








In- 
m- 
to 


en 


be 








RESPONSE TO LONGITUDINAL CONTROL IN LANDING APPROACHES 550 


All three airplanes responded to the up-elevator de- 
flection in that the descent was stopped and the lost 
height was regained. It took approximately twice as 
much time, however, for airplane B to stop its descent 
or to regain its lost height after up-elevator movement 
as it did for airplane A. In addition, airplane B lost 
about 50 per cent more height and traveled some 1S0 ft. 
further before the descent was stopped than did air- 
plane A. It is possible that a pilot may become aware 
of such differences when making height correction in 
short-time periods. Further, presuming the airplanes 
were retrimmed at the bottom of the descents, the 6-ft. 
increment in height between them could represent as 
sent as much as, say, 100 ft. in the touchdown point; 
such an increment might be of concern. Airplane C 
responded more quickly than did airplane B but not so 
quickly as airplane A. 

Fig. lb shows the variations of angle of attack and 
angle of pitch for the motions just discussed. Of pri- 
mary importance in changing height is, of course, the 
lift and how fast it can be changed. The angle of at- 
tack of airplane A changes at a rate about twice as fast 
as does airplane B. Because of the difference in wing 
loading and C,,, airplane A changed its normal acceler- 
ation about three times as fast as did airplane B, and 
it is this acceleration that the pilot feels. What the 
pilot does is to pitch the airplane, and what he sees is 
the motion of the horizon. The rate of change of the 
pitch angle for airplane A is also about twice as great 
as that of airplane B. These differences account for 
the differences in response as measured by the height 
changes previously discussed. Here, also, the results 
for airplane C fall between those for the other two air- 
planes. Although not shown here, because the values 
of angle of pitch presented are increments from the 
trimmed values, the initial trimmed pitch or attitude 
angle of airplane A was about 4°, whereas that for air 
plane B was well over 20° and that for airplane C was 
about 16 Such differences in attitudes may influence 
a pilot's feelings with regard to an airplane’s charac- 


teristics. 


Effect of Increasing Total Elevator Effectiveness 


In order to demonstrate reasons for these differences 
in height response, the total elevator effectiveness of 
airplanes B and C was increased to approximately that 
of airplane A by making the total elevator effectiveness 


parameter, 


(Crq/2ue) (Cmz /2uk,?) Ad, 


previously discussed, the same for all airplanes. The 
total elevator effectiveness as used herein refers to the 
effectiveness of the elevator in causing an initial rate of 
change of flight-path angle, neglecting the effect of the 
change in lift due to elevator deflection. Making this 
parameter the same is done arbitrarily by assuming the 
factor to be increased by increasing the full-up-elevator 
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deflections, based on the same slopes of the pitching 
moment with elevator deflection previously used. The 
maximum full-up-elevator deflection of airplane B 
was increased from —30° to about —74°, whereas for 
airplane C the deflection was increased from —20° to 
about —30°. It was realized that the increase in ele- 
vator deflection on airplane B was excessively large 
and beyond the linear range of Cm, 

The results of calculations with the total elevator 
effectiveness, the same for all three airplanes, are shown 
on Fig. 2. Here there is a great deal of improvement in 
the response of airplane B, Fig. 2a, the time for the air 
plane to respond to up-elevator deflection is now only 
slightly larger than that of airplane A, and the incre- 
ment in height between them has been reduced to only 
about 2 ft. There was only slight improvement in air- 
plane C, but, of course, less improvement and less cor 
rection was required. 

The variations of the angle of attack and angle of 
pitch are shown in Fig. 2b. Here the rate of change of 
angle of attack, as indicated by its slope with time, is 
more nearly the same for airplane B as for airplane A 
than before. Similarly, airplane B pitched much more 
rapidly than before when only its normal full-up ele 
vator was used. The changes on airplane C are only 
slight in that the required change in elevator deflection 


was much less extreme. 
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Effect of Eliminating the Lift Due to Elevator Deflection 


The lift due to elevator deflection is a force in oppost 
tion to the desired force. This force is larger on short 
coupled airplanes, such as airplanes B and C, than on 
conventional airplanes, as was previously discussed. 
It was reasoned that eliminating this factor may show 
some improvement on the height responses of airplanes 
B and C. 


neglected and with the total elevator effectiveness in- 


The results of calculations with this force 


creased in amounts as previously discussed are shown in 
Fig. 3. 

The height changes and control deflections are shown 
on Fig. 3a. The response of airplanes B and C now 
compare favorably with those for airplane A. The im- 
provement due to neglecting the lift due to elevator de- 
flection was greatest for airplane C. 

The time to respond to the up-elevator movement 


that is, the time to stop the descent —are all about the 
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same, although the actual changes in height are some. 
what different. 

The variations in the angle of attack and pitch angle 
with time for airplanes B and C, Fig. 3b, are also similar 
to those of airplane A, accounting at least in part for the 
Thus, the 
differences in short-time responses to elevator contro] 


similarity in height response just discussed. 


for two widely different types of airplanes (airplane A 
representing conventional straight-wing types and air 
planes B and C representing swept-wing types having no 
horizontal tails) are primarily the results of differences 
in total elevator effectiveness and the increment of lift 
due to elevator deflection. 


Effects of Other Parameters 


It has been reasoned that other parameters than those 
previously discussed also may influence the short-time 
responses of the airplanes. It was felt that the great 
differences in the L/D ratios of airplane A and _air- 
planes B and C could affect the results in that airplanes 
B and C would tend to trim in much steeper glide path 
angles following the up-elevator movement (see refer- 
A comparison of the 1 D ratios of airplanes 
Airplanes B and C 


have much smaller L/D ratios for the landing condi 


ence 1). 
A, B, and C is shown on Fig. 4. 


tions, but of even more importance the variations of 
L/D with angle of attack from the trimmed angle of 
attack are much larger than for airplane A. The re- 
sults previously discussed, however, indicate that these 
differences in L/D ratio do not influence the short-time 
responses, primarily because there was no appreciable 
The 
effect of the increased drag causing the larger variations 


variation in the velocity during the short times. 


in LD on airplanes B and C is to reduce the velocity 
more rapidly on airplanes B and C than on airplane A. 
This is primarily a long-period effect, however, and 
may be expected to affect only the changes in height 
after long periods. Fig. 5 shows the relatively long 
time variations of height for airplanes A, B, and C with 
the airplanes in the same conditions previously shown 
in Fig. 5. The total elevator effectiveness of airplanes 
B and C was increased, and the lift due to elevator de- 
flection was eliminated for airplanes B and C for the 
calculations of Fig. 5, as it was in Fig. 3. The change 
in height for the first 4 sec. in Fig. 5 is the same as that 
shown on Fig. 3. From this time until the end of the 
motion shown, there is developed a great difference in 
the height for the three airplanes, airplane A gaining 
the most altitude. All airplanes appear to be in a long 
period or phugoid motion and, at the end of the motion 
shown, are descending because of the increased drag at 
Airplane C, hav- 
ing the least LD ratio, is descending most rapidly. 


the angle of attack of maximum lift. 


Other parameters that may be considered as im- 
portant to the motion are the stability derivatives C,,, 
and Cryo. 
damping derivative, C,,,, on airplane B was studied, 


The effects of increasing the value of the 


Increasing the 


but the results are not shown here. 
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TABLE 4 


Comparison of Parameters Affecting Longitudinal Motion (Parameters for Airplanes B and C Given Relative to Those for Airplane A 


l (a) V2C ms, I b) VCs, Zum (c V2Cms ‘Cu 
2uk,” 2uk,? Que 2 ye 2 2uk,,? 2 ui 
( I ) Coma, (S5:*) ( I ) Chie ( VCLs, Cla VCLe ViCms. VCte 
Airplanes Quk,?] 4 (Cmie)a Quky? Ja 2ue/ 4 (Cia) 4 2ye ), (Cie), ( 2 ui ), ( Quk,? Wl ) 
B 2.06 0.29 0.60 1.15 1.71 1.97 0.62 071 43 
Cc 3.94 0.18 0.71 1.88 1.50 2.82 0.57 1.07 0.76 


Parameter indicates elevator effectiveness in causing a rate of change of pitch angle 


Parameter indicates change in lift due to elevator deflection 
Parameter indicates elevator effectiveness in causing a rate of change of flight-path angle (neglecting change in lift due to ele- 


yator deflection 





ness) for airplanes B and C was about one-quarter and 
two-thirds, respectively, of that of airplane A. This 
indicates why the impractically large increase in up 
deflection, previously referred to, was re- 


value of C,,, did not appreciably influence the time at 
which the descent, following the pushover and then 
pull-up, was checked. The primary effect was to re- 
duce the rate of change of the pitch angle and angle of — elevator 
attack, and thus in the pushover and subsequent pull-up — quired for airplane B. 
the height changes were reduced, with increasing Cm: The parameter VCr,,A6,/2ué (Table 3), of course, 
No specific calculations were made to determine the. ,. , , 
i , ; mm : : a indicates the change in normal acceleration due to con 
effect of changes in the static stability derivative, C,,_. . : pe , 
, «trol deflection. Considering the up-elevator deflections 


From the simplified analytical considerations, Pre- available, this parameter is about the same for airplanes 
viously discussed, however, the parameter in which 4 and B. The value for airplane C was about two 
Cn, is dominate (assuming the lift due to elevator de- and one-third times as great as for either of the other 
flection, C.;, = O) is airplanes. 


In general, the time it takes an airplane of a specific 


Cn,/Sp7hy”) + (V?Cn,/2uky*)] : : : : ; 
configuration to respond in height, for a given elevator 


q 


“¢ La 
The values of these parameters for airplanes A, B, and deflection, depends on the rate at which it can change 
C are nearly the same. The motions for all three air- jts flight-path angle. It has been shown that, for a 
planes were also nearly the same when the totalelevator — chort -period motion, the rate of changing the flight- 

effectiveness of airplanes B and C was increased to that path angle depends mainly on the magnitude of the 
of airplane A and the lift due to control deflection, C15, | weight, moment of inertia (/,), Cm, Cra, CLs,, and Aé,. 
was eliminated. It may be reasoned, however, that re- An increase in weight, J,, or Crs,, or a decrease in Cy, 
ducing the static stability, C,,,, would allow an air- . — . , age 
reais é Siig fiat : Cm; , or Ad, will tend to decrease the airplane's abilitv 
plane to pitch faster, to increase its angle of attack , 

: ; to quickly respond in height. 

faster, and thus to respond more quickly to elevator ; 

deflection. Increasing C,,, would, of course, have the 











opposite effect. CONCLUDING REMARKS 
General Discussion The investigation reported in this paper was under- 
| etd : tl , Ived i taken to see whether there were inherent differences in 
?xaminé of some of the parameters involved 1n . : 
| pony — " : : ; r tH the response of the flight-path angle to control move- 
ongitudina otions may give an understanding of the , ‘ . . 
a am _ dia | 6 : ments between swept-wing airplanes having no hori- 
results just shown. Some of these parameters for air- : , am ¢ : ; 
; 4 | relati | f aj zontal tail and using trailing-edge flaps for longitudinal 
janes B and C are compared relatively to those of air- : 
, : x of = ; ; mplified Ivtical control and conventional airplanes that are known to 
jane . Table 4. From the simplified analytica cr er : 
— . a , _ ; , ith ; have satisfactory response characteristics. Particular 
consideration, it is indicated that these parameters are a : . . , 
me ” — ; : ite emphasis is placed on the final few seconds of the land- 
indicative, comparatively, of the motion to be expected. 
: ' ae ing approach. 
at the time of control movement. The parameter ; ; 
Ta 2 tn: : This study shows that the only difference in the 
V2Cm; ‘2uk,? indicates the effectiveness of the elevator - ; ; 
cl ‘ , le of tl ‘tcl \j short-time response of the flight path between the two 
on the rate of changing the angle of the pitch. <Aiur- “i a ; a‘ ‘ 
- a Fecti a : 4 7 types of airplanes is in a time lag in the flight-path 
anes B< * were less effective than airplane A. 1e : : . : oa 
_— - — ‘i “9 ; , bili F , ‘tched angle response of the swept-wing airplanes. This time 
yaramete ../ 2ué indicates the ability, when pitched, ‘ m . . 
parameter 'Cz,,/<ué indicate : P ; lag was found to result from two factors: first, a rela 
o change lift refore the al acceleration anc ; ; 
to change lift ind therefore the normal accele ‘ tively small value of the parameter 
flight path. This parameter was smallest for airplane 
B and had about the same value for airplanes A and C. 
The product of these parameters and the available 
elevator deflection (called the total elevator effective- 


c 


(VCz,,/2uc) (V*Cmg, 2uk,*) Ad, 


which is a measure of the effectiveness of the control in 
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changing the flight-path angle, and, second, a relatively The importance of the differences found in the 
large value of the parameter sponse characteristics between the two types of air 


Tay “\AS lanes studied can only be evaluated by flight exper 
(J C15, /2ue) Ad, I : 4 : gs xperi 


ence. Other factors such as range of vision, contro) 


which is a measure of the amount of undesirable change _ feel, the pilot’s reaction to the relatively large nose uy 
in lift accompanying the control deflection required to attitudes of the low aspect ratio swept-wing airplanes 
change the angle of attack. It had been suspected that and psychological influences associated with new type 
the relatively large change in drag with angle of attack airplanes may be of equal or greater importance. 

of low aspect ratio swept-wing airplanes having no 
horizontal tails would be an important factor in the 
flight-path angle response but this was found not to 
be the case for the short time intervals under consider- 
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Supersonic Flow over Conical Bodies 
Without Axial Svmmetry 


ANTONIO FERRI,+ NATHAN NESS,? ann THADDEUS T. KAPLITA** 
Polytechnic Institute of Brooklyn 


SUMMARY 


Conical flow around cones without axial symmetry is analyzed 
by the superposition of linear conical solutions to a nonlinear 
ixially symmetric conical flow. The aerodynamic properties of 
ones with elliptical, drop-shape, and triangular cross sections are 
determined from these linearized solutions. It is shown that 
some gain on the shock drag can be obtained by using bodies 


having cross sections different from circular 


INTRODUCTION 


_ PRACTICAL DESIGN of any fuselage for airplanes 
or missiles is obtained normally from an engineer 
ing compromise of aerodynamic, structural, and in- 
stallation requirements that are peculiar to each par- 
ticular design. One of the parameters obtained from 
such a compromise is the shape of the cross section of 
Since the aerodynamic properties—e.g. 
are affected by the cross-section contour, 


the fuselage. 
lift and drag 
methods that permit the analysis of a supersonic flow as 
afunction of this parameter are of practical interest. 

The pressure distribution around slender bodies with 
out axial symmetry can be obtained by the linearized 
characteristic method discussed in reference 1. How- 
ever, this analysis can be performed only when the flow 
field near the apex of the slender body is known. This 
flow field, which is of the conical type, is analyzed by 
the method outlined in reference 1. The present paper 
is an extension of this method, with special emphasis 
placed on its practical application to conical bodies 
without axial symmetry. 

This method consists of the superposition of linear 
ized flow solutions upon the solution of the nonlinear 
flow around a circular cone. Because the same linearized 
flow solutions can be used regardless of the difference 
in cross-section shape between the circular cone and any 
of the other bodies considered, this procedure is espe- 
cially effective for the analysis of the effect of cross sec 
tion on aerodynamic properties. The equations that 
define the velocity components of the linearized flow 
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fields are solved by a step-by-step numerical procedure 
(i.e., the hodograph method). 

The present procedure permits the determination of 
the shock wave created by nonaxially symmetric bodies 
and, hence, gives solutions also at high Mach Numbers, 
where methods based on linearized théory (superposi- 
solutions on a uniform flow field) 
the latter case the Mach 


tion of linearized 
cannot be used because in 
cone at high Mach Numbers falls inside the body. 

In this paper, the flow fields around conical bodies of 
elliptical, drop-shape, and triangular cross section for a 


wide range of Mach Numbers are considered. 


List OF SYMBOLS 


r,v, 0 = spherical coordinates 
Un, @ = velocity components in spherical coordinates 

a = speed of sound 

as = entropy 

R = gas constant 

y = ratio of specific heats 

V; = undisturbed velocity, referred to limiting ve 
locity 

p = free-stream density 

3 = defined by Eq. (8) 

y y = coefficients of the Fourier series defining non 
axially symmetric shock surface 

Yur, Um, = coefficients of the Fourier series defining non 
axially symmetric conical body surface 

H = radius of hodograph diagram 

VW, = free-stream Mach Number 

( pressure coefficient 

iz = specific heat at constant volume 

V local resultant velocity, referred to limiting ve 
locity 

a/b = ratio of semidiameters of ellipse 

h = distance along body axis from apex 

rab = cross-sectional area of ellipse 

ar = angle of attack 


= lift, drag, and horizontal force coefficients re 


ferred to cross-sectional area normal to body 


axis 


CN = normal force coeflicient referred to cross-sectional 
area in the meridian plane perpendicular to 
normal force component .\ 

\ normal force component 

H = horizontal force component 

L/D lift-drag ratio 

dc; da lift curve slope 

x tan”! (a/h) 


Subscripts 
a properties of the basic axially symmetric conical 
flow 


n,m = properties of the superposed linearized flow fields 


908 
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—AXIS OF THE POLAR 
COORDINATES 


Fic. | Phe polar coordinate system and velocity components 
s = properties at the shock wave 
b = properties at the body surface 


DESCRIPTION OF THE METHOD 


Consider a conical flow field. Assume a spherical 


coordinate system 7, Y, J whose center is at the apex of 


the conical flow. Denote by v, the velocity component 


along the radius 7; v,, the component normal to the 
= constant; and w, the 


radius in the meridian plane J 


component normal to the meridian plane (Fig. 1). By 


combining the laws of conservation of mass, momen 
tum, and energy, the following equations can be ob 
tained for a nonviscous perfect gas without heat con 


duction :* 


rT we? WU, Q, Ow 
F (1 — ) - ( ; +. ) eQ@ (2) 
sin Y Ov a a” \sin Wy Ov oy 


a’ os ; Ow Ov 
=v, sin y -v 
yR dv oy ov 
Or, 4 : y + ! ) 
Un wv’, Sin wz, cosy (2) 
ov 
a* OS Ov; _ Ow w Ov, 
yR ov Ow dy | sin Ov 
v7, — weecoty (3 


where a is the speed of sound, S is the entropy, R is the 
gas constant, and y is the ratio of specific heats (as- 
sumed constant). 

A conical body traveling at supersonic speed creates a 
shock front that is conical with apex at the apex of the 
conical body. This shock front separates the upstream 
undisturbed air from the downstream nonuniform flow 
field and may be expressed in spherical coordinates by 


Y%=t+¥,2+ 2 ,cosnd + D> yp, sin md (4) 


" 0 m 1 
In the classical linearized theory the shock surface is 
coincident with the Mach cone and has a circular cross 


section; therefore, in the classical theory, the terms 


Yns and Wms 


is equal to the Mach angle. 


are considered equal to zero, while the 


angle y,, In this analysis 
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the terms y,, and y,,, are considered small but different 
from zero, and only first-order terms of y,, and y,,. are 
retained in the differential equations and in the bound 
ary condition. The term y,, is assumed different from 
the free-stream Mach angle. 

From physical considerations across the shock, the 
velocity components immediately downstream of the 
shock front (1.e., at Y = y,) are obtained. ! 


(v,)y I", cos x 5 
(Vy )y 
¥ (? 1’)? cos? ¥ 1)? sin? y ant) 
y+ 1 I, sin y 
I; sin yy, sin? 8B (6 
(vw) y I, sin y, sin B cos B +4 
Y I 7 , > Swe —” 
; ; 1? cos? y I")? sin? y, sin? B) X 
Y 
tan > 
Sear (7 
I", sin y, 


where all the velocity components are nondimensional 
ized with respect to the limiting velocity, I; is the free- 
stream velocity (also referred to the limiting velocity), 
and 3 is defined by 


tan 8 = Oy,/sin p,dv (8 


As only first-order terms of y,, and y,,, are being con 
sidered, any velocity component behind the shock 
created by the nonaxially symmetric body (at y = 
y,) may be expressed as a function of the same velocity 
component at the shock created by the axially sym 
metric body (at Y = y,,) and of the derivative of the 
velocity component with respect to Y multiplied by the 
difference between y, and y,,. In this way, velocity 
components consistent with the shape of the shock [de- 
fined by Eq. (4)] and with the shock relations [defined 
by Eqs. (5) to (7)] are obtained in the form 


v, = Urq + - k WasUrn cos nd + >. Y, 5 Ury 
n 0 m l 


sin md (9 
Uy = Ung Do Wn nn CoSNI + YS Yny tn, Sin md (10 
n 0 m l 


— 


w= > ny,.W, sin nd + » MY Wy, COS mI (11 
1 ) 


m ( 

The quantities v,,, Ungar Ura» Urmy Una» Unm» Wn,» and w 
are functions only of Y and are independent of 3; v,, and 
v,, are the radial and normal velocity components of a 
conical axially symmetric flow field and may be obtained 
At the angle y = y,,, the velocity 


»Unny Unms Wn, and w,, are independent 


from reference 3. 
components 7,,,, Urm 
of m and m and are given by the following expres 
sions: 

(12 


(z,, ) = (v,,,) = — [Il;sny+2,, 


Yd. a 








where 
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77? = (7 ) = 
Un) Ve m/ va 
7 l < 1 — V," cos? y OVna 
= cos Y 2] = = 57° ¥ 
7 1 sin? y OW Jia 
(13) 
(Wi)yg, = (Wmv. = Vi + (ng/sin ¥)y, (14) 
where 
Craly ry COS Ya. ( 15) 
y-—1 ( — JV," cos? *) : 
= = — (16) 
ong Y + I V, sin p y 


a 


(~") | Ura 2 Una cot y 17 
= —|o,,+ : (17) 
Oy /y, 1 — (_/Qa) va, 

When the velocity components, as defined by Eqs. (9) 
to (11), are substituted into Eq. (1)‘and only first-order 
terms of y,, and y,,, are retainea, tne toliowing three 


independent equations result: 


ay Wne = 
Pe (2 — :) + v,, cot Y + (1 _ ‘) = 0 (18) 


Qa oy \ a.” 
OVnn z n-w 
n a n 
(: oor l — = — Agr, — Barn 
oy Qa sin W 
(19) 
OUnm . a 
Ve + 1 — = 
Oy a 
MW m 
—- Adm — Bing —-—,_, 
sin y 
where 
Ye 
2 y-—I1 
a 
1, = coty + ( + Vue) 
Cy a," 
Vig + Ung cot ¥ 
l1— io? ‘a?) 
Y — 1] U%q + Ung cot 
B, T v a “na 


a," Ll — (Un_°/Aa? 


Eqs. (18), (19), and (20) are independent of J. Eq. 
(18) is the ordinary nonlinear differential equation for 
conical axially symmetric flow whose numerical solutions 
are given in tabulated form in reference 3. Eq. (19) 
is an ordinary linear differential equation independent 
of Y,,, and Eq. (20) is an ordinary linear differential 
equation independent of y,,,.. Therefore, the flow field 
considered is represented by a superposition of linearized 
, Wn, m [the pro- 


a 


disturbances proportional to 7,, 
portionality coefficients are the shock-front coefficients 
(Yn, m)s}] upon the nonlinear axially symmetric conical 
flow defined by Eq. (18). 

By combining Eqs. (2) and (3) and the momentum 


” 
m? “Nn, wm 


equation in the radial direction, the relation 


v, sin ¥(0.S/Oyw) —w(dS/dd) (21) 


OVER 


CONICAL BODIES dHS 
can be obtained.” From Eq. (2), for the approximation 
accepted, 0S/O# is of the order y,,; therefore, when < 

is at ieast of one order larger than y,,, Eq. 21 indicates 
that 0.S/Oy is of the order of (y,,)? and can be neglected. 
(The identical reasoning applies when the subscript 7 
is considered; henceforth, in this section, onlv the sub 


script m will be used.) Then Eq. 3 gives 


Ung = OW _/OW, Uny Ov,,,/O¥ (22 
while Eq. 2 results in 
Ge* - ” “s y Ow, 4 or 4 
yR One sin oy Urqlrn 7 
VnaUnn + Urq, SIN W + Ung W, COS P (23) 


Eqs. (22) and (23) are independent of 3; S’ is a con 
stant obtained from the shock equation at the shock 
front ~,,—1.e., 


S’ = (dS/dW)y, 
Eq. (23), with the use of Eq. (22), can be written as 


a" sy \ ae 
- S = 2,,(0;, + Wn, Sin yy) 4 
yR 


When Eq. (24) is multiplied by the integrating factor 


1/(> 1) 
(a Vna' 


—Unq sin p) 
the following expression is obtained: 


Un + WwW, Sin = — (=) a” oF. 
. - yR/ * 


3)/ (4 


y <7 

r 2 Qa - — 

(—Una sin y) 2 . a\? dy -") 
J ta. Ung Ung SN ¥ 


a 


7) 


where the constant of integration is zero because at ¥ 
¥., [from Eqs. (12) and (14)] 


Wy, = —?,,/sin 


26 

In Eq. (19) all the coefficients of the perturbation 
velocities and in Eq. (25) the right-hand side of the 
equation are known for any value of ¥; therefore, the 
quantity 


ma /- 
Urn T 


(Ov, OW) 


which is the radius //, of the hodograph diagram (for 
example, see reference +), can be evaluated at any y¥. 
‘ony Ung» and w,, for any value of m may 


The quantities ¢ 
and the 


then be determined between the shock front ¥ 
surface of the circular body y,,, (which created the shock 
front y,,) with the use of Eqs. (19) and (25) and the 
following step-by-step process: 
-Ay) + (Hn — Urndy X 

97) 


sin (— Ay 27 


COS | 


(Uny)y 


UT, — Ven) x 


+ 


cos (— Ay) + (1), 28 


) ma 


Urnly Ay Cnn 


)y sin (— Ay) 
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BOUNDARY CONDITIONS 


The conical body that produces the shock front defined by Eq. (4) can be expressed in the form 


YY = Ya, + pm VY», COS nd? + 7 Vimy Sin med (99 
n=0 


The boundary condition requires that the body surface be a stream surface for the conical flow field. 


bodies, this condition can be expressed as 


m=1 


For conical 


Un l oY, l ° s . _ ss 
== =— [— Xny,, sin nF +Z my, cos md] 20) 
w/y, sin p,ov sin Wp 
But at the body surface Wp, (v,)y, may be expressed in terms of a Taylor expansion about the basic circular cone 
Vas —t.e., 
Ov, 1/O07v (Wx — Way)? 
la na ) ab . i . 
Cy, = (6. oF ( (Yo — Wary) + + TWnz (nny. COS AI + Thing (Unmy. X 
a/¥q, , ) 2), 9 1p ' Yap 
“s oY Yap = oY Yay ic 
; wa {ome = OV'nm ; 
sin md + | Lyn, > cos 20 + Lyn, > sinmd | (Wy — va) t+... (3 
¥ Vay v “ah 
(wW)y, = UWns (W)¥,, n sin 20 + Lng (Wm) vq, m cos may + 
up « by 
om Ow, ; ™ OW in 
LW n sin nd + Lym, m cos md | (Wp — wa) +... (32 
oy Yay OY Ya, 


Therefore, when Eqs. (31) and (32) are substituted into the boundary condition [Eq. (30)], relations defining 


the coefficients y,, and y,,, can be obtained. 


These relations will differ for different assumptions of the relative 


order of magnitude of the body surface coefficients y,,,, ¥,,, With respect to the shock surface coefficients y,,,, Yng. 


For high Mach Numbers or for large cone angles, the shape of the shock follows closely the shape of the body, and 
the coefficients Y,,, Ying, Yass Wms are of the same order; therefore, terms of the order of Wy Wry, Vanes Vno?r mg? OF 
higher can be neglected in the boundary conditions, as well as in the differential equations of motion, and linearized 


boundary conditions can be used. 


oY 


Ov 
\ Na . . bd 
(Cray, i ( (Ly, cos nv + Ly,,,, sin ms) 
, Ya, 


If, at a given Mach Number, y,, is the shock angle 
produced by the circular cone y,,, then (v,,,)y.. = 0, and 
: Vap 


(OVng Wy, as defined by Eq. (18) is 


(OVng oy My = 


ap 


= ae 
o(U rq Ya, 


Therefore, the boundary conditions as obtained from 
Eq. (33) are 


y _ 2y (v ” ) 
nes ~— —¥nb \"ra/ “nn va 
b 
= 9 ” 9 
Wing aa 2Ymp (2 ra/ Unm y \ 


These equations relate the shock coefficients y,,, 


ap 


Wms to the radial velocity of the conical axially sym- 
metric flow v,, (which can be obtained in tabulated form, 
for example, from reference 3), to the normal pertur- 
bation velocities 2,,,,, U,,, [which are obtained from a solu- 
tion of Eqs. (19) and (20)], and to the body coefficients 
Ynor Ymp Which can be calculated once the body shape is 
prescribed. 

from Eq. (4). 


Therefore, the shock shape can be obtained 


For these linearized boundary conditions, Eq. (80) simplifies to (v,,) 


= Qor 


vr 


ZY me (Une) 


mie 


= — Ly,,, (Uny vay COS NO — sin mo (33 
; 


at 


For sinall cone angles and small Mach Numbers, the 
values of y,,, and y,,, are larger than the values of y,,, and 
Y»., and, therefore, higher order terms of y,, and Wn 
can be retained in Eq. (30). If higher order terms of 
Yao» Wmpy are retained in the boundary condition, a sys- 
tem of linear equations is obtained where each equation 
contains the unknowns y,,,, W,,,; therefore, in this case, 
the determination of the shock coefficients is much more 
involved. In this paper all the results are based on the 
linearized boundary conditions as given by Eq. (34) 

If the flow field has a plane of symmetry and the axis 
of the conical flow is in this plane of symmetry, the 
shock equation (as well as the equations for the velocity 
components) may be expressed in terms of only one sum- 
mation. If the plane of symmetry is the plane 3 = 0, 
If the flow field 
has two planes of symmetry, then only the even coefti- 


vd = , the coefficients y,,, are zero. 
cients Yo,, Yrs... - » Wns, Where 7 is an even integer, ap- 
pear in the summation. Furthermore, if the value of 
Was corresponds to the value of y,, for axially symmetric 
flow, Wo, is zero. 
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SUPERSONIC FLOW 


The particular case of the flow around circular cones 
at small angles of attack, where the plane of sym- 
metry is the plane 3 = O, 3 = 7, is given by » = 1, 
m = 0. This analysis was considered first in reference 
i. The present method is a generalization of this anal- 
ysis to cones without axial symmetry. The analysis 
of the nonaxially symmetric cones is obtained by means 
of numerically determined linearized solutions pre- 
sented in tabulated form.’ 


APPLICATION OF THE METHOD 


In order to obtain aerodynamic results from the pre- 
ceding analysis, several linearized flow fields have been 
determined for different circular cone semiapex angles 
y,, and free-stream Mach Numbers 1/;. The linearized 
solutions have been determined for the basic flow solu- 
tions shown in Table 1. For each basic flow solution, 
ten linearized solutions have been determined corre- 
sponding to values of ” or m equal to 1 to 10.* These 
values of the linearized flow field (which will be pub- 
lished shortly in tabulated form’) were used to obtain 
the flow field around bodies of different cross section 
ind to determine the effect of cross-section contour on 
the shock shape, on the pressure distribution, and on 
other aerodynamic properties. 

The pressure distribution is obtained from the ex- 


pression 


‘ 79 y¥/(y¥-1 
2 P=) gn ASlely— 0) | 
°~ yMYL\ — V2 


where I is given by 


[7 = [( + TWas (Urn)y 


cos nd + 


ral 
a/Yap ab 


ps sin md]? + [Xnyp,,.(w,)y,, sin nd + 


(Ure) y 
' Yab 


LMY », (Wm) cos md}? 


ab 

The values w, (and w,,) are defined by Eq. (26), 
which is valid along a stream line. The increase in 
entropy AS between the shock front and the surface of 
the body is considered constant along the body surface 
between meridian planes where the w velocity com- 
ponent is equal to zero, and its value is that obtained 


* The linearized solution for m = 1, which gives the effect of a 
small angle of attack on a circular cone must correspond to the 
linearized flow field given in the M.I.T. tables’ when both flow 
fields are referred to the same axis. The velocity components of 
the M.I.T. tables are proportional to the body parameter yy, 
while these components are referred to the shock parameter 
vi, in the present analysis, therefore the latter values must be 
multiplied by 6/e of reference 5 in order to make a comparison. 

A systematic error, however, has been found in the values 
of reference 5 due probably to an error in sign in the evaluation 
of the integral of Eq. (24) of that report [which corresponds to 
Eq. (25) of this report]. This error affects all the velocity 
components and is most pronounced for the w component, es- 
pecially at high Mach Numbers and large cone angles. This error 
probably also affects all the values of reference 8 which are deter- 
mined from the linear solution. 
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SECTION A-A 


Variation of shock wave with free-stream Mach Number 
for an elliptical cone at zero angle of attack 


Fic. 2. 
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SECTION A-A 
Variation of shock wave with free-stream Mach Number 
for an elliptical cone at zero angle of attack 


Fic. 3 


in the meridian plane # = constant where the stream 


lines diverge. 
Elliptical Bodies 


At zero angle of attack only the even values of n 


(xn = 2,4,..., 10) are necessary, while, because of the 
linearized boundary conditions, only » = | is required 
to express the effect of angle of attack and only m = | 


is required to express the effect of angle of yaw. 

A preliminary investigation of the accuracy of the 
linearized procedure for purposes of comparison with 
the experimental results that were available for 1/,; = 
1.81 and a/b = '/; indicated very good agreement.! 
In the present paper, the investigation of elliptical 
Mach 
Numbers for ellipses with semidiameter ratio a} 


'’; and for different values of h +v/ rab 


bodies is extended to cover a wide range of 


2/3, */o, and 


where /i is the cone height and zab is the cross-sectional 


TABLE 1 
Yo, = 5 7.5 105 12.5 20 25 
1.467 
1.455 1.871 
1.596 1.881 2.4438 1.825 
M, = 2 685 2.779 2.891 3.369 3.926 
5.800 5.1038 5.422 5.860 5.546 5.123 
9.100 1] ».593 8.096 


8.492 8.059 091 ¢ 

















568 JOURNAL OF THE AERONAUTICAL SCIENCES—AUGUST, 1953 
TABLE 2 _ 
Coefficients Defining Shock-Wave Cross Section [See Eq. (4) | 
Ellipse Ellipse 
a/b = '/. a/b = '/. Circle Drop-Shape Triangle 
Shown in Fig(s). 2 2,18 18 18 18 
VW, 1.60 9.10 9.10 9.10 9.10 
Was 0.693 0.219 0.219 0.219 0.219 
Wie 0 0 0 —0.0134 —().0100 
2s 0.00162 0.0399 0 0.0196 — 0.00560 
33 0 0 0 0.00600 0.0198 
V4, 0.00001 0.00487 0 0.00166 —0.000521 
V5. 0 0 0 0.000573 — 0.000775 
Vos 0 0.00065 0 0.000322 0.00176 
v7; 0 OQ i) 0.000166 —0.00174 
vs, 0 0.00010 0 0.000059 —0.000119 
V9, 0 0 0 0.000021 0.000155 
Yr0 0 0.00002 0 0.000012 0. 0000002 
TABLE 3 
Coefficients Defining Conical-Body Cross Section |See Eq. (29) | 
— Ellipse Drop 
a/b ='/2 a/b=1/3 a/b = '/2 a/b =? Circle Shape Triangle 
Shown in Figure 3 2 ‘ 13 14 
Semiapex Angle of Basic Circular Cone 25° 10° 10° 10° 10° 10° 10 
Yap 0.456 0.175 0.175 0.175 0.175 0.175 0.175 
vin 0 0 0 0 0 —().0150 —(0.0112 
v2» 0.128 0.0879 0.0577 0.0344 0 0.0284 —(). 00809 
Vp 0 0 0 0 0 0.0122 0.0401 
Wap 0.0281 0.0325 0.0142 0.00506 0 0.00485 —0.00158 
Wop 0 0 0 0 0 0.00239 — (0.00324 
Wen 0.00687 0.0133 0.00387 0.000826 0 0.00190 0.0104 
Vin 0 0 0 0 0 0.00138 —0.00014 
vs, 0.00177 0.00565 0.00110 0.000141 0 0. 00068 —(). 00136 
vo» 0 0 0 0 0 0.00032 0.00241 
Yio, 0.000470 0.00246 0.000323 0.000025 O 0.00025 — 0. O00004 
‘ area. Some of the more interesting results are pre- 
a | me er ar euamaehoats y sented in the figures. 
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Fic. 4+. Variation of pressure coefficients around an elliptical 
cone with free-stream Mach Numbers at zero angle of attack 
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Fic. 5. Variation of pressure coefficients around an elliptical 
cone with free-stream Mach Numbers at zero angle of attack. 





In Figs. 2 and 3 are shown the shock shapes for two 


cones with the same elliptic cross section (a/b = |)» 
both at zero angle of attack, but with different ratios of 
h’~V/ rab. The coefficients defining the cross sections 
of the various bodies considered and the coefficients of 
the cross sections of some of the resulting shock shapes 
It will be noted that, 
for small Mach Numbers and small cone angle, the 


shock is close to a Mach cone, but for high Mach Num- 


are tabulated in Tables 2? and 3. 


bers or large cone angle the shock shape is far from 
circular. To realize the limitations of the linearized 
theory, consider, for example, the body in Fig. 2. For 
Mach Numbers of the order of 3, the Mach cone would 
intersect the ellipse; therefore, the linear theory or any 
other higher order methods that start from linear solu 
tions cannot give results in this range of Mach Num- 
bers. 

In Figs. 4 and 5 the pressure distribution for different 
Mach Numbers around the two cones of Figs. 2 and 3 
are presented. 

In Fig. 6, the effect of the variation of the parameter 
ab of the ellipse on the pressure coefficients at a con- 
stant Mach Number and with approximately constant 
h/V/ rab ratios is presented. (These conical flows were 
obtained from the same basic axially symmetric conical 
flow corresponding to a circular cone of 10° semiapex 


angle at J; = 2.78; the cross-sectional areas of the 
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SUPERSONIC FLOW 
elliptical cones are not exactly the same, but the vari- 
ations are small. ) 

In order to determine the effect of the variation of the 
 bratio on the shock drag, the drag coefficient of each 
dliptical cone is divided by the drag coefficient of a 
circular cone of exactly the same cross-sectional area. 
The results are given in Fig. 7. As expected, the shock 
drag decreases as the ratio a/b approaches zero. 

The effect of angle of attack (or angle of yaw), within 
the accuracy of the linearized boundary condition, is 
considered by superposing upon the linearized flow solu- 
tion of the nonaxially symmetric cone at zero angle of 
attack (or vaw), the flow solution for a circular cone 
(from which the nonaxially symmetric body is pertur- 
bated) at the angle of attack (or yaw) at which the 
actual body is being analyzed. 

In Figs. S and 9, the pressure distribution at an angle 
of attack or vaw of 6° and J, = 2.78 for the cones con- 
sidered in Fig. 6 is presented. In Fig. 10, the lift, drag, 
horizontal, and normal force coefficients are shown; the 
lift, drag, and horizontal force coefficients are referred 
to the area of the cross section which is normal to the 
body axis; the normal force coefficient is referred to the 
cross-sectional area in the meridian plane which is per- 
pendicular to the normal component N. Also given 
in Fig. 10 is the normal force coefficient cy, as obtained 
from linearized theory, of a triangular flat plate of 
plan-form area equal to the respective ellipse. It can 
be seen that L/D of the body increases rapidly when 
a/b increases. 

In Fig. 11 the drag coefficient as a function of Mach 
1/,. but for dif 

The value of 


Number of elliptical cones with a/b = 
ferent values of h/+/ rab, is presented. 
¥a of the basic circular cone corresponding to each 
value of h/\/ rab is also given. In Fig. 12 the effect 
of Mach Number and /,/+/ rad on the ratio between the 
drag of the elliptical body to a circular cone of equal 
cross-sectional area and height h is presented. The 
ratio is less than unity for all Mach Numbers and / 
V rab ratios and changes slightly when 1/,; or h/V rab 
changes, indicating that some gain on the shock drag 
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Fic. 6. Variation of pressure coefficients around different 


elliptical cones at the same free-stream Mach Number at zero 
angle of attack 
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Fic. 7. Variation of drag coefficient of elliptical cone Cp, as 
function of a/b ratio of ellipse. (Drag coefficient is referred to 
circular cone of same height and cross-sectional area as ellipse. ) 
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Fic. 8. Pressure coefficient around elliptical cones of different 
shape at 6° angle of attack 
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Fic. 10. Lift, drag, normal, and tangential coefficients for ellip 


tical cones with different a/b ratios at 1\/; 2.78 and a = 6°. 
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Fic. 11. Drag coefficient as a function of Mach Number 


for elliptical cones of equal cross-sectional areas for different 
values of h/V mab 
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Fic. 12. Variation of the ratio of the drag of an elliptical 


cone to a circular cone of same height and cross-sectional area 
as a function of Mach Number for different values of 1/V ab. 
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SECTION A-A 
ea 
Fic. 18. Variation of shock wave with free-stream Mach 


Number at zero angle of attack for a cone with drop-shape cross 
section (parabola tangent to portion of a circle). 


can be obtained by using bodies having cross sections 
different from circular. 


Drop-Shape and Triangular Shape Bodies 


In Figs. 13 and 14, the shock contour at different 
Mach Numbers for two other conical bodies is presented. 


The cross section of the drop-shape cone of Fig. 13 is 
obtained by fairing a portion of a circle to a parabola. 


this cone is about the same as the corresponding ratio 
of the elliptical body shown in Fig. 2. The value oj 
h/V/area for this cone is 3.16; the corresponding ratio 
for the elliptical cone of Fig. 2 is 3.10, while this ratio 
The 
ratio between the maximum width to maximum height 
1.34. 

The cross section of the cone of Fig. 14 is a triangle 


for a circular cone of 10° semiapex angle is 3.21, 
at any cross section is | 
whose corners are faired by circular arcs. The ratio of 
the height to the cross-sectional area of this body js 
3.14, practically the same as the other bodies mentioned 
above. Its ratio of maximum width to height is 1.26 
In Figs. 15 and 16 the pressure distribution around the 
two bodies is presented. 

It has been found that the drop-shape body at .\/, = 
2.78 has zero lift at a = —1°, its minimum drag at 
a = —2”, and a lift curve slope (dc,/da),~» 1.72 
rad.; similarly, it has been found that the triangular 
shaped body at the same Mach Number has both zero 
lift and minimum drag at about a = 0° and a lift curve 
slope (dc, /da), =) = 2.41 rad. 

In Fig. 17, the drag coefficients of the circular, ellip- 
tical, drop-shape, and triangular bodies as a function of 
Mach Number are presented. Curves for the drop- 
shape body at a = 0 and a = —2° are shown, while the 
other bodies are at zero angle of attack. 

The value of the drag coefficient has been divided by 
the drag coefficient of a circular cone of the same area 
and height in order to eliminate the effects of small 
It will be noted 
that the drag ratio remains practically constant for a 


differences due to the ratio h/+V/ area. 


wide range of Mach Numbers and is less than unity 
for all bodies other than circular. Therefore, as in the 
case of elliptical bodies, some gain in shock drag may 
also be obtained from the drop-shape and _ triangular 
bodies. Furthermore, the drag ratio for the drop-shape 
body is less at a = —2° than at zero angle of attack for 
In Fig. 18, 
the shock shapes for the various bodies are shown for 
AT, = 9.10, a = 0°. 

It is important to mention that the linearized ap- 


the range of Mach Numbers considered. 


proximation in this analysis is not the same for all bodies 
considered because the intensity of the linearized flow 
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Fic. 14 of shock with free-stream Mach 


Variation 
Number at zero angle of attack for a cone with faired triangular 


Wave 


The ratio of the height / to the cross-sectional area of 


cross section (corners faired with circular ares) 
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jelds varies with different boundary conditions. These 
simplified boundary conditions introduce different ap- 
proximations for the various body shapes. However, 
the linearized flow fields are all small, and it is expected 
that the trend of the variations given in this report is 


correct. 
CONCLUSIONS 


The method of superposing linear solutions to a non- 
linear nonuniform flow field has been applied to the 
analysis of conical flows without axial symmetry. 

For elliptical cones, the results indicate that, at a 
constant free-stream Mach Number, the shock drag 
ratio (i.e., the ratio of the drag coefficient of the ellipti- 
cal cone to the drag coefficient of a circular cone of the 
same cross-sectional area and height) decreases as the 
ratio of the minor to major axis of the ellipse decreases. 

Furthermore, for a given elliptical body at zero angle 
of attack, the drag ratio remains practically constant 
and less than unity for the range of Mach Numbers 
considered. 

The effect of angle of attack is also analyzed, and the 
lift coefficient is determined. The lift coefficient is not 
far from the lift coefficient of a triangular flat plate of 
the same plan form. 

Cones having drop-shape and triangular cross sec- 
tions are also analyzed. The drag ratio of these bodies 
also remains practically constant and less than unity 
fora wide range of Mach Numbers. Therefore, the re- 
sults indicate that some gain on the shock drag can be 
obtained by using bodies having cross sections different 


from circular. 
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Fig. 15. Variation of pressure coefficients around a drop-shape 
cone with free-stream Mach Number at zero angle of attack 
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shaped cone with free-stream Mach Number at zero angle of 
attack 
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Fic. 17. Ratio of drag coefficients of conical bodies to the 


circular body of same height and cross-sectional area as a function 
of Mach Number 
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On the Nonexistence of a Finite Critical Time 
for Linear Viscoelastic Columns 


Joseph Kempner and Frederick V. Pohle 


Assistant Professors, Department of Aeronautical Engineering and 
Applied Mechanics, Polytechnic Institute of Brooklyn, Brooklyn, 
i 


May 7, 1953 
SYMBOLS 
E = Young’s modulus 
I = moment of inertia of the cross section with respect to a 
principal axis 
L = length of column 
\J’ = bending moment 
P = axial compressive end load 
Pe = xr*EI/L? 
Co = eccentricity 
n = (P/P)>1 
t = time 
“u = x/L 
x = axial coordinate measured from one end of column 
y = deflection due to loading 
é = t/r 
rT = relaxation time 


INTRODUCTION 

I A PAPER PRESENTED IN 1946, Freudenthal! concluded, for 

columns fabricated from a linear viscoelastic material, that a 
finite critical time existed. Since that time, the same conclusion 
has been repeated by Freudenthal in his book? and apparently has 
been accepted by Rosenthal and Baer.’ Libove‘ stated that no 
finite critical time existed for such columns and referred to an 
unpublished letter by the present authors. Hilton® referred both 
to the letter and to Freudenthal’s work. The purpose of this 
note is to show briefly that Freudenthal’s conclusion is incorrect 
and at the same time to make available the contents of the com- 


munication referred to above 


DISCUSSION 

Freudenthal analyzed a simply supported column fabricated 
from a Maxwell material and with a constant eccentricity ¢ 

The basic differential equation that governs the deflection is 
EIO(d?y/dx*?)/ot = —(OM/dt) — (M/r) (1) 
The method used by Freudenthal to solve Eq. (1) was that of 
successive approximations. The first approximation to ./ was 
assumed to be Peo. The differential equation for the first ap- 


proximation y, to the deflection y is 


EI0(d?y,/dx?)/ot = —Pco/r (2) 


of 


9 


The solution of Eq. (2) for a simply supported column under 
going purely elastic deformation at zero time is 


M1 = Co(w*/n) (1/2!) (1 w*)(1+¢ 
where u = x/Landé = t/r. 
A correction to the assumed moment is thus A./ = /’y Sin 


Eq. (1) is linear, the increment in moment when substituted int 
the right-hand side of Eq. (1) yields the following different; 
equation for ye, the correction to y;: 


O(d?y2./dx?)Ot = —Co( w?/n)? (1/2!) (1 u*)(E) (2 + £ 
The procedure used to obtain ye can be applied repeatedly t 


obtain the deflection in the form y = y,; + yo + y 
The result is 


m l 
where 
g(u) = u u? 
g(u) = u 2u* + ut | 
g(u) = 3u 5u3 + 3x! ue f 
gs(u) = 174 28u? + 14u tu? + uw 
and 
f(t) = 1 +4 
f(t) = 1 + 2 + (1/2)E | 
ff) = 1 + BE + (3/2)E° (1/6 7 
f(z) = 1 + 4¢ + 382 


The recursion formula for f,,,(£) is 


re 
fn(i) = { fm We) dé + fm—ali a & Tt: : l S 
J0 


The relation (8) was not stated by Freudenthal, but it will be 
essential for the determination of the general term of Eq. (5 





The maximum bending moment occurs at x = L 2 and is 
Mmer. = Peo{1 + 1.2337(1/n)f,(é) 4 268301 / 2 fol & 
1.2727(1/n)f3(¢) + 1.2728(1/n)*f,( 5) + Hy 
A Fourier analysis of the polynomials in « = x / which occur it 
Eq. (6) shows that the numerical coefficients of (1 1 )™f,,(&) it 


Eq. (9) are given by 


(4/mr) [1 — (1/3)2"*1 + (1/5)? { 
Thus the numerical coefficients have 4/7 = 1.27524 as an upper 


bound. The corresponding coefficients stated by Freudenthal 
were 1, 1.22, 1.26, 1.26, 1.26, ...; these coeflicients were not 
computed accurately, nor was any proof given that the sequence 
of coefficients possessed an upper bound 

Clearly, the convergence or divergence of the infinite series im 
Eq. (9) depends entirely upon the properties of 
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YS (1/n)"fn(8) (11) 


m=1 


g(é) = 
The function ¢(£) can be evaluated in closed form in the following 
way: If both sides of Eq. (8) are multiplied by (1/n)" and sum- 
mations are carried out over m from 1 to ©, the left side of the 
result will be y(), but the right side will contain an integral of 
{t), Differentiation of the result with respect to £ yields the 


differential equation 


dg/d& — [1/(n lj], = 1/(n — 1) (12) 
The solution of Eq. (12), with g(0) = 1/(m — 1)is 
e(£) = [n/(n — 1)le® ee ae (13) 


It is easily verified that Eqs. (13) and (11) are identical if the 
exponential in Eq. (13) is replaced by its infinite series expansion 
and like powers of (1/m) are collected. 

Since Eq. (13) is a closed expression for ¢(£), the convergence 
question is decided at once: 


¢(£) converges for all finite & = 
t/r, provided that m = Pe/P > 1. Thus no finite critical time 


exists. 
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The Minimum Weight Design of Wings for 
Flutter Conditions 


E.P. MacDonough 

Chief of Structures, Chance Vought Aircraft Division, United 
Aircraft Corporation, Dallas, Tex 

April 29,1953 


SUMMARY 


A method 
wing to obtain maximum values of primary bending and torsional frequencies 
for a given amount of structural weight With this method, the essential 
information for a preliminary flutter analysis can be provided early in the de- 
sign. This information includes mass distribution, primary bending, and 
torsional mode shapes and the relation between primary bending and tor 
It is necessary to determine what frequencies are needed 
The material requirements can then be de 


s presented for defining the best distribution of material in a 


sional frequencies 
to ensure freedom from flutter 
termined directly The procedure as presented is applicable only to the 
Primary bending and torsional modes 


DISCUSSION 


7 DESIGN OF THE PRIMARY STRUCTURE of a wing is the re- 
sult of compromises between the requirements for bending 


and torsional strength and bending and torsional stiffness. The 
methods for determining the optimum spanwise distribution of 
material to meet strength requirements are well known. It has 


generally been the practice to calculate the flutter speed that re- 
sulted when material was distributed according to the strength 
and aileron effectiveness With the increased 
speeds and lower thickness ratios that are being encountered, 
flutter requirements are becoming more critical and must there- 
fore be given more consideration in establishing the wing design. 


requirements. 
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Fic. 1. 


For a given wing plan form and weight distribution, the flutter 
speed is a function of the bending and torsional frequencies. It 
is, therefore, desirable to have a procedure by which the maxi- 
mum frequencies can be obtained for a given amount of structural 
weight. The procedure presented below is applicable only to 
the primary bending and torsional modes of the wing. 

The minimum strain energy for a given loading condition is 
obtained when the structural material is distributed in such a way 
that the energy per unit volume is constant throughout the struc- 
ture. For vibrating systems it can also be shown that the mini- 
mum weight of structure to attain a given frequency is approached 
when the energy per unit volume is constant under the loading 
associated with the primary mode. The condition of constant 
energy per unit volume exists when a constant stress is used. 
With a constant stress, the deflection shape of the primary bend- 
ing and torsional modes is a function of the effective depth of the 
wing primary structure. The acceleration at any point along 
the span can then be expressed as a function of the displacement. 
The spanwise distributions of acceleration and wing weight may 
be used to determine moment curves corresponding to the as- 
sumed primary mode shapes. The flange area and skin thick- 
ness can then be determined to resist the moments at any point 
along the span and provide the natural frequency that is neces- 
sary for the required flutter speed. The calculations for this 
work can be set up in a simple tabular form by dividing the wing 
into several chordwise strips. A detailed description of the pro- 
cedures is presented below. 


(a) Torsional Frequency 


Consider a box beam, as shown in Fig. 1, where spanwise sta- 


tions, y, are measured along the structural axis 


Flange area, An. = tc. 
Torsional shear stress, fp = 77/2At = T/2An.h, 
Angle of twist, @ = T//GJ = yil[(1/h.) + (1/c)], where y = 


shearing strain. 
At any station y, 
6 yf 1 ] 
= + dy 
Y 0 he Cc 
The torsional moment resulting from the inertia of a chord- 
wise strip, where the inertia is J, is 
AT = I,wd 
where w is the desired torsional frequency in radians per second. 
The torsional moment at any station is then 
l 
>, AT; 


1 Vv 


T = 


(Note that the torsional moment will be expressed as a function 


of shezring strain, y.) 
The required flange area at any station can then be determined 


directly from the expression 


An. = T. 


2Gh.y 


(b) Bending Frequency 
The bending stress, f, = M/h,An 
The radius of curvature of the deflected beam is R = h,/2e and 


the deflection at any station y is 
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y y "y y 1 7 
6 = I (2e/h,)dy | dy = 2e ( )as | dy 
0 0 J 0 0 h, 


The load resulting from the inertia of a chordwise strip, where 
the mass of the strip m = w/g, is 
Ap = mw*d 


where w is the desired bending frequency in radians per second. 


The shear at any station is 


and the moment at any station is 


a 
/ S dy 
/7y 


(Note that the bending moment will be expressed as a function 


M = 


of bending strain, e.) 
The required flange area at any station can then be determined 
directly from the expression 


An. = M/Ehe 


In order to make effective use of this procedure, methods 
must be available to estimate the primary bending and torsional 
frequencies necessary to prevent flutter. The procedure then 
provides a systematic method for studying the influence of vari- 
ous parameters on the primary frequencies of a wing and pro- 
vides the information needed to attain the required frequencies 


most efficiently. 


Gap Effect on Slender Wing-Body Interference 


H. Mirels 
Cornell University, Ithaca, N.Y. 
April 9, 1953 


Te LIFT EFFECTIVENESS of wings that are hinged, rather 

than solidly connected, to a cylindrical body was discussed 
in a previous note.' This arrangement introduces a gap, between 
the body and wing root, which becomes more pronounced as the 
wing is differentially reflected in respect to the body. The pres- 
ence of the gap reduces the lift of the wing-body combination. 
An estimate of this loss in lift has been published by Bleviss 
and Struble,? who replace the body by a reflecting wall and thus 
reduce the problem to one involving a wing and its mirror image. 
When the wing and body are at the same angle of attack, how- 
ever, the loss of lift can be found exactly within the limitations of 
slender-body theory. The solution is presented herein. Unless 
otherwise specified, the notation is the same as that used in refer 
ence 1. 

Consider the wing-body combination of Fig. 1. 
and body are at angle of attack @ with respect to a free stream of 
The body is cylindrical of radius a and the gap is of 
According to slender-body theory, 
two 


Both wing 


velocity U. 
constant width (yy, — a) 
the cross flow in each yz plane may be considered as a 
dimensional problem. The complex velocity in the yz plane can 
then be shown to equal 

—ial x 


al) = 


UW + 
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Fic. 1. Slender wing-body combination with gap 


Eq. (1) satisfies the following required boundary conditions 


(a) The velocity normal to the wing and body surface is w = 
-aU and v; = al’ sin 9, respectively. 

(b) The perturbation velocity potential is zero in the plane of 
the wing for a? < 6?= y? < y,?. 

(c) The complex velocity behaves like 1/¢* for large ¢. I 
Eq. (1), A’ and E’ are complete elliptic integrals of the first and 
second kind with modulus 
i i C¥a* — wi*) (9°? — a*) 

\ 2)? 
The other symbols are defined in either Fig. 1 or reference | 
The lift of the portion of the wing-body combination shown 


au 2/ Vo" 


in Fig. 1 (i.e., excluding forebody lift) is then 


Laa g 


2rqa 


Designating by Lea the lift for the no-gap case (y, = a but 4 
unchanged), 
| 2rqga = [yo — (a? yo) ]? 3 


The loss of lift due to the presence of the gap is 


Lea g Lea = | 2(E’ 


where 


yor + a (* “) | 
Vil Yo rT a yo a 


From Eqs. (4) and (5) it is seen that the lift is a function of & 


For k?< 1, Eq. (4) becomes 


only 


a 4 a 2\2 E - 
= (yo? + m7)] 1 4 21 Ve ; 4a? 
V7 V2? . V2 K 





LIFT RATIO — Laag /Lea,ere | 


Lee 


Eq. (4 
is used 
initial | 
large 


value 


If th 
no lift 


where 
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Eq. (4) is plotted as the solid line in Fig. 2. A logarithmic sca 
isused. Because of the logarithmic nature of A’ for small k, the 
initial loss of lift, for even the most minute gap sizes, is relatively 
large. Thus, for k = 10~, the lift is 87 per cent of the no-gap 
value 

If the body is considered to act as a reflecting surface, carrying 


no lift and inducing no upwash, Eqs. (4) and (6) become 


tai —2 (2) +40 | 


2 [ ky? ky? 
In (4/ko 2 4 In (4/k 
where the elliptic integrals are of modulus ko’ and 
ko’ = V1 — he? ( ‘ 
ko = (y; a)/(ye a) 


Eq. (7 corresponds to the case of an isolated slender wing at 
angle of attack a with a central semigap width of y,; — a and isa 
result presented in reference 2. It is noted that Eq. (7) has the 
same functional dependence on ky that Eqs. (4) and (6) have on 
k. These parameters differ, at most by a factor of 2. Namely, 


k/2Z<kjp<k (9 


If the effectiveness of the wing as a control surface is desired, it 
is necessary to solve for the case where the wing is at angle of 
attack a and the body is at zero angle of attack. The flow about 


such a configuration was presented in reference 1 for the no-gap 


? Lift on wing-body combination including gap effect. 


case. When the gap effect is included, the solution appears to 
involve untabulated functions, even for the slender-body case, 
and recourse to numerical methods is required. It was shown in 
reference 1, for the no-gap case, that the lift of this configuration, 
| ee 
| ene 


present, then 


is equal (within 4 per cent) to the geometric mean of Lg and 
If the same relationship is assumed to apply when a gap ts 


=) CE) 
— (10 
Wie) \ ha 
Eq. (10) depends on both ky and & and is plotted in Fig. 2 for the 
limiting cases ky) = k/2 and ky = k 
If there is no afterbody and the flight Mach Number is suffi 
ciently above 1, slender-body theory tends to overestimate the 
lift on the body. In this case it may be more fruitful to follow 
the approach used in reference 2 
come identical with those presented herein for ky = k 


The results of reference 2 be 


The presence of a boundary layer will tend to modify the re 
Viscous effects might pos 
sibly be incorporated into Eqs. (4), (7), and (10 
k and ky by the ratio of the estimated mass flow through the gap 
by the 

At any rate, with 


sults for the case of very small gaps 


by multiplying 
(based on viscous considerations) divided theoretical 
value of mass flow based on potential theory. 
increasing gap size, the equations become more realistic and a 


large loss of lift is indicated 
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The Propulsive Efficiency and Best Division of 
Power in an Augmented Jet Engine 


John Mitchell Stephenson 


Research Associate, Division of Engineering, Brown University, 
Providence, R.1. 


May 4, 1953 


(1) INTRODUCTION 


Fev ENGINE USED FOR PROPULSION consists of two basic 
elements: (a) a gas generator, which converts the chemical 
energy of the fuel burnt into the mechanical energy of a hot 
stream of gas; 
this energy into axial kinetic energy and, thence, into work done 
by forward thrust. 
without the augmenting effect of a power turbine driving a screw 
propeller or fan. 


and (b) a propulsive device, which converts 


The device consists of a jet nozzle, with or 


Controversy over the best type of thrust conversion for a par- 
ticular aircraft duty has touched most people connected with 
aeronautics in the past decade, and the subject is still unsettled. 
Even 5 years ago it was widely accepted that pure iet engines 
that the 
but that, for extremely long 


were only suitable for short-range military aircraft; 
turboprop gave a better range; 
routes, the old-style piston engine held the field. Today, exactly 
the opposite arrangement is often recommended. It is even less 
clear (except to their proponents) where the supersonic propeller 
and the by-pass engine fit in the economic scheme of aircraft op- 
eration. 

The extra weight, complexity, and cost of an engine with 
augmentation of the jet is only justified by its higher propulsive 
efficiency in the required region of operation. Therefore, one 
step toward a solution of the problem is an analysis of the pro- 
pulsive efficiency of engines with two separate sources of thrust. 

An excellent discussion of the ideal propulsive (‘‘Froude’’) 
efficiency of propellers and jets was given by A. V. Cleaver and 
D. H. Mallinson in Aircraft Engineering for June and September, 
1945. 
in the conversion process. 


But neither author considered the friction losses involved 


Then, in the issue of Aviation Week for November 15, 1948, 
there appeared an article on turboprops by Ivan H. Driggs which 
took the conversion efficiency into account. It was shown that 
there is an optimum division of power between propeller and jet, 
which depends not just on the forward speed but on a new ‘“‘num- 
ber’’—the ratio of the ram power to the available engine power. 
Thus it is not possible to say, for instance, that ‘pure jets are 
the best above 500 m.p.h.”” It depends on the altitude and the 
The 
analysis was based on the assumption that the propulsive effi- 


engine top temperature, as well as the propeller efficiency 


ciency of a propeller is constant. 
In the present note it is only assumed that the conversion 


that 
is, that the friction coefficients of all the components are inde- 


efficiencies of propellers, fans, and jet nozzles are constant 
pendent of scale. It then appears that for any value of Mr. 
Driggs’ number ¢, the optimum division of power between the 
two jet streams depends on the division of mass flow or “by-pass 
ratio.”” If this ratio is extremely large, as for a turboprop, the 
result is similar to that of Mr. Driggs. Curves of propulsive 
efficiency against ¢ are shown in Fig. 1 for by-pass ratios of 0, 1, 
10, and ©, assuming a jet efficiency of 0.95 and a fan efficiency 
(including the power turbine) of 0.80. 

Many results of interest can be calculated from these curves 
For example, the least efficiency of a supersonic propeller —below 
which the by-pass engine is better 
range of ¢. 
knots at the tropopause, a reduction of jet-engine top temperature 


is about 0.60 in the usual 
Another result is that for an aircraft flying at 500 


from 1,150° to 900°K. gives approximately the same improve- 
ment in propulsive efficiency as augmenting with a by-pass ratio 


of 1. 
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(2) THE PROPULSIVE EFFICIENC\ 

Let G, be the power generated, per unit mass flow, in an isen 
tropic expansion from the pressure and temperature at the gas 
Let G, 
If the forward speed is 


generator exit to atmospheric pressure. be the net thrust 


power per unit air-mass flow. v, the ram 
power of the entering air is (1/2)v*. Then the propulsive efi- 
ciency is defined as the ratio of the thrust power to the power 


added to the air stream by the engine 


np = Gz/|Gr — (1/2)0?] 


For comparison between engine types, it is sometimes better 
Cp(1/2)v? 
ternal drag coefficient of the engine installation, based on the up 


to write the numerator as [G, |, Where Cy, is the ex 
stream entry cross-sectional area of the flow stream tube. But 
the drag will not be considered here 

Consider a system where the gas from the generator is exhausted 
through a propulsive jet, at a velocity v;, but a fraction a@ of the 
available power is tapped off to drive a fan or propeller. The 
final exit velocity from the fan is v7. Then, if the ratio of the 
mass flows through the fan and gas generator is \—-the “by-pass 


ratio’’—we have 


Gz, = ov; — v) 4 


Arve — V) 
Hence, writing (1/2)v?2/G, as ¢, Eq. (1) becomes 
g 1 
1+A 2 


np = [2¢/(1 — $)] [(v;/v) + (dv,/v) - 


Now we define the conversion efficiencies of the jet nozzle (»;) and 


of the fan or propeller (7,) as the ratio of the axial kinetic energy 


produced to the energy input. Thus, 
nj = (1/2)0;7/(1 — a)G, = $(2;/v)2/(1 — a 
ny = (1/2)A(v;? — v?)/aG, = AP[(v7/v)? — 1]/a 


Substituting in Eq. (2), we find the propulsive efficiency to be 
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2 (1 — a) Aany 
io i 4. alas 4 —Y — (1 + a) (3) 
l|-*©¢ g \ ti) 
There are two extreme cases of this propulsive system: firstly, 
the pure jet engine, when a = d = 0), and 
io [2 | ¢)) (v Ni ob P 1) (4) 
Secondly, the turboprop, where A is extremely large, and 
» | 
<> l — a) “nf " 
= ( a + —— == (5) 
1-—@¢ o 26 


(3) THe Best DIvIsION OF POWER 


We can now find, for any value of ¢, the value of a which gives 


the highest propulsive efficiency. For 0n,/0a@ = 0 when 


1 — (¢n;/n;*) : 
a= (6) 
1 + (;/Any) 
This shows that there is no gain in propulsive efficiency at all by 
augmenting when @ > ny?/n;. 
For the turboprop, A is extremely large, and Eq. (6) becomes 
(1 — a) = ¢n;/n;/? (7) 
and substituting in Eq. (5), the best efficiency of a turboprop is 


9 
= <o ns 4 n) = (8) 


dial 9 9 
1—@¢ \2¢@ 2n; 


(4) PROPELLER EFFICIENCY 
The efficiency of a screw propeller is usually quoted in the form 
of thrust power divided by the torque—.e., in such a way as to in- 
clude the propulsive efficiency, as well as the friction and eddy 
losses. Thus, including the gear and power turbine losses as well, 
we define 
7 = Av( ve — v)/a G, 
and Eq. (8) becomes 


26 I l—a )n an . 
Np = = + | (8a) 
l—@¢ rs) 26 


Then if we assume that 7 is constant, differentiation of Eq. (8a) 


yields the best division of power as 
(1 — a) = gnj/n? (Va) 
Substituting in Eq. (8a), we have the best efficiency of a turbo- 


prop as 


2 n n 
1p = -— — } (8a) 
l - 26 2n 





However the assumption on which this latter analysis rests is 
certainly not true for a set of fans and may be doubted for a wide 
range of propellers. 

I have recently learned that P. C. Gupta, of Bristol, England, 
has arrived at a more complete result. In particular, he takes 
into account the intake efficiency n; and the outlet efficiency mo 
of the fan duct, so that the fan efficiency is written 
ny = (Ad a)-[(v, v)? no — nl 
The expressions for the propulsive efficiencies of the pure jet and 
the turboprop are, of course, unaltered, but the best division of 
power in the by-pass engine takes a form slightly different 
from Eq. 6, 

1 — (pn;n;/n;?n0) 

ie (n; Anyno) 

This refinement of the analysis is valuable, for it shows that the 
by-pass engine is extremely sensitive to the quality of the fan 
duct 
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On the Role of Three-Dimensional 
Disturbances in the Stability of Supersonic 
Boundary Layers* 


D. W. Dunn and C. C. Lin 

Department of Mathematics, Massachusetts Institute of Technology, 
Cambridge, Mass. 

May 4, 1953 


| i A PREVIOUS NOTE,!' it has been pointed out that, although 
the equations for three-dimensional disturbances in a com- 
pressible fluid transform formally to two-dimensional ones, it 
cannot be concluded from the properties of the transformation 
that three-dimensional disturbances necessarily have a higher 
critical Reynolds Number of instability (as in the case of an in- 
compressible fluid?). Calculations have now been carried out to 
demonstrate this point with numerical results and to show the 
exact nature of this effect of three-dimensional disturbances 


The case considered is a boundary layer in a gas with a free- 
stream Mach Number 1/ = 1.6. For convenience of calcula- 
tion, the viscosity coefficient is taken as proportional to the abso 
lute temperature. Estimation of errors involved shows that 
this approximation would not influence the main results to any 
appreciable extent 

Present calculations show that two-dimensional disturbances 
will be completely stabilized* at a wall-to-free-stream temper- 
ature ratio of 1.073. The critical Reynolds Number of wavy 
three-dimensional disturbances is obtained as a function of the 
angle 6, which the direction of propagation of a given wave makes 
with the free stream. For convenience of application, the Reyn- 
olds Number is represented in terms of the downstream distance 
in Fig. 1. 

It is seen that the oblique disturbances have a critical Reyn- 
olds Number of the order of 10°, leading to the conjecture that 
the transition Reynolds Number would be of the order of 107 
108, with turbulence level prevailing under ordinary experimental 


conditions. Thus, by such cooling processes, one can materially 


* Work done under Contract N5ori-07872 of the Office of Naval Research 
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delay the transition process but cannot achieve complete stabiliza- 
tion. 

Studies of small disturbances also lead to the conjecture that 
surface roughness is not always so effective in producing transi- 
tion in the boundary layer in the supersonic case as in the sub- 
The disturbance produced by a stationary rough- 
ness (including surface roughness and other stationary obstacles ) 


sonic case. 


is a supersonic disturbance‘ and has the nature of a Mach wave 
outside the boundary layer. The energy of such disturbances is 
easily carried away from the boundary layer in the nature of 


sound waves. 
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Apparent Recovery Factors on the Leading Edge 
of a Flat Plate in Supersonic Laminar Flow 


Irving Korobkin 
U.S. Naval Ordnance Laboratory, White Oak, Silver Spring, Md. 
April 16, 1953 


W 7 IND-TUNNEL MEASUREMENTS on a flat plate with a super 
sonic laminar boundary layer consistently result in the de 
termination of recovery factors higher than the theoretically pre 
dicted value, Pr®*, and the experimental values obtained on 
cones.'~> The leading edge of the flat plate model is usually 
beveled away from the measuring surface forming a wedge, and 
it will be shown that the high recovery temperature measurements 
are the result of steady-state heat conduction in the wedge 

An approximate calculation of recovery temperature can be 
made for a thin flat plate with its measuring surface oriented 
parallel to the flow and with a solid metallic wedge of small angle 
which is insulated from the rest of the model as the upstream sec 
tion. In Fig. 1 the origin of the orthogonal coordinate svstem 
used in the calculation is located at the tip of the model with the 
x-axis lying along the wedge-angle bisector. 

The measuring surface of the plate (Surface 2) has a local free 
stream Mach Number equal to the Mach Number upstream 
from the model, but the other surface (Surface 1) has a lower 
free-stream Mach Number because of the presence of a shock 
wave produced by the flow deflection. It will be assumed that 
each surface has a different constant recovery factor equal to 
Pr®-® based on conditions outside the boundary layer adjacent to 
each surface, with the corresponding recovery temperatures being 
T, and 7». 

The contour of the wedge can be specified by 

y = Kx (1) 


The expression for the local heat-transfer coefficient on an iso 
thermal flat plate with laminar boundary layer is given in refer- 
ence 6 as 


Nu = 0.29Pr'/*Re®-5 (2) 


where Nu is Nusselt Number, Pr is Prandtl Number, and Re is 


Reynolds Number. Therefore on each surface 
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where “1 and & are constants of proportionality and h is the loca 
heat-transfer coefficient. Assuming the plate to be thin and th 
wedge angle small, the variation in plate temperature in the 4 
direction may be neglected. The following differential equatior 


then describes the steady-state heat-transfer process in th 


d (,, dT ™ = ‘- 
de N29 Be Ge) = (Te — Td + tel T 2 — 1 { 


wedge: 


where &,, is the thermal conductivity of the wedge and 7°, is the 
wedge temperature. 


After substituting y, 4, and he, the following expression is ob 


d ( <<) 
v cx0-8 (7 
dx dx 


where c = (A + B)/2k,A and tr = (AT, + B7;)/(A +B 


Eq. (5) is now reduced to the same form as used in reference 7 


tained: 


for heat transfer to a wedge-shaped fin in natural convection, th 
resulting solution being 


: r= Esi\(4u Vc) + F¥i(4x 44 Vc) 6 


where £ and F are integration constants and Jy and }¥» are Bessel 
functions of first and second kinds of order zero. When x equals 


zero, then 


Vo(4x (ti Vc) = ‘ 


Since the temperature is finite along the wedge, F must be zero 
and, therefore, 


1 r= EJS (4x41 Ve 8 
The insulation at x = 1 means that 
dT ./dx|\ r=; = 0 4 
Since 
dT’, (dx == —Eiv cl ‘I(4L t4/ ¢ 0 
then 
E=0 lo = 7 11 


The above analysis indicates that the entire wedge, under the 
assumed conditions, is at a constant temperature that justifies 
the use of Eq. (2) for the local heat-transfer coefficient on an iso 
thermal flat plate 
The apparent recovery factor measured on the wedge is given 
by 
yr =(r — Toao2)/(1o Tas2 12 


where 7) is stagnation temperature and 7’. is local free-strean 


temperature along the measuring surface. 
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TABLE | 


Wedge Angle 
(20 of Fig. 1 


Apparent Recovery Factor 


0 0.869 
5° 0.874 
10° 0.880 
15° 0. 887 
20 0. 897 


Table 1 is based on the foregoing analysis and is computed at a 
Mach Number of 2.4 for a supersonic wind tunnel with a 1- 
atmosphere supply pressure and a 300°K. supply temperature, 
where the theoretical value of recovery factor based on free 
stream conditions is 0.869. 

It can be seen from Table 1 that the wedge is at a temperature 
higher than the value obtained from Pr*-® when calculated for 
the measuring surface of the flat plate. Consequently, the 
boundary layer in contact with the wedge is warmer than under 
zero heat-transfer conditions. The effects of this high temper- 
ature are felt downstream along the main part of the model, sinc« 
the plate now attains thermal equilibrium with a boundary layer 
of higher temperature resulting in the measurement of erroneous 
recovery factors. This phenomenon must be considered when 
dealing with all flat-plate wind-tunnel models with laminar 
boundary layers. Though the above analysis assumes a thin 
plate with a solid metallic wedge of small angle which is insulated 
from the rest of the model, the same effects will be present in flat 


plates with different details of construction. 
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Effect of Gas Radiation in the Boundary Layer 
on Aerodynamic Heat Transfer 


ohn W. Smith 
Aerodynamicist, Bell Aircraft Corporation, Buffalo, N.Y 


May 7, 1953 


HE USUAL CALCULATION of skin temperatures in very-high- 
4 onke supersonic flight involves the balancing of two major 
heat conducted into the skin from the boundary layer, 
The latter keeps 


terms: 
ind heat radiated by the skin into outer space 
skin temperatures at values less than adiabatic wall temperature 
and also lowers temperatures in the interior of the boundary 
layer. Nevertheless, at critical flight conditions, temperatures 
many times higher than skin temperature can occur in the bound- 
ary layer. This raises the question whether radiant emission and 
ibsorption by the gas molecules in the boundary layer, hitherto 
neglected, can appreciably influence the aerodynamic heating. 
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Skin temperatures could be affected in two ways: (1) radiant 


heating or cooling of gas elements in the boundary layer can alter 
the temperature profile, causing a change in the conduction heat- 
transfer coefficient; and (2) radiant emission from gas molecules 
causes an additional quantity of heat to be transmitted to the 
skin. These effects will be briefly considered here. 

The chief constituents of atmospheric air are nitrogen and 
oxygen molecules, which are practically inert in the infrared fre 
quency range. The greatest contribution to radiative heat trans 
fer in atmospheric air is due to relatively minor constituents, 
such as water vapor and carbon dioxide. The chemical composi- 
tion of the boundary layer, on the other hand, may differ appre 
ciably from atmospheric air because of the high temperatures, 
with the result that such gases as nitric oxide may contribute ap 
preciably to the radiative heat transfer. For the present it will 
be assumed that there exists a certain density concentration 
of radiating gases in the boundary layer and that the emissivities 
of these gases are comparable to that of water vapor. Let the 
emissivity ¢ of a thin slab of gas be defined as the total energy 
emitted across a unit area of one face in unit time divided by the 
black body intensity o7‘*(o = Stefan-Boltzmann constant; 7 = 
gas temperature). For slabs of small optical thickness npy (p = 
total density; y = thickness), the emissivity is expressible in the 
differential form, 

de = o(T)np dy (1 
where ¢(7) is independent of y—i.e., reabsorption can be neg 
lected. Consider now the radiative energy balance for an in- 
finitesimal slab of gas in the boundary layer with faces parallel to 
the solid surface. The incident radiant intensity transmitted 
across the lower surface of the slab is approximately o7 
(7, = skin temperature). Since absorptivity and emissivity are 
related by Kirchhoff's law, one obtains 


de de 


=0¢ (Tw*t — BT*) = of( T)np(T 04 sT* (2) 
dt dy 


where de/dt is the rate of change of energy in the slab per unit 


8 =2. The instantaneous 


volume and the geometrical factor 
radiative cooling is therefore given by 
dT 1 de od( 7 )n 

= - (Te* BT* (3) 
dt pc; dt Cy 


For the numerical example cited below, this can assume negative 
values of about 10*°R. per sec. 

In order to calculate the effect of radiative cooling on the tem 
perature profile, it is necessary to solve the boundary-layer energy 
equation subject to an additional term expressed by Eq. (2 
The modified energy equation (in the ordinary notation) is 


Ry) acts op od oT 
pu (cpT) + pv (Cc, T) =u t k T 
Ox oy Ox oy oy 
(oy de 
Mt t ) 
ov dt t 
Let 7 = 7) + 6, where 7» satisfies the energy equation neglecting 
the radiation term. Then @ satisfies the equation 
re) ro) 2) O80 de 
pu > pf) tT pv - (¢ pf) = k t 5 
Ox oy oy ov dl 


Suppose, for the moment, that p, u, and 7’ depend only on the 
“similarity variable” » = y/+V/ x, and let Eq. (5) be transformed 
from the (x, y) plane to the (x, 7) plane. One can eliminate pi 


by use of the continuity equation 


= 1 *y 
pv = 4 f pudy = ( { pu dn n pul ) (6) 
ox J 2V x \v0 


Eq. (5) then becomes 








580 JOURNAL OF THE 


1 f. F dé ld k dé 4 de 
es pu dn a 7 = ¥ 7) 
2x J0 dn x dn dn dt ( 
Because of the radiation term, x does not drop out, and the as- 
sumption of similar profiles cannot be maintained. 


This result greatly complicates the general solution of Eq. (5). 
The approach adopted here assumes that the convection terms are 
small in comparison to the conduction term for points sufficiently 
far downstream. This leads to a simple solution for @, which in 
turn is used to justify the given assumption a posteriori. Neglect- 
ing convection terms, Eq. (5) becomes an ordinary differential 
equation, 


0 (8) 


(d/dy) [k(d0/dy)] + (de/dt) 


Asa first approximation, one can set 7 = To in the radiation term 
to obtain the integral, 


Tr f° dy'k (Te) f dy"$(T)e(8To* — Tw*) + 
. y . 
Kf k-\T>) dy + Kz (9) 


The integration constants A; and K» are determined by the con- 
dition that @ vanish at the boundaries y = 0 and y = 6 (6 = 
boundary-layer thickness). 

We are now in a position to calculate the change in heat con- 
viz., 


aie E (2 o7e) -x 
: oy oy a = ' (10) 


duction 


On the Flow at the Rear of a Two-Dimensional 
Supersonic Airfoil with Thickness* 


Ronald F. Probstein 


Assistant Professor, Aeronautical Engineering Department, 
Princeton University, Princeton, N.J. 


May 4, 1953 


AN TREATMENT OF THE PROBLEM Of the flow deflection immedi- 

ately rearward of a two-dimensional supersonic airfoil was 
given by Kahane and Lees! in 1948. Their analysis was carried 
out on the basis of oblique shock and Prandtl-Meyer wave 
theory, utilizing the fact that the dynamic and static pressure 
change through isentropic and shock waves can be expressed as 


* This note is an outgrowth of work supported by the Office of Naval Re 
search and Office of Scientific Research, ARDC, U.S. Air Force, under Con 
tract No. N6-onr-270 


(—3y? + 17y 


12 12 


AERONAUTICAL 


SCIENCES—AUGUST, 1953 


as well as the radiative heat transfer from the boundary layer to 
the skin, the latter being expressible as an effective reduction jn 
the emissivity of the skin by an amount Ae. According to Eq, (1), 


this becomes 
Pens Eg 
Ae = —n f, (T)p T., dy (11 


A numerical calculation was carried out for the following condi- 
tion: flight Mach Number = 20; altitude = 150,000 ft., ratio 
of wall to free-stream temperature = 4; m (active gas concentra- 
tion) = 1/10; and x (distance from the leading edge) = 10 ft. 
Laminar 7)-profiles were taken from Van Driest,' and ¢(T 
was taken from the data of Hottel and Egbert? for water-vapor 


emissivities. The results are, summarily, 


(0 T) max. = —2.2 X 10 3 
(Aq/q) = =] X 10-° 
Ae = —5.3 X 1074 


indicating effectively no change in any of the given quantities. 
It should be noted that the net effect, though small, tends to lower 
the equilibrium skin temperature. 


REFERENCES 
1 Van Driest, E. R., Investigation of the Laminar Boundary Layer in Com- 


pressible Fluids Using the Crocco Method, NACA TN No. 2597, January, 1951 


2 Hottel, H. C., and Egbert, R. B., Radiant Heat Transmission from Water 
Vapor, Trans. Am. Inst. Chem. Eng., Vol. 38, No. 3, June, 1941. 


a power series in the flow deflection angle. Since, to the writer’s 
knowledge, this is the only place where the dynamic pressure 
coefficients have been published, it was thought worth while to 
point out that one of the coefficients in the expansion for an isen- 
tropic wave (also appearing in the shock expansion) appears to 
be incorrect. 


The authors expressed the change in dynamic pressure through 
an isentropic wave as 
Aq = a(+ G6 + G.6? + G:63 ++...) 
and through a shock wave as 
Aq = q(+ G60 + G.6? + (G; — H)e® +...) 
where g = pw?/2(p = density, w = flow velocity); @ is the angle 


of detlection, positive in a counterclockwise sense; and q is the 


dynamic pressure ahead of the wave. The corrected value for 


G; should be 


— 36) (21 — 4y) 14 4 


Ms + —_—— M* — 8M* + - - 





where ./ is the flow Mach Number ahead of the wave. 
Determination of the flow angle immediately rearward of a 
flat-plate airfoil at angle of attack and for a thick airfoil at angle 
of attack, where it was assumed that the configuration was such 
that four oblique shock waves were attached to the nose and tail 


(M2 — 1)'/72 


» « o 


of the airfoil (see Fig. 1), was carried out. Fortunately, the cal- 
culations did not involve the use of the incorrect coefficient G 
directly. In spite of this, however, for the thick airfoil shown 
in Fig. 1, the analytic formulation of the flow deflection angle 8 to 


first order appears to be incorrect. 
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The thick airfoil 


Fic. 1. 


Recalculation of 8 gives 


D 
B= 3c [(83 + a)? + (6, — a)? — (& — a)® — (82 + @)8] 4 
«Vl 
E _— 
2C [( 5s t a)! + (6; a)! (84 a)* - (be | a)*] + 
“(1 
H/2|(d2 + a)* — (&: — a)4] — 
[((dD/dé) + DG,}) 
[(63 + a)4 — (6; - a)*] 


2C, 


where D, Ci, and E are the coefficients of @ in the following ex- 
pansion for the static pressure change through an oblique shock 
wave: 

3 —-D)e?+(CG,— Ej*+...] 


Ap = gq, [+Ci0 + C.6? + (6 


The required coefficients for the determination of 8 are functions 
only of the free-stream Mach Number, 1/;, and the ratio of spe- 
Their values (other than £) are given in reference 
The computation 


cific heats, y. 
FE) may be found in reference 2. 
E is easily carried out by use of the results of 


1, while (C, 
of C, and, 
Miles. 


hence, 
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Interaction of a Plane Flame Front with a Plane 
Sinusoidal Shear Wave* 


George H. Markstein 


Head, Combustion Section, 
Inc., Buffalo, N.Y 
May 6, 1953 


Cornell Aeronautical Laboratory, 


_— INTERACTION OF A WEAK PLANE DISTURBANCE with 
plane shock wave has been studied recently by several in 
vestigators.'~? These linearized treatments were carried out 
as a first step toward the solution of problems of interaction 
between turbulent flow and shock waves. An analogous ap- 
proach to the problem of interaction between turbulent flow and 


flame fronts appears to be of interest, particularly in connection 


* This work was conducted under the auspices of Project SQUID, jointly 
sponsored by ONR, ARDC, and OOR under Contract N6-ori-105, T.O.1II 
Subcontract 1 


FORUM 581 


with observations and analyses of ‘“flame-generated turbu- 
lence.’’* § 

Foliuwing the procedure of references 1 and 3, the general case 
of interaction between a combustion wave with an oblique sinu- 
soidal shear wave can be transformed into a steady-flow problem 


tangential to the flame front that 


by superposing a velocity | 
makes the resultant flow velocities parallel to the shear waves. 
This is shown in Fig. 1, which is analogous to Fig. 2 of reference 
1. The case of the combustion wave differs from that of the 
shock wave, however, in the following respects: 

(1) Since the burning velocity S, is extremely small compared 
to the speed of sound (the Mach Number is of the order of 10~8 
for hydrocarbon-air mixtures), compressibility effects can be 
neglected, except for almost glancing intersection of the waves 
(y and y nearly 90°), where the fictitious sweep velocity V and 
the resultant velocities W, and W» may approach and surpass 
sonic values. 

(2) For subsonic resultant velocity W,, a perturbation poten- 
tial flow must be considered upstream, as well as downstream, of 
the combustion wave. 

(3) While plane shock waves are always stable, combustion 
waves may be unstable,® at least for wave lengths of the disturb- 
ance larger than a critical wave length \,. This circumstance 
obviously reduces the range of validity of linearized analysis. 

In the present note, only the case of normal intersection of 
shear wave and combustion wave, ¢ = ¥ = 0, will be discussed 
in detail. For this case, the linearized equations of incompres- 
sible and nonviscous flow are satisfied by the perturbations 


Uy = (Ag +. ekY (1) 
pias .= rete = - fA 
upstream (x < 0), and 
= k 
uy/Sy = (Are 74 As)e*¥ (2 
. F = © ki +iy) «) 
—1 Up ae => Po Pura = 1o¢ ore, | 
downstream (x > 0) of the flame front. Here Ao, As are the 


As» those of the 
S,, and k 


amplitudes of the respective shear flows and A), 
potential flows, all referred to the burning velocity .- 
27/X is the wave number. 

The problem consists of determining A,, Ao, 
from the boundary conditions at the flame front 
displacement xp of the latter, which is caused by the flow dis 


A; for given Al 
The sinusoiday 


turbance, can be written 
k , 
kxp = Aye’ , (3) 


with a further unknown amplitude factor Ay. The simplest 
boundary conditions would be obtained by assuming constant 
density ratio p,/p» = € and constant burning velocity S,. How 
ever, these assumptions imply instability of the flame front for 
all wave numbers‘ and would thus invalidate the linearized anal 


ysis completely. While it seems reasonable to retain the as 
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sumption p.u/pp) = 
and chemical kinetics will normally cause perturbations of the 
burning velocity. Because of our present ignorance about the 
details of these effects, they are lumped into a phenomenological 
parameter y. The first-order perturbation of burning velocity is 
thus introduced, as in previous work, in the form 


Su/Su = wL/R = —pl(d*xp/dy?) = oAge™ (4) 
with 
o = kpl = 2aul/r (5) 


where R is the local radius of curvature of the distorted flame 
front and L is a length of the order of the thickness of the com- 
bustion zone. 

With these assumptions, the boundary conditions for velocity 
components and pressure lead to the equations 


Ay =——_ Ag = —Ay 
A, + As —eo Ay = 0 (6) 
hi Me —~(e-14,= 0 . 


| 
oS 
—) 


Ay — As + 2a(« = 1)Ag ed 
The determinant of the coefficients of this system vanishes for 
o = (e — 1)/2e, a value corresponding to the wave length of neu- 
tral stability 

(7) 
The linearized analysis breaks down for \ > X,, since for these 
A lower limit is ob- 


Ae = 4rpLe/(e — 1) 


wave lengths the flame front is unstable.® 
viously also imposed on \ by the finite width of the combustion 
zone. 
Substituting Eqs. (5) and (7), one arrives at the following solu- 
tions of Eqs. (6): 
A, = Ao(d/de + em! — 1)/(1 — d/de) 
Az = Ao(X/re — €72 + 1)/(1 — A/de) 
Ale! = RRA = Xr) 
2A o(A/Ac)/[(e — 1) (1 — A/A-)] 


Ay = 


| (8) 
A, = 
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These results are plotted in Fig. 2 in the form of amplitudes of 
the various perturbations relative to that of the upstream shear 
flow. 

Even though no detailed conclusions regarding the much more 
complicated case of flame-generated turbulence can be drawn from 
this specialized result, certain features would seem to be of general 
the flame 


significance. Thus, amplification of shear flow by 


occurs only for 


X/Ne > (1/2) (e! + 1) 





although the flow perturbations immediately downstream of the 
flame front always exceed the upstream shear flow in ampli 
tude. This distinction between shear flow and local potential] 
flow that decays to extremely small amplitude within a few wave 
lengths up- and downstream of the flame would appear to be 
of equal importance in the case of actual turbulence. 

Another prominent feature of the Eqs. (8) is the ‘‘resonance”’ 
that occurs when A approaches A,._ This effect will, in reality, be 
suppressed to some extent by the leveling influence of nonlinearity 
on amplitude. Moreover, the amplitudes tend toward infinity 
only for normal intersection; in the oblique case one finds that the 
coefficient determinant vanishes only for a complex value of <, 
and thus no true resonance occurs for real values of . In spite 
of these factors, one would still expect a different response of 
flames to turbulent approach flow, depending on whether \, js 
large compared with the scale of turbulence or not. This expec- 
tation seems indeed to be borne out by observations on lean and 
rich turbulent butane flames.’ 


A fuller account of this work will be published later. 
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Errata—‘‘An Experimental Study of Natural 
Vibrations of Cantilevered Triangular Plates” 


P.N. Gustafson, W. F. Stokey, and C. F. Zorowski 
Carnegie Institute of Technology, Pittsburgh 
May 20, 1953 


, pw OMISSIONS OF SYMBOLS appeared in a paper by P. N 
Gustafson, W. F. Stokey, and C. F. Zorowski.! These are 
corrected as follows: 

(1) The last term on the right-hand side of Eq. (18), page 336, 
should read 


1 E’ph’1 — p*) 
a! Ep’h*1 — p’*) 
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(2) Eq. (20) on page 336 should read 
fr LBL Ep (1 — 2) a 
E ath E p’ (1 — pw”) ‘ 
REFERENCE 
1 Gustafson, P. N., Stokey, W. F., and Zorowski, C. F., An Experimental 


Study of Natural Vibrations of Cantilevered Triangular Plates, Journal of the 
Aeronautical Sciences, Vol. 20, No. 5, pp. 331-337, May, 1953. 


Transformation Between Compressible and 
incompressible Boundary-Layer Equations 


Nguyen Van Le 
AiResearch Manufacturing Company, Los Angeles 


April 10, 1953 


I IS PROPOSED TO SHOW that the boundary-layer equation of 
compressible flow can be reduced to that of incompressible 
flow. Such work was initiated by Stewartson! and by Rott and 
Crabtree.2 In the following some of the restrictions imposed by 
references 1 and 2 will be removed, and it will be shown that the 
transformation from compressible boundary layer to incompress- 
ible boundary layer can be applied to the laminar, as well as 


turbulent, case. A direct method will be used for this purpose. 


DERIVATION OF THE COMPRESSIBLE BOUNDARY-LAYER EQUATION 

The following assumptions are made (see Fig. 1): (1) two- 
dimensional flow; (2) Prandtl Number is taken equal to 1 (Pr = 
1); (3) in the case of turbulent boundary layer, the time average 
value of the density, velocity, . . . etc., are considered; (4) the 
usual boundary layer assumption of 0p/Oy = 0; and (5) no heat 
is transferred across the boundary-layer wall. 

In the control volume ABCD, application of the momentum 


theorem gives 


d é d ad dp 
pur dy — U, pudy = — — 7 (1) 
dx Jo dx Jo dx 
Outside the boundary layer there is isentropic flow 
T ds = 0 = dh — (1/dp,)dp, .dh = (1/p:)dp; (2) 
From the energy equation we have 
dh, + U, dU; = 0 (3) 
Combining Eqs. (2) and (3) yields 
(1/p1) (dp, /dx) = —U, dU; (4) 


Since 0p/dy = 0, p = pi, and dp,/dx = dp/dx, relation (4) gives, 
then, 


(1/p:) dp, = U, dU,/dx 


which can be written 


. apy 
6 = 
dx dx 


d ” - " - . 
_ aUerdy + Uy pi; dy (5) 
dx Jo dx Jo 


With the Eq. (5), the momentum balance from Eq. (1) yields 


d “6 d ~6 
(pi I y,? - pu?) dy - Ui (pi I A — pu) dy = TT) 
dx Jo dx Jo 


(6) 


d ‘ _ ad ; 
aU? + U; pi U, 
; dx 


The momentum and displacement thicknesses are now defined 


so that the energy equation can be taken into consideration in the 


boundary layer 
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With Pr = 1, satisfying the energy equation in the boundary layer 


has, as consequence, 


(8) 


Upon introducing the stagnation temperature 7) as a reference, 


Eq. (8) becomes 
u2 
— 1] 1 -—- — 
U? 


T , 4 (7 
T; fi 
In addition the assumption 0p/d0y = 0 shows that in the bound- 


ary layer the density and temperature are related by 


(9) 


pT = pi T; (10) 


With relations (9) and (10), one can evaluate the different terms 


of Eq. (6) 
. y ; To . u\p 
(mU, — pu) dy = pli = l . dy + 
Ti Jo Ui} pr 
u(? i) ‘(1 >) up fica iase a 
as T, A U,) U, ——" = ali 
K -1  f{K -1 
1 + 9 M,? 6) + pill 1 9 M0; ( l l ) 


5 F) 
[ (p,U,;2 — pu?) dy = Ui(oi.Ui — pu) dy + 
J0 0 


6 2 6 
p u u 
( a ~) dy = 
0 pi Ui U? 0 


Ui, X 
(p.U; — pu) dy + pr Ui? 


0 


Pi U;? 


(12) 


The momentum equation (6) then becomes 


d d : 
nh + = (mU; — pu) dy — 
dx dx Jo 


é 
_ @ ; 
U; (mU, — pu) dy = ro 
dx 0 
or after integration by parts 
d i dl . : ; 
pr Ui, + (mU, — pu) dy = ro 
dx dx Jo 


and, with Eq. (11), we obtain the compressible boundary-layer 


equation, 
d aU 
aU, + Ui ‘xX 
dx dx 
K -1 K -1 
| 9 M,70; + 4 (: + = us) | = 7 
or 
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dU, 1 dp, dU, 
2 U,"6 U 
dx 7 pi dx 77 dx 


K-1 — 0 
(2 + = us) A+ (: + r= ut)a| = : (13) 
2 é 1 


But for isentropic flow outside the boundary layer, combination 
of the isentropic relation and the energy equation gives 

— 1 dp, dl A 

sd 


= - UV, M;? 
pi dx dx 


Ui 


(14) 


Putting Eq. (14) in Eq. (13) we obtain the final equation for com- 
pressible boundary layer, 
do; 
dx 
0, dU, 
U, dx 


a2 (15) 


As derived, the equation does not make any assumption on the 
viscosity in the boundary layer. The time average values of the 
velocity and density are considered for the case of turbulent 
boundary layer, although there may be possibility of treating 
their time components as in reference 3. Until measurements 
show the large contribution of the time components, the previous 
equation can be applied to laminar, as well as turbulent, boundary 
layer because of the nature of its derivation. For M, = 0, it re- 
duces to the well-known incompressible case, 


dé, + (2 4 *) 4 dU; = To 


dx a; U; dx aU? (16) 


TRANSFORMATION FORMULAS 


A set of transformation formulas modified from reference 1 
is now adopted as follows: 
, a - 
Y = dy 
0 Po 


x \ (3K —1)/(K—1) 
=] 
0 (* 


= aoM, = (do ‘ay) Ui 
(ao /a,)u 
where 


sound velocity 7) 
sound velocity 7; 
Cp/Co 


a= 
q= 


| ae 


Define further the following integrals: 


5: 


A 


which can be interpreted as the displacement and momentum 
thicknesses of an incompressible flow where the characteristics 
are as shown in Fig. 2. Using the transformation formulas, one 
can show that 
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Fic.. 2. 


(pi/po) (a1/ao) (61) = 52 


which gives, further, 


(“)’ d (Vm) = ye 
ao/ dz Se aa 


a, 4 7 d J 4 
— Fh a Vs (20) 
ao dz 


Upon introducing Eqs. (19) and (20) in the Eq. (15) for com- 
pressible flow, one obtains 


a,\id ay _ av; T0 
Vi702 + 2Vi = 
ao/ dz do dz pi 


d0. be (4K 
+ (2+ 
dz Oo 


(ao/a,) 
If a modified wall shear stress is now defined 


r = 1 (ao/a) 4K -DMK—-V) 


Oe dv, To 2)/(K—1) 


(21) 
Vi dz 


po Vi? 


(22) 


the value of 7» being obtained from reference 4, for instance, 


Eq. (21) is reduced to 


d0» ri (2 ‘ *) A dV, J 
dz * wt é& wt 


which is exactly the equation of an incompressible boundary layer. 


(23) 


Thus, in order to obtain the solution for a compressible boundary 
layer, turbulent as well as laminar, one chooses a temperature 
profile from experimental data just in the same way as for the 
velocity profiles. Then one carries out the transformation pre- 
scribed in Eqs. (17) and (22) and solves the incompressible bound- 
ary-layer equation (23), which is more familiar. 

For comparative purpose it is reminded that the Polhausen 
method generalized to the compressible boundary layer required 
a solution for a system of two nonlinear differential equations. 
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